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FOREWORD 


This  report,  "A  Method  for  Predicting  Three-Degree-of-t reedom  Store 
cep.tr<ation  Tra3ectories  at  Speeds  up  to  the  Critical  Speed,"  describes 
* i.wkbined  theoretical-experimental  program  directed  toward  developing  a 
computer  program  for  predicting  the  trajectory  of  an  external  store  dropped 
: com  an  aircraft  of  the  fighter-bomber  type  at  speeds  up  to  the  critical 
sceea.  The  work  was  carried  out  by  Nielsen  Engineering  & Research,  Inp., 
8aO  Maude  Avenue,  Mountain  View,  California  94040,  under  Contract  No. 

F3361 5-69-C-1337 . The  contract  was  initiated  under  Project  8219,  Task 
811902,  of  the  Air  Force  Flight  Dynamics  Laboratory.  The  Air  Force 
Project  Engineer  on  the  contract  was  Mr.  Jerry  E.  Jenkins,  AFFDL/FGC . 

The  report  number  assigned  by  Nielsen  Engineering  & Research,  Inc.  is 
NEAR  TR  26. 

The  authors  wish  to  thank  Mr.  Jenkins  for  his  assistance  in  solving 
seme  of  th«  technical  and  administrative  problems  which  arose  during  the 
course  of  the  investigation.  Also,  they  would  like  to  thank  Mr.  John  C. 
Marshall  and  Mr.  Willard  E.  Summers  of  the  4T  Projects  Branch,  Propulsion 
wind-Tunnel  Facility,  Arnold  Engineering  Development  Center,  for  the  timely 
performance  of  the  experimental  test  program. 

The  work  documented  in  this  report  was  started  on  December  1,  1968  and 
w;*3  effectively  concluded  with  the  submission  of  this  report.  The  report 
was  released  by  the  authors  in  January  1971. 

This  technical  report  has  been  reviewed  and  is  approved. 


E.  H.  FLINN 

Control  Criteria  Branch 
Flight  Control  Division 


ABSTRACT 


This  report  is  the  final  technical  report  which  describes  a combined 
theoretical-experimental  program  which  has  been  conducted  with  the  aim  of 
developing  a computer  program  to  predict  three-degree-of-freedom  trajec- 
tories of  stores  when  dropped  from  fighter-bomber  type  aircraft  at  speeds 
up  to  the  critical  speed.  Both  single  store  and  multiple  store  installa- 
tions are  treated.  The  report  first  describes  the  mathematical  models 
used  to  represent  the  various  aircraft  components.  Then  the  calculation 
of  the  flow  field,  accounting  for  primary  interference  effects,  using 
these  models  is  described  followed  by  the  method  of  calculating  the 
normal  force  and  pitching  moment  from  this  flow  field.  A method  of 
accounting  for  additional  interference  between  the  wing,  pylon,  and 
store  is  next  presented.  Comparisons  between  the  calculated  results  and 
available  experimental  data  and  aata  obta.ned  during  the  present  program 
are  included  for  all  aspects  of  the  analysis.  Good  overall  agreement  is 
obtained.  A comparison  between  one  sample  calculated  trajectory  s"~'  an 
experimental  trajectory  is  presented.  The  agreement  is  quitr  gooc.. 
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everywhere  except  for  section  4 and  for  Cm  of  figures  36 
and  37,  where  it  is  taken  equal  to  store  length,  £s 

store  length 


number  of  control  points  on  pylon 
store  mass 


free-stream  Mach  number 

sum  of  pitching  moments  acting  on  store  taken  about  store 
center  of  gravity,  positive  nose  up 


number  of  control  points  on  left  panel,  or  direction  normal 
to  body  surface 

normal  force 

normal  force  acting  on  store  due  to  primary  interference 

normal  force  acting  on  store  due  to  additional  interference 
vortex  system 

static  pressure  or  number  of  control  points  on  right  panel 
static  pressure  on  lower  surface  of  store 


static  pressure  on  upper  surface  of  store 


dynamic  pressure 

strength  of  k point  source 

Qk/4-rrf^V^  for  bodies,  Qk/2TrcVoo  for  airfoils 


radial  distance  in  y-z  plane 


u,v,w 


u,v,w 


^ 4>  4> 

1 " IT" 


radial  distance  to  i external  vortex 

wing  semispan;  in  section  7,  distance  from  pylon  to  right 
wing  tip 

empennage  semispan 

store  cross  sectional  area,  va 2 

area  of  store  normal  to  sidewash  velocity  V 

reference  area  taken  equal  to  body  frontal  area,  ira^ax 

time  or  wing  thickness;  also,  in  section  7,  distance  from 
pylon  to  left  wing  tip 

perturbation  velocities  in  x , y , z directions, 
respectively  “ 

perturbation  velocities  in  x,  y,  z directions,  respectively 
(fig.  23) 

u/v®»  v/Vw,  w/V^ 


u.  ,v.  ,w. 
1*  l’  X 


u 


us’vs’ws 

Up , Vp  , Wp 

U,  V,  w 


U,  V,  W 

u*,v*,w* 


values  of  u,  v,  w associated  with  vortex  r.  including 
its  image  system 

perturbation  velocity  in  direction  normal  to  line  of 
constant  percent  thickness  of  swept  wing 

values  of  u,  v,  w due  to  store  alone 

sums  of  u^,  v^,  w^  over  all  external  vortices 

total  velocities  in  x_,  y , z_  directions,  respectively 

s s s 

total  velocities  in  x,  y,  ~z  directions,  respectively 

V/Vw,  w/vM 


u ,v  ,w 

s’  s’  s 


us»vs«ws 


VB  * wb 


total  velocities  in  x , y , z directions,  respectively, 
as  seen  from  a point  on  tne  store 

W W WS^.S 

values  of  v and  w for  body  alone  of  conical  wing-body- 
pylon  combination 


P(s) 


P(k) 


value  of  v on  pylon  due  to  store  alone 

value  of  v to  be  cancelled  by  vortex  system  at  k 
control  point  on  pylon 


th 


cro 


ssflow  velocity, 


V* 
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v„/v, 


c <»e 


V.  , V 
hor’  vert 


V_ 


horizontal  and  vertical  components  of  store  velocity, 
respectively 

radial  velocity  in  y -z  plane 

S 5 


v” 


V /V 

jy  oo 


store  ejection  velocity 


3D  (k) 


V 


V„ 


value  of  V associated  with  at  k pylon  control 

point 

aircraft  free-stream  velocity 
ejected  store  free-stream  velocity 


xv 


wL(g) 

WR(j) 

Ww(s) 

w(0 

Ws(a) 


th 

velocity  w to  be  cancelled  by  vortex  system  at  g 
control  point  on  left  panel 

velocity  w to  be  cancelled  by  vortex  system  at  j 
control  point  on  right  panel 

value  of  w on  wing  due  to  store  alone 

complex  potential  for  flow  around  store  in  planes  parallel 
to  ys“zs  Plane 

complex  potential  for  store  alone  perturbation  velocities 


Wr(o) 


x,y,z 

x',y\z’ 

x,y,z 
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★ 
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cs«ys>2s 
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s,  o 
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ZC(X) 


zth(x) 
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complex  potential  associated  with  N external  vortices  and 
their  images  including  the  center  vortex 

Cartesian  coordinate  system  in  compressible  space 

Cartesian  coordinate  system  in  incompressible  space;  related 
to  x,y,z  by  equation  (2) 

Cartesia;.  coordinate  system  with  origin  at  store  nose 
(see  fxg.  23) 

x/£  for  bodies,  x/c  for  airfoils 
x-location  of  kth  point  source 
xk/fR  for  bodies,  x^/c  for  airfoils 

Cartesian  coordinate  system  with  origin  at  nose  of  ejected 
store  (see  fig.  22) 

separation  location  on  store  body 
side  force 

side  force  on  store  alone  due  to  viscous  crossflow 
z/c 

wing  camber  distribution  in  compressible  space 

wing  thickness  distribution  in  compressible  space  measured 
from  the  camberline,  positive  upwards 

angle  of  attack 

a lgle  of  attack  of  fuselage 
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combined  angle  of  attack  a«-  the  axial  location  of  the  tail 
centroid;  that  is,  angle  between  local  velocity  vector  and 
store  axis 

local  angle  of  attack,  tan-1 (W/U) 

parent  aircraft  angle  of  attack;  that  is,  angle  between 
wind-tunnel  longitudinal  axis  and  fuselage  axis 

angle  of  attack  of  store 

angle  of  attack  of  wing 

VT  - or  angle  of  sideslip 

body  or  airfoil  slope  at  the  j point  on  the  surface 

flight  path  angle  of  fuselage 

flight  path  angle  of  store 

strength  of  vortex  at  o = 

store  pitch  angle,  positive  nose  up 
store  pitching  rate,  dt?/'dt 
store  pitching  acceleration,  d^/dt2 
s_/a_  or  expression  given  by  equation  (91) 

t/  L 

sweep  angle  or  expression  defined  by  equation  <113) 

coordinates  of  store  center  of  gravity  relative  to  fuselage 
nose,  see  figure  44 

velocities  of  store  center  of  gravity  relative  to  fuselage 
accelerations  of  store  center  of  gravity  relative  to  fusel  •>>ie 
mass  density 

sidewash  angle,  tan”1  (V/U)  , or  complex  variable  in  y,.-/ 
plane,  - » yg  + izs 

location  of  ft  vortex 

potential  function  in  compressible  space 

complete  potential  for  airplane  including  store  ar.d  pylon 
under  consideration 
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s, cg 
SB 
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th 
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perturbation  potential  of  store  alone 

part  of  potential  0^  due  to  presence  of  store  and  pylon 

d20/dx2 

a20/dy2 

a20/az2 

complete  potential  for  aircraft  except  for  store  and  pylon 
under  consideration 

stream  function  or  algebraic  expression  given  by  equation  (90) 
^/^RV<»  f°r  bodies,  if//cVm  for  airfoils 

Subscripts 

in  the  attached  position 
body 

due  to  buoyancy 
due  to  viscous  crossflow 
two-dimensional  (airfoil) 
of  the  empennage 
leading  edge 

ejected  store  coordinate  system 
store  center  of  gravity 
cv  to  slender-body  theory 
t»i.l  centroid 
tnickness 

of  the  tail  alone  (two  tail  panels  joined  together) 

trail inq  edge 

total 

vortex  lattice 
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A METHOD  FOR  PREDICTING  THREE-DEGREE-OF- 
E‘ RE  EDOM  STORE  SEPARATION  TRAJECTORIES 
AT  SPEEDS  UP  TO  THE  CRITICAL  SPEED 


1 . INTRODUCTION 

This  report  is  the  final  technical  report  describing  a combined 
theoretical-experimental  program  which  has  been  carried  out  with  the 
objective  of  developing  a computer  program  which  will  predict  the 
trajectories  of  external  stores  dropped  from  aircraft  of  the  fighter- 
boniber  type.  The  speed  range  is  limited  to  speeds  below  the  critical. 

The  stores  ca  . be  placed  either  under  the  fuselage  or  the  wing  and  may 
be  single  stores  on  pylons  or  groups  of  stores  on  racks.  Only  store 
motion  in  the  vertical,  or  pitch,  plane  is  considered. 

The  three  principal  tasks  in  the  prediction  of  a store  trajectory 
are:  first,  the  determination  of  the  nonuniform  flow  field  in  the 
neighborhood  of  the  ejected  stores;  second,  the  determination  of  the 
forces  and  moments  on  the  store  in  this  flow  field;  and  third,  the 
integration  of  the  equations  of  motion  to  determine  the  store  trajec- 
tory. The  first  two  tasks  are  the  more  difficult  because  they  are 
complicated  by  interference  effects. 

The  general  methodology  adopted  in  the  present  work  is  that  the 
overall  difficulty  of  the  interference  problem  and  the  probable  length 
of  the  computer  program  make  it  desirable  to  use  the  simplest  possible 
methods  consistent  with  accuracy.  Thus,  even  though  more  precise  methods 
may  be  known  for  handling  some  of  the  aerodynamic  problems,  simplified 
methods  are  used  if  they  are  deemed  sufficiently  accurate  for  the  task 
at  hand. 

At  each  stage  of  development  of  the  prediction  method,  comparisons 
will  be  made  with  experimental  data  to  assess  the  accuracy  of  the  method. 
It  has  been  found  that  there  is  a general  lack  of  data  on  external  stores 
suitable  for  a critical  check  of  a theory  for  predicting  forces  and 
moments  on  a store  in  a nonuniform  flow  field.  Consequently,  an  experi- 
mental program  has  been  conducted  as  part  of  the  present  investigation 
with  the  aim  of  obtaining  such  data.  Data  have  been  obtained  for  single 
stores  and  for  triple  ejection  rack  (TER)  and  multiple  ejection  rack  (MER) 
groupings.  A coordinated  set  of  data  has  been  obtained  giving  the  flow- 
field  velocity  components  where  the  store  is  to  be  located,  forces  and 
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moments  on  the  store  in  this  location,  and  in  most  cases,  pressure  distri- 
butions and  thus  load  distributions  on  the  store  in  this  location.  In 
order  to  isolate  effects  due  to  adding  a component  to  the  parent  aircraft, 
the  parent  aircraft  has  been  built  up  component  by  component  and  the  tests 
repeated  for  each  parent  aircraft  configuration. 

The  data  obtained  during  this  program  are  not  presented  in  this 
report  but  are  tabulated  in  reference  1.  The  complete  test  program  is 
described  in  that  reference.  Representative  data  from  this  program  are 
also  presented  in  references  2 and  3. 


2.  GENERAL  METHOD  OF  APPROACH 

The  accurate  calculation  of  store  separation  trajectories  requires 
that  the  forces  and  moments  acting  on  the  store  be  known  at  every  point 
in  the  trajectory.  If  these  forces  and  moments  are  to  be  calculated  by 
theoretical  means,  then  the  velocity  field  in  the  vicinity  of  the  ejected 
store  relative  to  that  store  must  first  be  determined.  This  field 
depends  not  only  on  the  store  free-stream  velocity,  but  also  the  pertur- 
bation velocities  induced  by  all  of  the  aircraft  components  and  any 
stores  still  attached  to  the  aircraft.  The  determination  of  this  velocity 
field  requires  that  models  for  the  aircraft  components  be  developed  and 
that  the  important  interferences  between  components  be  accounted  for.  In 
the  basic  method  presented  in  this  report,  models  are  developed  for  the 
fuselage,  wing,  and  stores  and  the  following  flow-field  components  are 
accounted  for  in  the  primary  interference  method; 

(1)  Fuselage  volume 

(2)  Fuselage  angle  of  attack 

(3)  Wing  thickness 

(4)  Wing  angle  of  attack,  camber,  and  twist 

(5)  Wing- fuselage  interference 

(6)  Store  volume 

(7)  Wing-store  interference 

Pylon  and  rack  effects  are  not  included  in  the  primary  interference. 

Fuselage  volume  is  taken  into  account  by  approximating  the  shape  by 
an  axisymmetric  body  and  representing  this  body  by  a series  of  point 
sources  placed  on  the  body  axis.  The  volumes  of  all  of  the  stores 
present  are  accounted  for  in  the  same  manner.  Fuselage  angle  of  attack 
is  taken  into  account  as  a simple  Beskin  upwash.  Mutual  interference 
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between  wing,  pylon,  and  store  are  accounted  for  in  an  additional  inter- 
ference method. 

Wing  thickness  is  accounted  for  by  representing  the  airfoil  section 
immediately  above  the  store  being  ejected  by  a series  of  two-dimensional 
point  sources.  Simple  sweep  theory  is  then  used  to  modify  the  velocities 
induced  by  this  source  distribution  to  account  for  wing  sweep  and  taper. 

The  remaining  three  items:  wing  angle  of  attack,  twist,  and  camber; 

wing-fuselage  interference;  and  wing-store  interference  are  accounted  for 
in  the  following  way.  The  wing  is  represented  by  a vortex  lattice  and 
the  velocities  induced  by  the  fuselage  and  stores  at  the  wing  are  computed 
and  treated  as  an  induced  wing  camber  which  is  added  to  any  geometric 
twist  and  camber.  The  wing  model  then  accounts  for  this  induced  camber 
and  the  velocities  induced  by  the  wing  on  the  store  to  be  ejected  are 
computed  with  the  store  in  its  initial  position.  After  the  store  is 
ejected  and  moves  relative  to  the  wing,  the  velocities  induced  by  it  at 
the  wing  will  change.  To  account  precisely  for  these  changing  velocities 
would  require  calculating  a new  wing  model  at  each  point  in  the  trajec- 
tory and  would  greatly  increase  computation  time.  Since  the  initial 
motion  of  the  store  is  the  most  important,  it  is  felt  that  any  change  in 
the  wing  model  due  to  movement  of  the  store  is  a second-order  effect, 
and,  hence,  it  is  not  accounted  for  at  the  present  time.  This  assumption 
can  be  relaxed  if  proven  necessary. 

The  wing  also  changes  the  flow  field  in  the  vicinity  of  the  fuselage 
and  stores  in  such  a way  as  to  cause  them  to  operate  in  a nonuniform 
flow  field.  Since  the  basic  fuselage  and  store  models  are  derived  for  an 
assumed  uniform  flow  field,  the  existence  of  the  nonuniform  flow  causes 
additional  changes  in  the  wing  loading  not  accounted  for  by  the  inter- 
ference effects  discussed  thus  far.  At  the  present  time,  the  modification 
of  the  wing  loading  due  to  this  effect  is  not  included  in  the  primary 
interference  but  is  included  in  the  additional  interference. 

The  fact  that  one  store  can  induce  velocities  in  the  vicinity  of 
another  store  such  as  to  cause  that  store  to  be  in  a nonuniform  flow  field 
is  neglected  for  single  stores,  but  included  for  multiple  stores.  For 
single  stores,  each  mounted  on  its  own  pylon,  the  distance  between  stores 
is  such  that  this  effect  is  felt  to  be  a second-order  effect.  For  stores 
grouped  together,  as  on  TER  or  MER  racks,  experiments  show  that  this  <_ifect 
is  important.  For  the  cases  tested,  the  store  volume  effects  are  dominant 
and  have  been  included  in  the  theory. 
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The  general  method  being  described  in  this  report  allows  the  trajec- 
tory of  a store  ejected  from  under  the  fuselage  or  the  wing,  either  as  a 
single  store  or  one  of  a cluster  on  a TER  or  MER  rack,  to  be  calculated. 

In  the  preceding  discussion,  the  method  of  accounting  for  the  pylon  or 
rack  was  not  described.  In  order  to  investigate  the  importance  of  mutual 
interference  between  the  wing,  pylon,  and  store,  an  interference  model 
has  been  developed  for  the  special  case  of  a single  store  mounted  on  a 
pylon  beneath  the  wing.  This  interference  calculation  has  not  been  carried 
out  for  MER  and  TER  configurations,  but  the  basic  approach  is  applicable 
to  any  of  the  store  locations  and  groupings  of  interest.  The  magnitude 
of  the  additional  interference  for  three-degrees-of-freedom  as  calculated 
is  generally  negligible.  This  may  not  be  the  case  for  six-degree-of- 
freedom  motion. 

The  determination  of  the  trajectory  of  the  ejected  store  requires  that 
the  forces  and  moments  acting  on  the  store,  including  damping,  be  calcu- 
lated at  each  point  in  the  trajectory.  This  is  done  by  calculating  the 
three-dimensional  velocity  field  in  which  the  store  is  operating  as  seen 
by  the  store.  The  free-stream  velocity  and  all  of  the  velocities  induced 
by  the  wing,  fuselage,  and  other  stores  are  summed  up  together  with  that 
due  to  the  angular  velocity  at  the  particular  point  on  the  store.  The 
velocity  field  over  the  store  length  is  calculated  in  this  manner.  The 
axial  distribution  of  the  forces  a.id  moments  due  to  buoyancy,  potential 
flow,  and  viscous  crossflow  are  calculated  from  this  velocity  field  and 
then  integrated  over  the  store  length.  At  the  present  time,  only  normal 
force  and  pitching  moment  are  calculated  although  the  information  for 
calculating  the  side  force  and  yawing  moment  can  be  readily  generated 
by  the  computer  program. 

3.  COMPRESS IBILITY  CORRECTION 

Consider  a wing-fuselage  combination  moving  at  Mach  number  and 

angle  of  attack  of  the  fuselage  a_.  and  assume  that  the  flow  field  for 

D 

compressible  subcritical  flow  is  governed  by  a perturbation  potential  ; 
which  satisfies 


(1  - M2)0  + 0 +0  - 0 (1) 

xx  yy  zz 

where  the  x,y,z  coordinate  system  is  fixed  in  the  fuselage  (see  fig.  1). 
The  approach  to  be  used  in  solving  for  the  wing-fuselage  flow  field  is  to 
transform  the  wing-fuselage  combination  to  an  equivalent  incompressible  one 
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determine  the  incompressible  flow  field  for  the  equivalent  combination,  and 
then  transform  the  flow  field  back  to  the  compressible  speed. 

The  transformation  of  the  compressible  space  (x,y,z)  to  the  incom- 
pressible space  (x',y',z')  can  be  done  in  a number  of  ways.  For  the 
present  work  let  us  use  the  following  transformation: 


x' 


y , 


z ' 


z 


(2) 


The  coordinates  (x,y,z)  and  (x'jy'jZ1)  so  related  are  termed  corresponding 
points.  Using  this  transformation,  the  configuration  cross  sections  in 
corresponding  planes  are  not  changed,  but  the  configuration  is  stretched 
in  the  x direction  to  make  it  more  slender.  Consequently,  the  surface 
slopes  in  planes  parallel  to  the  x-axis  are  reduced  by  the  factor  (3, 
where  £ = vr  - M^.  That  is, 


dz ' r dz 
dx'  “ 0 dx 


(3) 


Let  us  now  determine  the  boundary  conditions  on  the  wing  and  body. 

At  any  streamwise  section,  the  wing  is  defined  by  a thickness  distribution, 
z.vJx),  and  a camber  distribution,  z (x) . The  linearized  wing  boundary 
condition  is  then  for  z = 0 


v4-aB  + 


dz 


th 


dx 


(4) 


The  body  boundary  condition  will  be  written  with  the  aid  of  the  following 
sketch: 
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If  the  body  radius  in  any  plane  perpendicular  to  the  body  axis  is  a(x)> 
then  the  linearized  body  boundary  condition  is 

V»  ("«B  Sin  9 + dt)  = H = U cos  e + H sin  e <-•/ 

Let  us  now  write  ihe  potential  and  boundary  conditions  in  the  incom- 
pressible (x',y',z')  plane.  If  0'  is  the  potential  function  which 
satisfies 

0'  , , + 0'  , , + 0‘  . , = 0 (6) 
x x y y z z’  v ; 

then  the  boundary  conditions  in  the  = 0 space  which  correspond  to 
equations  (4)  and  (5)  are 


If  we  now  let 


= Pac 


(9) 


then  the  incompressible  problem  is  completely  defined  by  the  boundary 
conditions  given  by  equations  (7),  (8),  and  (9),  together  with  the  condi- 
tion that  the  perturbation  velocities  vanish  at  infinity. 

We  must  now  show  how  the  potential  0 is  related  to  0'  at  corres- 
ponding points.  From  equations  (2),  (4),  (7),  and  (9),  we  find  on  the 
wing 


(Id) 
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On  the  body  we  have,  from  equations  (2) , (5) , and  (8) , 


( "“b  sin  0 + 


* v»  (-°b  sin  e' + 


Both  equations  (10)  and  (11)  can  be  satisfied  if,  at  corresponding  points, 


Since  the  relationship  given  by  equation  (12)  will  satisfy  both  the  wing 
and  the  body  boundary  conditions,  and  since  0 will  satisfy  equation  (1) , 
we  have  the  desired  solution  to  the  compressible  problem. 

The  velocities  in  the  compressible  space,  u,  v,  and  w,  are  related 
to  those  in  the  incompressible  space,  u',  v' , and  w' , as  follows: 

„ « M _ A.  30'  = 1 30'  3xc  = 1 30'  _ u* 

3x  8 3x  0 3x'  3x  0^  3x ' 02 

v - M - i 50'  _ 1,  30 1 3v ' 1.  30'  vl  l /•,-,» 

v“3y-p3y-33?r^=337r-0  / (13) 

w = M.  = I Mi  = k Ml  Ml  - 1 Ml  = ill 

* <3z  0 3z  8 3z'  3z  ~ 0 3z'  0 


The  prooedtire  used  to  calculate  the  forces  and  moments  acting  on  a 
store  in  the  presen  '.e  of  a wing- fuselage  combination  for  a compressible 
subcritical  Mach  number  will  now  bt?.  described  m detail.  Consider  the 
wing-fuselage  combination  shown  in  figure  2 and  a single  store  located 
relative  to  it  as  shown  in  figure  3(a).  The  store  geometry  is  shown  in 
figure  3(b).  The  first  step  is  to  transform  the  wing-fuselage-store 
combination  to  an  equivalent  one  for  incompressible  flow.  This  is  done 
by  using  a coordinate  system  whose  origin  is  at  the  fuselage  nose  as  the 
x,y,z  system.  All  of  the  configuration  components  are  located  in  this 
system.  The  transformation  given  by  equation  (2)  is  used  to  determine 
tho  equivalent  configuration.  This  results  in  a lengthening  of  the 
fuselage  and  store  by  a factor  1/0.  The  methods  to  be  presented  in 
section  4.1  are  then  used  to  obtain  source  distributions  to  represent 
the  volume  distribution  oi  these  equivalent  bodies. 
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The  transformation  also  results  in  lengthening  the  wing  chord  by  a 
factor  1/p,  which  results  in  a thinner  airfoil  section.  The  methods  to 
be  presented  in  section  4.2  are  used  to  obtain  the  source  distribution  to 
represent  this  thickness  distribution. 

The  determination  of  the  vorticity  distribution  to  represent  the 
v/ing  surface  is  also  done  in  the  equivalent  plane.  In  addition  to  the 
longer  chord,  the  tangents  of  the  wing  leading-  and  trailing-edge  sweep 
angles  are  increased  by  a factor  1/p.  If  the  wing  is  cambered,  the  slopes 
of  the  mean  camber  surface  are  reduced  by  a factor  p as  is  the  angle  of 
attack  of  the  wing.  The  vorticity  distribution  is  determined  by  the  methods 
to  be  described  in  section  4.3.  Interference  of  the  store  on  the  wing  is 
accounted  for  in  the  determination  of  the  vorticity  distribution. 

The  next  step  in  the  calculation  of  the  forces  and  moments  is  to 
calculate  the  velocity  field  in  the  fuselage  coordinate  system  in  the 
vicinity  of  the  store  for  the  equivalent  configuration.  The  perturbation 
velocities  induced  by  the  fuselage,  wing  thickness  distribution,  and  wing 
vorticity  distribution  are  calculated  in  the  region  that  the  incompressible 
store  would  occupy  if  it  were  present.  They  are  added  together  and  trans- 
formed back  to  the  compressible  space  using  equations  (2)  and  (13)  . 

Finally,  the  free-stream  velocity  components  are  added  to  the  compressible 
perturbation  velocities.  The  example  just  considered  is  a case  where  the 
wing  and  store  arc  at  zero  incidence  relative  to  the  fuselage  axis.  If 
they  are  at  incidence,  the  incidence  angles  are  reduced  by  a factor  of  p 
m the  incompressible  space. 

It  is  important  to  point  out  that  the  incompressible  perturbation 
velocity  field  must  be  calculated  and  resolved  into  velocity  components  in 
the  fuselage  coordinate  system  before  transforming  back  to  the  compressible 
space.  Since  the  transformation  of  the  configuration  was  done  in  this  coor- 
dinate system,  the  velocities  must  also  be  transformed  in  this  system.  Once 
in  the  compressible  space,  the  velocities  can  be  resolved  into  any  desired 
coordinate  system. 


4.  MATHEMATICAL  MODELS  FOR  DETERMINING  FLOW  FIELD 

The  methods  used  to  evaluate  the  forces  on  a released  store  require 
detailed  knowledge  of  the  flow  field  in  the  region  of  the  store  in  its 
absence.  Mathematical  models  of  the  aircraft  components  that  influence 
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this  flow  field  are  described  in  this  section.  The  methods  used  to 
account  for  mutual  interference  among  the  components  are  also  described. 

Models  for  the  thickness  or  volume  distributions  of  the  isolated 
aircraft  components  at  0°  angle  of  attack  have  been  developed  for  the 
wing,  fuselage,  and  stores.  The  wing  thickness  model  is  based  on  a dis- 
tribution of  two-dimensional  point  sources  which  represents  the  thickness 
distribution  of  the  airfoil  cut  out  of  the  wing  by  a vertical  streamwise 
plane  tnrough  the  store  axis.  Adjustments  to  this  model  are  applied 
which  generalize  it  to  three-dimensional  wings  with  sweep  and  taper. 

The  fuselage  and  stores  at  0°  angle  of  attack  are  approximated  by  ax i sym- 
metric bodies  which  are  represented  by  a series  of  three-dimensional 
point  sources  distributed  along  the  body  axes.  Body  angle  of  attack 
effects  are  accounted  for  by  doublets  whose  strengths  are  determined  by 
slender-body  theory. 

A model  composed  of  two  superimposed  planar  arrays  of  horseshoe 
vortices  is  used  to  account  for  the  effects  of  wing  angle  of  attack,  wing 
geometric  camber  and  twist,  and  wing  loading  changes  due  to  the  presence 
of  other  aircraft  components.  The  vortex  strengths  of  the  first  vortex 
array  or  lattice  are  proportional  to  the  angle  of  attack  of  the  wing  and 
account  solely  for  the  wing  angle  of  attack  effects.  The  second  vortex 
lattice,  assumed  independent  of  angle  of  attack,  accounts  for  the  geometric 
twist  and  camber  and  the  normal  velocities  induced  at  the  wing  by  the  fuse- 
lage and  stores , The  velocities  induced  by  these  bodies  are  treated  as  an 
"induced  camber"  which  is  added  to  the  geometric  twist  and  camber. 

4.1  Axisymmetric  Bodies  at  Zero  Degrees  Angle  of  Attack 

4.1.1  Fundamental  equations 

The  stream  function  f and  the  axial  and  radial  velocity  components, 

U and  Vr,  respectively,  due  jointly  to  a free-stream  velocity  aligned 

with  the  body  ax’s,  the  x-axis,  and  a series  of  N point  sources  distri- 
buted along  the  body  axis  are  given  in  reference  4.  With  reference  to 
figure  1,  which  shows  the  coordinate  system,  they  can  be  written  in 
dimensionless  form  as 


N 


= j r*p  - Qj*  <1  + j 
k=i  \ [ 


(x*  - x*) 


(x*  - x*)2  + r*2j 


rn 


(14) 
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N 


i di//*  r~»  Q£(x*  - x*) 

u*(x*,r*)  - ~ SJL  = 1 + ) -f — K—r 

r dr  [(x*  - x*)2  + r*2]3/2 


(15) 


N 


V*  (x*,r*) 


= kr 

r*  dX*  k^i  [(x*  " X*)2  + r*2]3/2 


(16) 


where 


x*  = 


c * 


r*  = 


R 


U*  = 


U 
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V£ 


V 

I 

V 


Q* 


JL 


4r  /- v 
R °» 


•^* 


= __iL 
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and  where  l is  a reference  length  which  for  convenience  will  be  taken 
as  the  length  of  the  body.  The  quantities  q£  and  x£  are  the  source 
strength  and  location,  respectively,  of  the  k*-*1  point  source,  and  the 
point  (x*,r*)  designates  the  field  point  in  cylindrical  coordinates. 

4.1.2  Conditions  used  to  obtain  source  strengths 

A mathematical  representation  of  the  flow  field  about  an  axisvmmetric 
body  is  obtained  by  imposing  three  conditions  on  equations  (13)  and  (16) . 
These  equations,  the  velocity  equations,  are  used  to  calculate  the 
source  strengths  rather  than  the  streamline  equation  because  they  were 

found  to  give  a better  shape  representation.  This  is  due  to  the  fact  that 
they  require  both  body  ordinates  and  surface  slopes  for  their  solution. 

For  the  first  condition  the  flow  directions  at  (N  - 2)  points  (xj»rj>» 
j = 1 , . . . , (N  - 2),  are  specified  by  angles  t3j  with  respect  to  the  positive 
x*-axis;  that  is, 


V*(x*  r*) 
tan  pj  = U*(x*,r*) 


Cor  3 = 1 N - 2 ) 


(17) 
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can  then  be  used  to  determine  the  velocities  at  any  point  in  the  flow  field 
around  the  body. 

4.1.3  Calculative  examples 

Several  arbitrary  choices  remain  in  the  application  of  equations  (19) 
through  (21)  to  particular  bodies.  In  order  to  obtain  adequate  corres- 
pondence between  the  shape  of  the  body  that  is  to  be  represented  and  its 
mathematical  representation,  which  is  the  surface  ^*(x*,r*)  = 0,  a 
suitable  selection  is  required  for  the  number  of  sources  and  the  (N  - 2) 
locations,  (Xj,r^)  of  equation  (19),  where  the  streamline  direction 
constraints  are  to  be  imposed.  Excellent  "fits"  of  several  bodies  have 
been  obtained  to  date  and  the  selection  of  these  parameters  for  two  of 
these  bodies  is  given  below. 

A comparison  between  the  surface  shape  of  a particular  body  (the 
fuselage  of  a wind-tunnel  model  described  in  ref.  5)  and  the  -.'/*(x*,r*)  = 0 
surface  is  presented  in  figure  4(a) . In  examining  this  figure  note  that 
the  r*  scale  has  been  expanded  over  the  x*  scale.  A criterion  that 
has  been  successfully  used  for  this  body,  and  for  all  reasonably  smooth, 
pointed  or  slightly  blunted,  slender  bodies,  is  that  the  distance  between 
adjacent  source  points  be  proportional  to  the  local  body  radius.  For  the 
specific  case  shown  in  figure  4(a) , the  distance  between  source  points 
was  1.2  times  the  local  radius.  For  porportionality  factors  of  the  order 
of  1.5,  irregularities  of  undesirable  magnitude  in  the  •:  * (x*,  r*)  = 0 
surface  tended  to  occur,  whereas,  no  noticeable  improvement  resulted 
from  factors  less  than  1.2.  The  most  forward  point  source  for  the 
example  shown  is  located  at  x*  - 0.002.  For  conically-tipped  forebodies, 
distances  of  this  order  are  required  for  satisfactory  fitting.  For 
values  of  the  first  source  location  significantly  larger  than  0.002,  the 
foresection  of  the  t'*(x*,r*)  =■  0 surface  tends  to  be  blunter  than 
desired,  whereas,  smaller  values,  once  again,  show  no  noticeable  improve- 
ment. For  bodies  whose  aft  sections  are  blunt,  such  as  the  example  in 
figure  4(a),  the  procedure  is  to  extend  the  aft  section  by  a boattail. 

This  procedure  is  based  on  the  fact  that  in  the  real  flow  the  body 
streamtube  separates  from  the  body  at  the  base  and  converges  to  a point 
in  the  wake.  The  extreme  aft  point  source  is  located  at  x*  = 1-118, 
and  the  number  of  sources  used  for  this  example  is  37. 
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The  axial  locations  of  the  (N  - 2)  positions  where  the  flow  velo- 
city directions,  or  surface  slopes,  are  imposed  have  been  chosen  at 
positions  midway  between  adjacent  source  positions.  Since  this  would 
specify  one  too  many  slopes,  a point  in  the  middle  of  the  body  has  been 
e liniinated . 

Results  for  a blunt  axisymmetric  body  are  presented  in  figure  4(b). 
The  body  is  a cylinder  with  ellipsoidal  ends  and  a length-to-maximum- 
diameter  r bio  of  10.  The  ( N - 2)  axial  locations  at  which  the  flow 
velocity  directions,  or  surface  slopes,  were  imposed  are,  as  for  the 
previous  case,  midway  between  source  locations.  For  this  case  no 
precise  criterion  of  source  location  was  applied  other  than  symmetry  with 
respect  to  the  body  center  since  the  fore  and  aft  sections  are  identically 
shaped.  Otherwise,  the  sources  were  concentrated  in  regions  of  high  body 
curvature.  The  total  number  of  point  sources  for  this  case  was  18. 


4.2  Wing  Thickness  at  Zero  Degrees  Angle  of  Attack 

The  flow  about  a symmetrical  airfoil  at  0°  angle  of  attack  is  the 
basic  flow  used  to  account  for  wing  thickness.  Its  mathematical  represen- 
tation is  obtained  in  much  the  same  manner  as  that  for  the  axisymmetric 
bodies.  It  is  then  extended  to  wings  with  sweep  and  taper.  It  should 
be  noted  that  the  flow  parameters  in  this  section  that  pertain  to  the 
symmetric  airfoil  at  zero  angle  of  attack  are  distinguished  by  the 
subscript  " 2d." 

4.2.1  Fundamental  equations 

The  stream  function,  and  t^1e  velocities,  ana  K.,^, 

(fig.  j)  due  to  a uniform  flow  of  velocity  V in  the  positive  x 
direction  with  a senes  of  N two-dimensional  point  sources  distributed 
along  the  x-axis  are  given  in  dimensionless  form  (ref.  4) . 


■Id'**'2*’  ' 


- t Qit 


- tan 


-•  r — a — \\ 

(X*  - X*)  j 


U*D(x*,z*) 


E «( 


(x*  - x*) 

(x*  - x£)r  + z*? 
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x*  = * z*  = 2 
x c » z c 


u*  - -EL  w*  = JL 
u v > w v 
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and  where  the  reference  length  c is  the  chord.  In  these  equations  q£ 
and  x£  are  the  source  strength  and  .location,  respectively,  of  the  kth 
point  source. 

4.2.2  Conditions  used  to  obtain  source  strengths 

The  flow  field  about  a symmetrical  airfoil  is  found  by  imposing 
three  conditions  on  equations  (23)  and  (24)  analogous  to  those  used  for 
the  axisymmetric  case.  First,  at  (N  - 2)  points  (x*,z*), 
j = 1 ,...,( N - 2),  which  correspond  to  airfoil  surface  locations,  the 
flow  direction  is  set  equal  to  the  corresponding  local  airfoil  surface 
slope,  tan  = dz*  (x*  ,z*)  /dx* . Thus  we  obtain 


tan  (3 . = 


-'V*3’2!)  £l 

j U5d(x*,z*> 


or  j = 1 , . . . , (N  - 2) 


The  resulting  expressions  are 


tan  g.  = £ Q£ 


z%  -•  tan  Pj  (x*  - x£) 


(X*  - X*)”  + z*p 


for  j = 1, . . . , (N  - 2)  (26) 


The  second  condition,  the  sum  of  the  source  strengths  equals  zero,  gives 


I -° 


14 


.vi  mcoz»sBiv^n 


The  stagnation  condition  U£D  =0  at  the  leading  edge,  x*  = z*  = 0, 
leads  to  the  result 


N 

I 

k=i 


SE.i 


(28) 


§ 
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The  resultant  set  of  N linear  equations  given  by  equations  (26)  through 
(28)  are  solved  for  the  N unknowns  Q£.  The  velocities  at  any  point 
can  then  be  calculated  using  equations  (23)  and  (24) . 

4.2.3  Calculative  example 

The  foregoing  method  has  been  applied  to  a NACA  65A006  airfoil  and 
the  results  are  shown  in  figure  6.  The  criterion  used  in  selecting  point- 
source  locations  is  that  the  distance  between  adjacent  sources  be  propor- 
tional to  the  local  airfoil  thickness.  This  is  analogous  to  the  criterion 
applied  to  the  first  axisymmetric  body  considered  previously.  The  extreme 
fore  and  aft  source  locations  are  x*  = 0.0006  and  x*  = 0.996,  respec- 
tively, and  the  total  number  of  sources  is  43.  The  axial  positions  at 
which  the  (N  - 2)  velocity  direction  conditions  ..re  imposed,  were 
intermediate  between  adjacent  source  locations,  omitting  the  position  near 
the  mid-chord.  The  comparison  between  the  actual  airfoil  shape  and  the 
calculated  shape  shewn  in  figure  6 shows  good  agreement.  The  velocity 
field  around  the  airfoil  can  be  calculated  using  aquations  (23)  and  (24)  . 

4.2.4  Calculation  of  the  velocity  field  about  wings  with 

sweep  and  taper 

The  calculation  of  the  flow  about  a symmetrical  swept  tapered  wing 
from  knowledge  of  the  flow  about  an  airfoil  following  the  method  of 
reference  6 is  briefly  discussed  in  this  section.  The  geometry  associated 
with  the  problem  is  illustrated  in  figure  7.  The  circular  symbols 
represent  planform  locations  of  points  where  the  flow  velocity  is 
desired.  Lines  of  constant  percent  chord  are  shown,  together  with 
normal  lines  that  are  perpendicular  to  the  lines  of  constant  percent 
chord  at  the  circular  symbols.  Airfoil  sections  along  the  normal  lines 
are  called  normal  sections. 

The  calculation  of  the  thickness-induced  velocities  for  an  infinite 
swept  v/ing  by  simple  sweep  theory  requires  a knowledge  of  the  thickness 
distribution  of  the  normal  airfoil  section  and  the  component  of  the  free- 
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stream  velocity  normal  to  the  leading  edge.  For  swept  and  unswept  wings 
of  finite  span,  in  regions  sufficiently  removed  from  the  wing  root  and 
tip,  the  lines  of  constant  pressure  tend  to  be  parallel  to  lines  of 
constant  percent  fhicKness.  Furthermore,  for  wings  of  constant  thickness 
ratio  and  constant  thickness  distribution  over  the  span  these  lines  are 
parallel  to  lines  of  constant  percent  chord. 

Referring  to  figure  7,  the  calculation  by  simpie  sweep  theory  of  the 
perturbation  velocities  and  w (w  is  the  vertical  velocity  below  the 

point)  requires  that  a source  distribution  be  determined  which  represents 
the  normal  airfoil  section  through  the  point.  Then  the  perturbation 
velocities  can  be  determined  for  the  component  of  the  free-stream  velocity 
parallel  to  this  normal  section,  V^cos  A where  A is  the  angle  between 
the  normal  section  and  the  local  chord.  The  velocity  V^sin  A perpendicular 
to  the  normal  section  does  not  enter  the  calculation  of  un»  Once  ufi  is 
evaluated,  then  the  perturbation  velocities  u and  v in  the  x and  y 
direction,  respectively,  are  simply  given  by 

u = uncos  A (29) 

v = unsin  A (30) 

Designating  and  -*TE»  respectively,  as  the  leading  and  trailing  edge 

sweep  angles , A is  given  by 

an  A * tan  A^  + (tan  A^g  - tan  /^£)  for  0 £ ~ ^ 1 (31) 

where  x/c  is  a fraction  of  the  loca]  chord.  For  locations  forward  of  the 
wing  leading  edge  (or  behind  the  trailing  edge),  A is  taken  as  A 

Lb 

(or  A ) . 

TE 

If  the  above  method  is  to  be  used  for  tapered  wings  and  velocities 
are  required  at  a large  number  of  points  under  the  wing,  then  a large 
number  of  source  distributions  would  be  required,  one  for  each  point, 
since  the  normal  airfoil  section  would  be  different  at  each  point.  A method 
for  calculating  the  perturbation  velocities  without  calculating  ? source 
distribution  for  each  airfoil  section  is  given  in  reference  6.  It  is 
assumed  there  that  the  flow  over  any  normal  section  is  related  to  that 
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over  the  streamwise  section  through  the  point  and  therefore,  only  one 
source  distribution  is  required. 

With  regard  to  evaluating  u^ , consider  the  two  constant  chord  wings 
in  figure  8.  Both  the  swept  and  unswept  wings  have  identical  streamwise 
chords,  streamwise  thickness  ratios,  and  thickness  distributions.  The 
thickness  ratio  of  the  normal  section  of  the  swept  wing  is  larger  than 
that  of  the  streamwise  airfoil  (of  either  wing)  by  a factor  1/cos  n.  it 
is  shown  in  reference  6,  that,  since  the  perturbation  velocities  depend 
linearly  on  thickness  (for  small  thickness  ratios),  the  effect  of  the 
larger  thickness  of  the  normal  airfoil  is  offset  by  the  smaller  imposed 
velocity,  V^cos  A,  in  the  normal  direction.  Thus,  the  perturbation 
velocities  of  normal  sections  of  the  swept  wing  are  equal  to  those  at 
corresponding  points  of  the  unswept  wing.  This  result  holds  approximately 
for  symmetrical  tapered  swept  w-ings  as  well.  Thus,  the  streamwise  and 
lateral  perturbation  velocities  for  the  swept  and  tapered  wing,  u and  v 
! of  equations  (29)  and  (30),  may  now  be  written  in  terns  of  the  perturbation 

1 velocities  of  the  streamwise  airfoil  section  as 


u*  = u*Dcos 

■ (U2D  - 

l)cos 

(32) 

v*  = u^Dsin 

•'  ■ 1U2D  - 

l)sin  A 

(33) 

w*  = 

: w5d  - "2d 

(34) 

For  tapered  wings  the  local  sweep  angle,  , as  given  by  equation  (31)  is 
used  in  equations  (32)  and  (33)  . 

In  summary,  flow  equations  for  the  velocity  components  applicable 
to  swept  and  tapered  wings  with  similar  streamwise  sections  in  the  span- 
wise  direction  are  given  by 


(x*  - x*) 


U* (x*  ,z*)  = 1 + cos 


(x*  - xj»r  + z* 


(35) 


"vW'S.'i'K' 
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V* (x* ,z*)  = sin  A z 


(X*  - X*) 


k^i  (x*  " xk>2  + z*2 


(36) 
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W*( 


(x*,z*)  £ Q* 


z* 


^ (X*  - x£)2  -r  z*2 


(37) 


In  these  equations,  the  source  strengths  Q£  and  their  locations  x£ 
correspond  to  the  local  streamwise  airfoil  section  of  the  swept  and 
tapered  wing  of  interest.  All  quantities  are  made  dimensionless  on  the 
basis  of  the  local  streamwise  chord  of  the  swept  and  tapered  wing  of 
interest. 


4.3  Wing  Angle  of  Attack;  Lifting-Surface  Theory 


4.3.1  Description  of  method 

The  model  used  to  account  for  wing  angle  of  attack,  wing  geometric 
twist  and  camber,  and  interference-induced  twist  and  camber  is  the 
method  described  in  references  7 and  8.  This  method  is  a lifting-surface 
theory  with  the  planform  of  the  wing  represented  by  two  superimposed  arrays 
of  horseshoe  vortices  representing  the  basic  and  additional  lift  distri- 
butions. In  contrast  to  a lifting-line  theory,  this  method  employs 
vortices  distributed  both  chordwise  and  spanwise.  The  wing  planforms 
handled  by  the  method  must  be  representable  by  a series  of  straight- 
line  segments.  The  camber  and  twist  distribution  of  the  wing  is  arbi- 
trary. The  details  of  the  method  will  not  be  repeated  here,  as  they  are 
contained  in  the  above-mentioned  references;  however,  the  following 
paragraph  will  briefly  describe  the  method. 


The  vortex  lattice  is  determined  by  dividing  the  wing  into  area 
elements  formed  by  lines  of  constant  percent  chord  and  lines  parallel 
to  the  root  chord.  A horseshoe  vortex  is  placed  with  the  bound  portion 
at  the  quarter  chord  on  each  of  the  area  elements.  The  two  semi-infinite 
trailing  legs  are  in  the  plane  of  the  wing  and  at  the  side  edges  of  the 
area  elements.  The  boundary  condition  of  no  flow  through  the  wing  is 
applied  at  the  mid-span  of  the  three-quarter  chord  of  each  area  element. 
Two  vorticity  distributions  are  determined,  one  which  is  a function  of 
angle  of  attack  and  the  other  which  is  not.  The  angle  of  attack 
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distribution  cancels  the  component  of  the  free-stream  velocity  normal  to 
the  wing  planform.  The  zero  angle-of-attack  distribution  accounts  for 
the  local  angle  of  attack  due  to  twist  and  camber  when  the  wing  root  chord 
is  at  0 angle  . attack.  The  local  angles  of  attack  at  the  wing  control 
points  due  to  geometric  twist  and  camber  are  input  to  the  computer  program 
and  the  induced  camber  due  to  interference  of  the  fuselage  and  stores  is 
calculated  by  the  program  and  added  to  the  geometric  camber.  Since  both 
vorticity  distributions  are  required  to  satisfy  the  boundary  condition  of 
no  flow  through  the  wing  at  the  control  points,  they  are  obtained  by  solving 
two  sets  of  simultaneous  algebraic  equations  which  differ  only  on  the  right- 
hand  side.  The  computer  program  of  reference  8 then  uses  these  vorticity 
distributions  as  input  to  calculate  the  velocities  induced  by  the  wing  at 
any  point  in  its  flow  field. 

4.3.2  Effect  of  the  number  of  vortices  on  the  calculated  flow  field 

Since  the  wing  is  represented  by  a finite  number  of  horseshoe 
vortices,  the  question  arises  as  to  how  close  to  the  wing  the  velocity 
field  can  be  calculated  within  prescribed  error  limits.  Stores  are  often 
sufficiently  close  to  the  wing  that  the  number  of  vortices  required  to 
obtain  accurate  flow  fields  is  quite  large.  The  purpose  of  this  section 
is  to  provide  a notion  of  the  lattice  size  required  for  reasonably 
accurate  results  at  a given  distance  below  the  wing. 

The  calculations  to  be  shown  have  been  made  for  the  swept  wing 
geometry  of  references  5 and  6,  as  shown  in  figure  S.  The  wing  is  taken 
to  be  at  6°  angle  of  attack.  The  calculations  are  made  for  downwash 
and  sidewash  due  to  angle  of  attack  neglecting  fuselage  and  wing  thickness 
effects.  The  flow  field  results  are  presented  in  terms  of  the  dimension- 
less perturbation  velocities  normal  to  the  plane  of  the  wing  (w/Vo)  and 
in  the  lateral  direction  (v/Vaj)  versus  the  chordwise  distance  (x/c)  at 
the  mid-semispan  of  the  left  wing  panel  for  the  two  distances  below  the 
wing  (z/c)  of  -0.10  and  -0.05.  The  coordinate  system  is  shown  in  figure  10. 
The  perturbation  velocities  u,  v,  and  w are  related  to  the  components 
U,  V,  and  W of  the  local  velocity  vector  as  follows: 


U = V»  CCS  aw  + u 
V = v 


W = V sin  a + w 

00  w 
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The  normal  and  lateral  velocities  at  z/c  = -0.10  for  two  combina- 
tions of  chordwise  and  spanwise  vortices  are  shown  in  figures  11(a)  and 
11(b).  Cyclic  chordwise  variations  evident  for  the  4 chordwise  by  9 
spanwise  (4x9)  array  are  imperceptible  for  the  10x9  array.  Curves  for 
a 7x9  array  were  indistinguishable  from  those  for  the  10x9  array  except 
near  the  wing  leading  edge.  It  appears  evident  from  the  rapid  convergence 
that  curves  from  arrays  with  an  increased  number  of  chordwise  vortices 
would  coincide  with  those  from  the  10x9  array. 

Velocities  at  z/c  = -0.05  for  the  same  two  array,  as  above  are 
presented  in  figures  12(a)  and  12(b).  For  this  case,  small  i nplitude 
fluctuations  persist  for  the  10x9  array.  It  is  evident  from  the  pattern 
cf  convergence,  however,  that  the  10x9  curve  clcsely  approximates  the 
ultimate  curve,  and  that  the  cyclic  variations  would  perhaps  be  impercept- 
ible for  a 13x9  array. 


I 


5.  CALCULATION  OF  FLOW  FIELD 
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5.1  Interference  Effects 

Methods  have  been  presented  for  determining  the  flow  fields  of  the 
fuselage  or  store  alone  and  the  wing  alone.  In  the  case  of  a wing-fuselage 
combination  with  stores,  the  flow  field  is  the  sum  of  the  fields  due  to 
the  various  components  plus  one  due  to  interference  among  the  components. 

In  terms  of  the  potential  functions,  we  have 
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t - potential  due  to  wing  thickness 


= potential  due  to  wing  twist,  camber,  and  angle  of  attack 


m 


■1>_  = potential  due  to  fuselage  volume  and  external  store  volumes 

ts*. 


= potential  due  to  fuselage  and  external  store  angle  of  attack 
i and  sidowash 


= potential  due  to  wing- fuselage  and  wing-store  interference 


20 


- whsm.  »■ fjn«s5*«5s  f^%s  -»v  *^^^*i*^s»***?eii*ek 


f %vf?> 


...*^»S!Kr  ^^^j^4-wrtyc!*'  />'** 


t 


< 

» 


) 

} 

t 

{ 

! 

I 

i 

} 

» 

! 

* 

t 

} 

» 

5 


i 


In  order  to  determine  <£u  , the  only  component  not  considered  so  far, 
we  cam  use  two  means  of  approach,  a noniterative  approach  and  an  iterative 
approach.  In  the  noniterative  approach  singularities  are  placed  on  the 
fuselage  and  the  wing,  and  the  strength  of  the  singularities  are  all 
calculated  simultaneously  using  a large  number  of  simultaneous  equations. 
This  approach  requires  large  computer  times  and  may  be  limited  in  accuracy 
by  the  number  of  simultaneous  equations  the  computer  can  solve.  In  the 
iterative  approach,  adopted  here,  the  computer  time  is  kept  reasonably 
short  and  the  accuracy  of  calculation  is  adequate. 

In  the  iterative  approach  to  determine  the  flow  under  the  wing  of 
the  king- fuselage  combination,  consider  first  the  wing  Mone  at  the 
airplane  angle  of  attack  . Now  let  the  wing  alone  have  velocities  induced 
normal  to  it  by  the  fuselage  flow  field.  The  wing  boundary  condition  is 
no  longer  satisfied,  but  a wing-alone  solution  producing  equal  and  oppo- 
site velocities  to  those  induced  at  the  wing  by  the  fuselage  can  be  found 
from  the  vortex-lattice  program.  This  solution  represents  the  first 
approximation  of  the  interference  solution  in  the  vicinity  of  the  wing. 

To  obtain  the  first  approximation  to  the  interference  potential  in 
the  region  of  the  fuselage,  consider  the  fuselage  alone,  and  then  let  its 
normal  velocity  distribution  be  modified  by  the  induced  flow  field  of  the 
wing  alone.  The  additional  fuselage  solution  required  to  cancel  these 
normal  induced  velocities  Is  the  first  approximation  to  the  inter- 
ference field  in  the  neighborhood  of  the  body.  It  is  clear  that  an 
iterative  scheme  can  be  set  up  to  account  for  higher-order  interactions. 

For  our  purposes  only  one  iteration  is  required,  as  comparison  between 
experiment  and  theory  will  show. 

Some  simplification  has  been  found  possible,  consider  the  applica- 
tion of  the  preceding  equation  to  the  flow  field  under  the  wing.  For  a 
high-wing  monoplane  the  velocities  associated  with  bg  and  both 

produce  velocities  normal  to  the  wing,  alii.ough  for  a midwing  configuration , 
bn.  produces  no  such  velocities.  By  considering  the  normal  wing  veio- 
cities  associar.eo  with  bg  and  0S  , we  can  calculate  interference 
potentials  bj_^  and  bj^  using  the  vortex  lattice  program.  These 
potentials  would  have  the  properties  that 
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where  n is  normal  to  the  wing  surface.  The  physical  argument  could  be 
mace  that  if  the  sum  of  0B^  and  0j_  produces  zero  normal  velocity  at 
the  wing  surface  and  dies  away  with  distance  from  the  surface,  it  will 
produce  negligible  velocities  under  the  wing  except  in  the  neighborhood 
of  the  wing-fuselage  juncture  where  the  first  approximation  to  the  inter- 
ference potential  is  not  accurate.  In  this  region  the  mutual  interference 
problem  must  be  solved.  The  same  argument  can  be  made  for  the  sum  of 

$>Bt  and  ^if 


A calculation  was  made  to  check  this  argumer.c.  The  normal  velocity 

at  the  wing  associated  with  0B  is  the  upwash  produced  by  the  doublet 

distribution  representing  the  body  at  angle  of  attack,  also  termed 

Beskin  upwash.  A twist  distribution  equal  and  opposite  to  that  due  to 

fie  Beskin  upwash  was  introduced  into  the  wing,  with  zero  twist  in  the 

region  blanketed  by  the  body,  and  the  flow  field  associated  Wxtn  the 

resulting  0 j.  potential  was  calculated.  The  downwash  and  sidewash  fields 

resulting  from  the  combined  effect  of  0B  and  was  generally  less 

than  a degree.  Accordingly,  the  assumption  has  been  made  that  the  combined 

effect  of  0B^  and  0^  can  be  neglected.  A precise  calculation  is 

possible,  but  is  not  warranted  at  this  time.  While  the  0R  + 0^  effect 

lias  been  neglected,  the  0R  + 0j  terms  have  been  included.  Both  fuse- 

c t x t 

lege  and  external  store  effects  have  been  included  in  these  terms. 


With  regard  to  the  flow  field  under  the  body,  the  sum  of  0Bfc  and 
satisfies  the  body  boundary  condition  exactly.  The  sum 
0 + 0W  + 0^  must  produce  zero  normal  velocity  at  the  body.  As  we 

move  downward  from  the  body,  the  induced  velocities  due  to  this  sum  must 
decrease.  Accordingly,  wc  have  neglected  tnese  induced  velocities 
everywhere  under  the  body. 

The  simplifications  and  approximations  introduced  in  making  the 
ving-bedy  flow-field  calculation  greatlv  reduce  the  complexity  of  the 
analysis  and  the  amount  of  computer  time  required  at  very  little  loss  in 
accuracy.  The  adequacy  of  the  above  procedure  is  to  be  judged  by  the 
agreement  between  experiment  and  theory  for  downwash  aid  sidewash  subse- 
quently to  be  presented. 
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5.2  Method  of  Calculation 

5.2.1  Single  store  under  the  wing 

The  first  step  in  calculating  the  flow  field  for  the  case  of  a 
single  store  under  the  wing  of  a wing-fuselage  conibination  is  to  obtain 
source  distributions  to  represent  the  volume  distributions  of  the  fuse- 
lage and  store  in  the  incompressible  space  as  well  as  the  wing  thickness 
distribution  in  this  plane.  The  next  step  is  to  calculate  the  vorticity 
distribution  to  represent  the  lifting  wing.  Tve  velocities  induced  at 
the  wing  surface  by  the  store,  and  also  the  fuselage  if  the  wing  is  not 
a midwing,  are  treated  as  an  induced  camber  and  are  added  to  any  geometric 
twist  and  camber. 

The  flow  field  in  the  vicinity  of  the  store  is  calculated  by  summing 
up  the  velocities  induced  by  the  fuselage  and  the  wing.  If  additional 
stores  are  present,  their  influence  on  the  wing  is  accounted  for  as  are 
the  velocities  induced  by  their  volume  distributions  in  the  vicinity  of  the 
ejected  store. 

5.2.2  Single  store  under  fuselage  centerline 

In  calculating  the  velocity  field  under  the  fuselage,  the  presence 
of  the  wing  is  ignored  for  the  reasons  discussed  in  section  5.1.  Thus, 
the  first  step  is  to  obtain  a source  distribution  to  represent  the  fuse- 
lage volume  distribution  in  the  incompressible  space.  For  zero  degrees 
angle  of  attack,  the  velocity  field  induced  by  this  source  distribution 
is  the  field  which  the  store  sees  unless  other  stores  are  present  in 
which  case  the  fields  induced  by  their  volume  distributions  are  additive. 
Fuselage  angle-of-attack  effects  are  treated  as  a Beskin  upwash  and  given 
by 


%Tm  sin  [l  - (jr)  ] l38> 

where  w'  is  the  velocity  normal  to  the  fuselage  axis  in  the  incompress- 
ible space  a distance  z'  from  the  axis,  a is  the  local  fuselage  radius, 
and  is  the  fuselage  angle  of  attack  in  the  incompressible  space . 

This  velocity  field  is  added  to  the  volume-induced  velocity  field. 
Equation  (38)  is  obtained  from  the  complex  potential  for  a circular 
cylinder  in  uniform  flow  given  on  page  29  of  reference  9. 
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5.2.3  Multiple  store  under  the  wing  or  fuselage 

The  flow-field  calculation  for  the  case  of  a TER  or  MER  installation 
is  performed  in  a manner  identical  to  that  for  the  single  store  with  the 
stores  remaining  on  the  rack  treated  as  additional  stores  in  the  flow 
field.  It  is  apparent  from  the  close  proximity  of  the  stores  on  a rack 
that  store-store  interference  may  be  important.  At  the  present  time,  only 
body-volume  effects  are  accounted  for  since  they  are  the  dominant  effects 
as  will  be  seen.  The  source  distribution  for  each  store  is  calculated  as 
if  the  store  were  in  a uniform  flow  by  itself. 

5.3  Comparison  of  Calculated  and  Experimental  Flow  Fields 

References  1 and  5 contain  the  results  of  extensive  experimental 
flew- field  surveys.  Those  presented  in  reference  5 were  taken  in  the 
vicinity  of  the  wing-fuselage  combination  shown  in  figure  9.  The  wing- 
fuselcige  combination  used  in  reference  1 is  shown  in  figure  2.  This 
model  is  a 40-percent  scale  model  of  that  shown  in  figure  9.  Reference  1 
presents,  in  addition  to  data  taken  under  the  wing- fuselage  combination, 
data  taken  with  a pylon  present  at  the  1/3  semispan  location,  a TER  rack 
attached  to  the  pylon,  and,  also,  with  stores  mounted  on  the  two  shoulder 
positions  of  the  rack.  By  building  the  configuration  up  component  by 
component,  as  was  done,  interference  effects  can  be  isolated. 

The  majority  of  the  data  comparisons  presented  in  this  section  will 
use  the  data  of  reference  5.  A few  comparisons  with  the  data  of  reference  1 
will  be  made  primarily  to  show  the  importance  of  knowing  the  flow  character- 
istics of  the  empty  wind  tunnel  when  comparing  with  wind-tunnel  data. 

Since  the  primary  interference  flow-field  calculation  method  does  not 
include  a model  for  the  pylon  or  rack,  comparisons  with  the  data  with  these 
components  present  will  not  be  made.  A comparison  will  be  made  with  the 
perturbation  velocities  induced  by  adding  the  two  shoulder  stores  to  the 
TER  rack. 

5.3.1  Comparison  with  data  from  NACA  RM  L56J19  (ref.  5) 

Reference  5 contains  extensive  experimental  flow-field  surveys  in 
the  vicinity  of  a swept-wing  fuselage  configuration  for  a Mach  number  of 
about  0.15.  The  experimental  wing-fuselage  model  is  that  shown  in  figure  9 
with  the  coordinate  system  and  angle  definitions  shown  in  figure  10.  The 
positive  directions  of  the  local  angles  of  attack  and  sidewash  are  as  indi- 
cated in  figure  10.  It  is  noted  that  positive  sidewash  angle  corresponds 
to  flow  coming  in  from  the  left  when  facing  forward  on  the  wing.  This 


sign  convention  agrees  with  that  used  in  the  present  theory.  The  sign  of 
the  experimental  sidewash  angles  taken  from  reference  5 has  been  changed 
to  conform  with  the  present  sign  convention. 

Comparisons  are  made  for  a very  small  angle  of  attack  at  locations 
below  the  wing  and  remote  from  the  fuselage  to  see  how  the  two-dimensional 
point  source  model  accounts  for  wing  thickness  effects.  Comparisons  under 
the  wing  at  angle  of  attack  serve  principally  to  assess  the  vortex- 
lattice  method.  Finally,  comparisons  under  the  fuselage  for  both  0°  and 
6°  angle  of  attack  serve  to  evaluate  the  methods  used  to  account  for  the 
fuselage. 

Far  all  experimental  data  used,  estimates  of  the  experimental 
accuracy  are  given  in  reference  5 as  +1.0°  and  +1.5°  for  local  angles  of 
attack  and  sidewash,  respectively.  Included  in  the  quoted  experimental 
accuracy  are  estimates  of  probe  alignment  errors  and  wind-tunnel  mis- 
alignment angles. 

For  all  computations  in  this  section  involving  the  vortex-lattice 
model,  a 10  chordwise  by  12  spanwise  array  (on  one  wing  panel)  was 
used.  All  the  flow  angles  were  measured  under  the  left  wing  panel. 

The  source  strength  and  location  distributions  that  have  been 
obtained  for  the  65A006  airfoil  (see  section  4.2.3)  have  been  used  in 
conjunction  with  the  method  of  extension  to  swept  and  tapered  wings 
described  in  section  4.2  4 tc  evaluate  the  flow  under  the  wing  at  0° 
angle  of  attack  for  the  swept  wing  configuration  illustrated  in  figure  9. 
The  calculated  local  angle  of  attack  and  sidewash  angle  results  are 
compared  with  the  data  of  reference  5 in  figure  13  for  three  distances 
below  the  wing.  The  angles  are  plotted  against  axial  distance  from  the 
local  leading  edge  in  percent  local  chord.  The  comparisons  are  at  the 
mid-semispan.  This  location  was  chosen  in  order  that  the  influence  of 
the  fuselage  which  was  present  for  the  measurements  would  be  negligible 
and  thus  the  accuracy  of  the  wing  thickness  model  could  be  assessed.  The 
small  influence  of  the  fuselage  is  illustrated  by  including  for  the 
z/c  = -0.17  case  of  each  figure  the  calculated  angle  including  and 
excluding  the  fuselage  effects.  As  can  be  seen,  they  are  very  small. 

It  should  be  noted,  also,  that  the  measurements  were  made  for  a wing 
(and  fuselage)  angle  of  attack  of  -0.2°.  The  effect  of  this  angle  of 
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attack  on  the  calculated  results,  shown  for  the  z/c  = -0.17  case,  is 
not  significant. 

From  figure  13  it  can  be  seen  that  the  computed  and  experimental  flow 
parameters  are  in  quite  good  agreement,  with  the  largest  differences 
tending  to  occur  ahead  of  and  behind  the  wing.  These  differences  could 
be  caused  by  any  combination  of  the  following  things:  inaccuracies  in 

the  modeling  of  the  three-dimensional  wing  thickness  effects,  viscous 
effects  which  have  been  ignored,  and  wind-tunnel  stream  angles. 

The  combined  effects  of  the  vortex  lattice  and  wing  and  fuselage 
thickness  models  were  examined  over  a free-stream  angle-of-attack  range 
from  -4.3°  to  12.3°,  and  comparisons  between  prediction  and  experiment, 
in  terms  of  the  local  angle  of  attack  and  the  angle  of  sidewash  versus 
chordwise  distance,  are  presented  in  figure  14.  These  comparisons  were 
made  at  the  mid-semispan,  17  percent  of  the  local  chord  below  the  wing 
and  show  how  the  accuracy  of  the  method  varies  with  angle  of  attack. 
Similar  comparisons  at  7,  17,  and  27  percent  of  the  local  chord  below 
the  wing  for  a free-stream  angle  of  attack  of  6.1°  are  given  in  figure  15 
and  exhibit  the  accuracy  as  a function  of  distance  from  the  wing. 

Finally,  comparisons  at  the  quarter,  half,  and  three-quarter  semispan 
locations  for  6.1°  angle  of  attack  at  17  percent  of  the  local  chord  below 
the  wing  are  presented  in  figure  16.  From  these  comparisons  the  accuracy 
as  a function  of  semispan  location  can  be  assessed. 

Figures  14  to  16  exhibit  agreement  between  theory  and  experiment 
which  is  felt  to  be  adequate  for  the  present  investigation.  The  largest 
differences  tend  to  occur  for  the  largest  angle  of  attack  of  12.3° 

(figs,  14(b)  and  14(c)),  the  smallest  distance  below  the  wing  of 
z/c  = -0.07  (figs.  15(a)  anc;  35(b)),  and  the  most  inboard  location  of 
y/(b/2)  = -0.25  (figs.  16(a)  and  16(b)).  This  behavior  is  as  might  be 
expected  since  the  assumption  of  linearity  becomes  less  accurate  as  the 
angle  of  attack  increases.  For  instance,  12.3°  is  beyond  the  linear 
portion  of  the  lift  curve  for  the  wing-fuselage  combination,  while  there 
is  an  assumed  linearity  in  the  vortex-lattice  model.  In  addition,  the 
superposition  of  wing  thickness  effects  and  wing  angle-of-attack  effects 
is  increasingly  in  error  with  increasing  angle  of  attack.  Some  error  is 
developed,  also,  by  excluding  fuselage  angle-of-attack  effects. 
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Some  error  arises  through  the  approximate  wing  thickness  model  used. 

In  figure  14  this  error  can  be  assessed  by  using  the  vortex  lattice  method 
together  with  experimental  thickness  effects  to  determine  the  local  angle 
of  attack  and  sidewash  angle.  These  results  are  shown  by  the  dashed  lines 
of  figure  14.  These  curves  have  been  obtained  by  adding  the  difference 
between  the  experimental  data  and  the  calculation  for  aw  = -0.2°, 
where  the  flow  angles  are  due  principally  to  wing  thickness,  to  the 
calculations  made  for  the  other  angles  of  attack.  As  can  be  seen,  the 
resulting  curves  agree  better,  in  general,  with  the  experimental  data 
than  the  purely  theoretical  solid  curve.  This  agreement  can  be  taken  as 
a measure  of  the  accuracy  of  the  vortex-lattice  method. 

Comparison  of  calculated  and  experimental  local  angle  of  attack  at 
various  distances  under  the  fuselage  axis  for  near  zero  degrees  (-0.2°) 
and  6.1°  free-stream  angles  of  attack  are  presented  in  figure  17.  As  was 
discussed  in  section  5.1,  the  flow  perturbations  due  to  the  wing  thickness 
and  the  vortex  lattice  are  not  included.  The  computations  for  -0.2°  angle 
of  attack  are  based  on  fuselage  thickness  alone;  that  is,  the  small  effect 
of  the  -0.2°  angle  of  attack  has  been  ignored.  The  fuselage  angle-of- 
attack  effects  for  the  case  of  6.1°  angle  of  attack  were  calculated  using 
equation  (38) . This  velocity  is  added  to  the  fuselage  thickness  velocity 
to  obtain  the  local  angle  of  attack. 

From  figure  17(a),  a reasonable  agreement  between  the  calculated  and 
experimental  results  can  be  seen  for  the  case  of  -0.2°  angle  of  attack. 

In  the  region  adjacent  to  the  cylindrical  portion  of  the  fuselage,  a 
fluctuation  in  the  computed  values  of  about  0.2°  amplitude  occurs  for 
the  case  z/a„_  * -1.16.  The  quantity  am=v  is  the  maximum  fuselage 

IH3X  H13X  ^ 

radius.  The  amplitude  is  reduced  for  z/a  = -1.55  to  less  than  0.1  , 
and  is  imperceptible  for  z/amax  ,=  -2.21.  The  origin  of  this  variation 
is,  of  course,  the  fuselage  thickness  model.  The  magnitudes  of  the 
point  sources  representing  the  fuselage  thickness  are  such  that  the 
calculated  body  shape  has  a slight  "ripple"  in  the  cylindrical  section 
(which  is  imperceptible  at  the  scale  of  the  = 0 sketch  given  in 
figure  4(a)),  and  the  adjacent  velocity  field  of  course  has  a corresponding 
ripple  which  diminishes  with  distance  from  the  fuselage.  The  field  points 
at  z/a  = -1.16  are  16  percent  of  the  fuselage  radius  from  its  surface 
in  the  cylindrical  portion,  and  a small  fluctuation  is  therefore  not 
surprising.  While  the  amplitude  of  the  present  variation  is  only  0.2°, 
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this  case  does  indicate  the  desirability  of  an  accurate  fuselage  repre- 
sentation when  stores  in  the  immediate  vicinity  of  the  fuselage  are  under 
consideration.  On  the  basis  of  the  source  distributions  that  have  been 
obtained  for  axisymmetric  bodies  with  cylindrical  sections,  fluctuations 
in  the  cylindrical  region  appear  to  be  characteristic,  and  the  best  pro- 
cedure is  to  minimize  their  magnitude  by  increasing  the  number  and  varying 
the  locations  of  the  sources. 

Considering  still  the  case  of  -0.2°  angle  of  attack,  figure  17(a), 
the  data  lie  slightly  above  the  computed  values  in  the  region  behind  the 
wing,  x/c  > 1.0.  This  is  due  to  small  deviations  between  the  analytic 
representation  of  the  fuselage  shape  and  its  actual  shape  in  this  region. 
Doth  the  data  and  calculated  values  will  commence  to  rise  in  this  region 
since  the  fuselage  radius  beginning  at  x/c  = 1.5  decreases  with  increasing 
x/c.  It  is  indicated  on  figure  4(a)  that  the  p*  = 0 surface  contracts 
more  gradually  than  the  actual  surface,  and  the  calculated  angle  of  attack 
will  be  low  as  a consequence. 

For  the  6.1°  angle-of-attack  case  (fig.  17(b)),  the  agreement  is 
about  as  good  as  that  at  -0.2°  except  at  the  most  remote  distance  of 

z/amax  = ~2-21- 

5.3.2  Comparison  with  data  from  NEAR  TR  24  (ref.  1) 

Figure  18  presents  a comparison  between  the  calculated  flow  field  and 
the  measured  flow  field  using  data  from  Volume  IV  of  reference  1.  The 
wing- fuselage  model  used  in  the  tests  is  shown  in  figure  2 and  the  data 
presented  in  figure  18  were  taken  under  the  left  wing  panel  at  the  1/3 
semispan  location  30  percent  of  the  local  wing  chord  or  2.12  inches  below 
the  wing.  No  pylon  was  present  on  the  wing.  The  Mach  number  is  0.25. 

The  region  of  the  survey  can  be  seen  more  clearly  by  examining 
figure  3(a).  The  purpose  of  the  survey  was  to  measure  the  velocity  field 
in  the  region  the  axis  of  the  store  shown  in  figure  3(a)  would  occupy  if 
it  were  inserted  one  store  diameter,  0.75  inch,  below  the  location  shown 
in  the  figure.  The  data  in  figure  18  are  pJotted  against  (x/£)  where 
i_  is  the  length  of  the  store  shown  in  figure  3 and  x_  is  measured  from 
the  position  the  store  nose  would  occupy  if  it  were  present.  The  directions 
of  the  positive  velocities  are  shown  in  figure  10. 


is  measured  from 


The  directions 
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Data  are  presented  in  figure  18  for  two  angles  of  attack,  ap,  of 
the  wing-fuselage  combination.  This  angle,  ap,  is  the  angle  between  the 
fuselage  axis  and  the  wind-tunnel  longitudinal  axis.  Figure  18(a)  is  for 


a. 


and  figure  18(b)  is  for  ap  = 6 


Two  theoretical  curves  are 


shown  on  the  figure.  One  assumes  an  angle  of  attack  equal  to 


ar 


and 


the  other  assumes  an  angle  of  attack  different  from  ap.  This  second 
case  takes  into  account  the  wind-tunnel  stream  angles  indicated  by  the 
tunnel-empty  survey  made  and  reported  in  Volume  IV  of  reference  1.  These 
data  were  analyzed  in  Volume  I of  reference  1 and  it  was  found  that  for 
aT 


0°  a tunnel-empty  upwash  angle  of  1.25°  existed  in  the  portion  of 


the  tunnel  where  the  survey  of  figure  18(a)  was  taken.  Likewise,  a 
tunnel-empty  upwash  angle  of  0.75°  was  indicated  for  ap  = 6°  (fig.  18(b)) 


The  difference  in  the  tunnel-empty  upwash  angle  for  the  two  values  of  cij. 
is  due  to  the  fact  that  as  ap  is  changed  testing  is  done  in  a different 
part  of  the  tunnel  test  section.  Since  tunnel-empty  surveys  were  only 
taken  at  the  one  spanwise  location,  a uniform  correction  was  made  to  the 
wing  angle  of  attack  used  in  the  vortex-lattice  calculation  rather  than 
using  a wing  twist  distribution  due  to  stream  angle  of  the  wind  tunnel. 


Figure  18(a)  for  ap  = 0V 


indicates  that  better  over-all  agreement 

with  the  experimental  data  is  obtained  when  the  tunnel-empty  stream  angle 

is  accounted  for,  particularly  in  the  region  below  the  wing  leading  edge, 

(x/f)  = 0.05.  This  same  conclusion  cannot  necessarily  be  drawn  from 

s 

figure  18(b)  for  ap  = 6°.  In  some  regions  the  agreement  is  improved  and 
in  others  it  is  made  worse.  On  the  basis  of  these  comparisons,  it  is 
felt  that  to  make  precise  comparisons  with  experimental  flow-field  data 
a knowledge  of  the  angle  of  attack  which  the  wing  actually  sees  in  the 
wind  tunnel  is  required.  The  spanwise  variation  should  also  be  taken 
into  account  by  imposing  a twist  distribution  on  the  wing  in  the  vortex 
lattice  calculation.  In  the  tests  reported  in  Volume  IV  of  reference  1 , 
tunnel-empty  tests  were  not  performed  anywhere  but  at  the  1/3  semi span 
location  nor  were  they  performed  in  the  region  the  wing  chord  plane  would 
occupy  if  it  were  in  the  tunnel.  If  these  data  had  been  taken  and  used 
in  the  theoretical  calculations,  the  overall  agreement  may  have  been 
improved.  The  agreement  exhibited  in  figure  18  is  still  quite  good  for 
the  upwash,  W/V^,  and  sidewash,  V/Vw,  velocities;  it  is  of  the  order  of 
+0.01.  In  terms  of  local  flow  angle,  this  is  +0.58°. 
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Figure  19  presents  a comparison  between  experiment  and  theory  for  the 
same  conditions  as  those  of  figure  18  except  that  the  Mach  number  is  0.70. 
Figure  19(a)  is  for  ap  = 0°  and  figure  19(b)  is  for  ap  = 6°.  In  the 
theoretical  calculations,  the  tunnel-empty  stream  angles  as  given  in 
Volume  I of  reference  1 for  these  test  conditions  have  been  included  in 
the  wing  angle  of  attack.  The  overall  agreement  between  experiment  and 
theory  is  not  quite  as  good  as  exhibited  in  figure  18  for  a Mach  number  of 
0.25.  In  this  case,  the  velocity  ratios  are  predicted  within  +0.02  or  in 
terms  of  local  flow  angle,  +1.2°.  This  difference  may  be  caused  by  not 
knowing  what  free-stream  flow  the  wing  actually  saw  or  by  inaccuracies  in 
the  method  used  to  account  for  compressibility  effects. 


Since  the  primary  interference  flow-field  calculation  method  does  not 

incorporate  models  for  a pylon  and  a rack,  comparisons  with  flow-field  data 

taken  with  these  components  present  will  not  be  made.  Their  presence  will 

be  accounted  for  by  slender-body  theory.  The  method  does,  however,  include 

models  for  other  stores  present  in  the  flow  field  as  in  the  case  of  a 

TER  or  MER  rack.  Consider  the  TER  grouping  shown  in  figure  20  which  was 

tested  in  reference  1.  The  pylon  and  TER  rack  are  not  shown  in  the  figure, 

although  they  were  present  during  the  tests.  The  wing-fuselage  combination 

is  shown  in  figure  2 and  the  store  geometry  m figure  3(b).  By  subtracting 

the  data  obtained  without  stores  2 and  3 present  on  the  TER  rack  from  the 

data  obtained  with  these  two  stores  present,  the  velocity  increments  due 

tc  the  addition  of  the  two  stores  can  be  determined.  These  increments  have 

been  determined  for  one  axial  survey  at  two  angles  of  attack  and  the  data 

are  plotted  in  figure  21.  The  survey  was  conducted  in  the  region  the  axis 

of  store  1 in  figure  20  would  occupy  if  it  was  one  diameter  below  the 

location  shown  ir.  the  figure,  that  is,  -2.83  inches  below  the  wing  chord 

plane.  The  velocities  are  plotted  in  figure  21  against  (x//)  where  f 

s s 

is  the  length  of  the  store  shown  in  figure  20  and  xg  is  measured  from  the 
position  the  nose  would  occupy  if  it  was  present. 


Also  shown  in  figure  21  are  the  results  of  theoretical  calculations. 
These  curves  were  obtained  by  subtracting  the  velocity  field  calculated 
under  the  wing-fuselage  combination  from  that  calculated  under  the  wing- 
fuselage  with  stores  2 and  3 present  including  wing  camber  induced  by  these 
stores.  By  doing  this,  perturbation  velocities  associated  with  wing  angle 
of  attack  and  fuselage  volume  are  subtracted  out  and  only  perturbation 
velocities  associated  with  the  volumes  of  stores  2 and  3 and  wing  camber 
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induced  by  these  stores  remain.  Since  these  are  not  function'’  of  angle 
of  attack,  the  theoretical  curves  for  both  figures  21(a)  and  21(b)  are 
the  same.  A comparison  of  the  data  plotted  in  the  two  parts  of  the  figure 
also  indicates  little  effect  of  angle  of  attack.  The  agreement  between 
experiment  and  theory  is  quite  good.  The  theoretical  lateral  velocity 
increment,  AV/VW,  is  zero,  as  is  borne  out  by  the  data,  since  stores  2 
and  3 induce  equal  and  opposite  lateral  velocities.  In  addition,  because 
of  the  large  distance  from  the  wing,  the  vorticity  distribution  in  the 
wing  which  cancels  the  small  amount  of  store-induced  camber  contributes 
nothing. 


m 
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6.  PRIMARY  INTERFERENCE  FORCE  AND  MOMENT 

The  calculation  of  the  trajectory  of  a store  being  ejected  from  an 
aircraft  requires  the  determination  of  the  forces  and  moments  acting  on 
the  store  at  any  point  in  its  trajectory.  To  calculate  tnese  forces  and 
moments  accurately,  the  nonuniform  nature  of  the  flow  field  must  be 
accounted  for.  In  the  previous  sections  of  this  report,  mathematical 
models  have  been  presented  for  the  various  components  of  the  aircraft 
which  allow  th<»  calculation  of  the  perturbation  velocities  induced  at 
any  point  in  the  flow  field  by  these  components.  The  effect  of  these 
velocities  in  changing  the  forces  and  moments  acting  on  the  store  are  the 
known  dominant  effects  and  will  be  termed  primary  interference.  The  term 
additional  interference  will  be  used  to  identify  the  additional  change  due 
to  mutual  interference  between  the  wing,  pylon,  and  store.  This  section 
of  the  report  will  describe  the  calculation  of  the  primary  interference 
force  and  moment.  Additional  interference  will  be  discussed  in  section  7. 
Generally,  it  is  the  aim  to  put  the  dominant  interference  effects  into 
primary  interference  so  that  the  additional  interference  can  be  neglected. 

By  summing  up  the  perturbation  velocities  due  to  the  aircraft  compo- 
nents and  adding  them  to  the  free-stream  velocity,  the  total  velocities  at 
any  point  in  the  field  where  the  store  will  be  located  can  be  calculated. 

By  performing  this  calculation  for  various  points  over  the  length  of  the 
store,  the  complete  nonunifotra  flow  field  can  be  calculated.  Pitch,  yaw, 
and  roll  damping  can  be  easily  accounted  for  by  adding  the  angular  velocities 
of  the  store  to  the  calculated  flow  field.  From  the  resulting  nonuniform 
velocity  field,  the  load  distribution  on  the  store  can  be  calculated,  and 
the  forces  and  moments  can  then  be  obtained  by  integrating  over  the  store 
length. 
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This  section  of  the  report  describes  the  methods  used  to  calculate 
zb?  forces  and  moments  from  the  known  velocity  field.  Store  forces  and 
moments  due  to  buoyancy,  3 lender-body  theory  {potential  flow)  , and  viscous 
crossflow  are  determined  as  well  as  the  force  and  moment  due  to  the  store 
empennage  if  present.  The  empennage  is  restricted  to  triform  or  cruciform 
tail  arrangements.  The  forces  and  moments  thus  do  not  depend  on  the  roll 
attitude  of  the  store  provided  that  vortex  separation  is  not  present. 

This  restriction  is  not  made  of  necessity,  but  because  it  includes  nearly 
all  of  the  cases  of  interest.  The  damping  due  to  the  pitching  motion  of 
the  store  is  accounted  for  by  adding  to  the  velocity  field  increments  due 
to  pitching  rats.  The  velocity  field  at  any  point  on  the  store  is,  thus, 
that  seen  from  the  store. 


Figure  22  shows  the  xs>ys>zs  coordinate  system  fixed  in  the  store  with 

the  origin  at  the  store  nose.  The  velocities  Us,  Vg,  and  Wg  are  positive 

m the  x , y , z directions.  The  store  velocity  relative  to  the  atmos- 
s s s 

phere  is  Vw  and  the  store  angle  of  attack  is  ag.  The  velocity  compo- 
nents, U , V„,  w are  calculated  at  points  which  would  be  on  the  store 
surface  or  axis  if  the  store  were  present.  They  are 
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US  = V»sCOS  “S  + UVl  + "th 


V = v . + v. . 
s vl  th 


W = V sin  a + w , + w.,  + 0(x  - x ) 

s «>s  s vl  th  ' s s,cg  j 


(39) 


The  quantities  VOT  cos  ag  and  VM  sin  ag  are  free-stream  components. 

The  terms  with  subscript  vl  are  velocities  induced  at  the  point  by  the 
vorticity  distribution  representing  the  wing.  Those  with  the  subscript 
th  are  the  sum  of  the  velocities  induced  by  the  thickness  distributions 
of  the  wing,  fuselage,  and  any  other  stores  present  on  the  aircraft.  The 
term  9(x  - x ) in  the  expression  for  Wg  accounts  for  the  pitching 

motion  of  the  store.  For  points  off  the  body  axis,  a 6 term  should  appear 
in  the  expression  for  U . This  term,  being  very  small,  has  been  neglected 
in  the  present  work.  Yawing  motion  can  be  accounted  for  by  adding  a 
similar  term  to  V . The  velocities  of  equation  (39)  are  made  dimension- 
less by  V to  obtain 
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Us 

u*  = _ = cos  ^ + u*x  + U£h 
s 


V*  = — — = v*  + v* 
s V Vvl  + vth 

oo 

S 


Ws  0(x  “ x ) 

= ■ ■ as  e i w /m  _l_  «.>4  _i_  . mL  _l 


Ws  = V = sin  as  + wtl  + "th  + 


s 


(40) 


It  xs  noted  that  0 is  positive  for  nose— up  motion.  These  expressions, 
along  with  the  work  of  the  preceding  sections  of  this  report,  allow  the 
velocity  field  at  any  point  which  would  lie  on  the  store  surface  or  axis, 
if  the  store  were  present,  to  be  calculated.  The  pitching  rate  of  the 
store,  0,  is  determined  as  a function  of  time  during  the  integration  of 
the  equations  of  motion. 


The  calculation  of  the  forces  and  moments  on  an  external  store  in  a 
nonuniform  flow  field  is  basically  a difficult  problem  which  can  be 
significantly  simplified  by  the  assumption  of  slenderness  valid  for  all 
sections  of  -the  store  except  the  fins  (or  blunt  noses) . The  method  used 
herein  is  similar  to  that  of  McKinney  and  Polhamus  (ref.  10) . For  the 
purpose  of  understanding  the  method  used  herein,  consider  the  variation 
of  downwash  and  sidewash  angles  along  the  centerline  of  a store.  Let  the 
store  be  cambered  in  both  the  vertical  and  lateral  directions  so  that  its 
axis  now  conforms  to  a streamline  of  the  flow  with  its  downwash  and  side- 
wash  components.  Any  lift  or  side  force  due  to  the  downwash  or  sidewash 
velocities  will  be  eliminated  in  this  fashion,  and  the  pressure  field 
acting  on  the  body  will  be,  to  the  first  order,  the  sum  of  its  own 
thickness  pressure  distribution  plus  the  static  pressure  field  acting  in 
the  nonuniform  flow.  The  resultant  forces  (excluding  drag)  and  moments 
of  the  doubly-cambered  body  will  be  due  solely  to  the  static  pressure 
field  of  the  nonuniform  flow.  These  components  of  the  total  forces  and 
moments  are  thus  due  to  buoyancy. 


The  assumption  concerning  the  camber  of  the  store  can  now  be  relaxed, 
and  downwash  and  sidewash  velocities  can  be  considered  to  act  over  the 
length  of  the  axisymmetric  store.  It  is  an  easy  matter  using  slender- 
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body  theory  and  a,  arent-mass  methods  to  evaluate  the  normal-  and  side- 
force  distributions  on  the  store  together  with  uheir  resulting  moments. 
For  stores  subject  to  a lar^.a  value  of  the  combined  angle  resulting  from 
downwash  and  sidewash,  separation  of  the  body  boundary  layer  may  occur. 
Beyond  this  point,  slender-body  theory  is  not  continued.  In  the  present 
calculative  method,  we  have  assumed  that  simple  viscous  crossflow  theory 
is  applicable  downstream  of  the  separation  location. 


6.1  Buoyancy  Force  and  Moment 

The  equations  for  the  buoyancy  normal  force  and  pitching  moment  will 
now  be  derived  for  compressible  flow.  If  the  pressure  coefficient  is 
defined  as 

P ~ P« 


then  it  can  be  expressed  in  terms  of  Mach  number  and  velocity  for  subsonic 
compressible  flow  as  (see  for  example  ref.  11,  page  30,  eq.  (3.33)), 


-[i- 


wf )]  ^ 


“ (us2  + Vf  + Ws2)]  + 


Only  the  first  Mach  number  term  is  shown  since  successive  terms  are  much 
smaller.  Since  in  the  present  work  M^  is  less  than  1.0  and  the  pertur- 
bation velocities  in  equation  (4C)  are  small  compared  to  the  free-stream 
velocity,  this  first  Mach  number  term  is  much  less  than  one.  Consequently, 
in  the  present  work  it  has  been  ignored  and  the  following  incompressible 
expression  used  for  the  pressure  coefficient 


CP  - 1 - K 
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Consider  the  following  sketch 


which  is  the  cross  section  of  the  store  at  any  xg.  The  quantity  u 

is  the  pressure  coefficient  at  the  location  of  the  upper  surface  of  the 

store  (0  * tt/2)  and  C . is  the  pressure  coefficient  at  the  location  of 

s P>  * 

the  lower  surface  (0g  «=  -7r/2)  of  the  store.  Let  the  assumption  be  made  that 
the  vertical  pressure  gradient  is  linear  and  that  there  is  no  lateral 
gradient.  Thus,  as  a function  of  0 

s 

“p  - (Cp,«  - Cp,p’Sin  ®S  l44) 

and  the  buoyancy  normal- force  coefficient  acting  on  an  element  of 
surface  area  is 
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From  equations  (44)  and  (45) 
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In  a similar  manner  the  pitching-moment  coefficient  due  to  buoyancy  is 
found  to  be 
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(48) 


These  two  integrals  are  to  be  evaluated  numerically  using  calculated 

values  of  C . and  C as  a function  of  x . 

P,£  P>u  s 

6.2  Angle  of  Attack  Force  and  Moment 

Let  us  now  derive  the  normal  force  and  pitching  moment  due  to  angle 
of  attack.  To  do  this,  use  will  be  made  of  the  method  presented  in 
reference  12.  In  that  work  it  was  found  that  for  axisymmetric  bodies  in 
uniform  flow  the  gross  normal  force  and  pitching  moment  could  be  calcu- 
lated quite  accurately  for  a wide  range  of  shapes  by  assuming  that  poten- 
tial flow  applies  over  the  forward  portion  of  the  body  up  to  a location 
where  viscous  forces  become  important.  This  location  did  not  change  with 
angle  of  attack  and  is  correlated  in  reference  12  as  a function  of  the 
position  of  the  maximum  negative  rate  of  change  of  cross  sectional  area. 
This  correlation  is  for  a body  in  uniform  flow.  No  correlation  of  this 
type  for  axisymmetric  bodies  in  nonuniform  flow  is  known. 

Using  the  method  of  reference  12  but  accounting  for  the  nonur.iform 
flow  field,  the  normal  force  and  pitching  moment  are  made  up  of  two 
terms,  the  first  calculated  by  slender-body  theory,  or  potential  flow, 
and  the  second  calculated  by  viscous  crossflow  theory.  Slender-body 
theory  is  assumed  to  apply  over  the  forward  portion  of  the  body  to  some 
point,  xg  , where  the  viscous  forces  are  important.  From  this  point  to 
the  end  of  the  body,  the  force  and  moment  are  calculated  by  viscous 
crossflow  theory.  Thus,  the  normal  force  due  to  angle  of  attack  C„  is 
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A similar  expression  can  be  written  for  the  pitching  moment.  The  expres- 
sions for  the  force  and  moment  will  now  be  derived. 

6.2.1  Slender-body  force  and  moment 

It  is  net  difficult  to  determine  the  forces  and  moments  on  a slender 
body  in  a nonuniform  flow  field  using  slender-body  theory.  We  have 
determined  the  buoyancy  contributions  for  the  case  wherein  the  body  is 
cambered  to  conform  to  the  streamlines  of  the  nonuniform  flow  field. 

The  flow  field  can  now  be  made  uniform  retaining  the  cambered  body,  and 
the  forces  and  mome  vts  associated  with  apparent-mass  effects  can  be  calcu- 
lated. Let  us  consider  the  shape  of  the  cambered  body. 

The  local  angle  of  attack  in  the  vertical  plane  is  W*,  and  the  local 
sidewa=h  angle  is  V*  within  t.ie  small-angle  assumption.  Let  xg  be  the 
distance  to  a point  on  the  body  axis  and  let  £g  g " Ys  g + *-zs  g be 
location  of  the  body  axis  at  distance  xg.  Then  we  have 
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The  flow  pattern  is  then  as  depicted  in  the  following  sketch. 


The  complex  potential  can  be  written  for  the  flow  following  the 
results  of  reference  9,  page  29 
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The  complex  force  Y + iN  can  be  very  simply  determined  from  che  fore- 
going complex  potential  by  use  of  the  following  result  from  reference  9, 
page  50.  (The  development  in  reference  9 assumes  a unit  free-stream 
velocity,  v^.  Hence,  the  1/V^  factor  in  the  first  term  of  the  following 
equation. ) 
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wherein 


S(xs)  = va5 
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and  a,  is  the  coefficient  of  the  £ * term  (simple  pole)  in  the 
Laurent  expansion  of  the  complex  potential. 

It  is  easy  to  show  that 
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With  these  relationships  equation  (53)  yields 
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so  that 
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in  simple  terms,  the  normal  force  up  to  axial  distance  xg  depends  only 
on  the  body  radius  and  local  angle  of  attack  at  that  station  with  a 
similar  result  for  side  force.  The  values  of  the  foregoing  quantities 
for  the  store  can  therefore  be  determined  easily  by  evaluating  equations 
(59)  and  (60)  at  the  separation  location.  It  is  of  interest  that  no 
coupling  exists  between  normal  force  and  side  force.  The  pitching-moment 
coefficient  and  yawing-moment  coefficient  associated  with  the  slender- 
body  theory  require  a knowledge  of  the  distribution  of  normal  force  and 

side  force  along  the  external  store  and  require  an  integration  along  the 
body  for  their  determination. 
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Similarly  for  yawing  moment,  we  have 
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The  integrations  in  equations  (61)  and  (62)  will  normally  be  taken  from 
the  store  nose  to  the  separation  location,  xg  Q. 

It  is  noted  that  a positive  normal-force  increment  in  front  of  the 
store  center  of  gravity  causes  a nose-up  pitching  moment,  and  a positive 
side  force  in  front  of  the  center  of  gravity  causes  a nose-right  yawing 
moment.  The  directions  of  these  moments  correspond  to  their  positive 
directions . 

A further  point  of  interest  is  worth  noting.  It  would  have  been 
possible  to  consider  that  the  flow  in  planes  perpendicular  to  the  axis 
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cf  the  store  obeys  Laplace's  equation,  and  that  the  upwash  and  sidewash 
velocities  vary  along  the  store  length.  If  the  potential  is  constructed 
on  this  basis  and  the  body  pressure  coefficients  determined  from  the 
unsteady  Bernoulli  equation,  integration  of  the  body  pressures  to  obtain 
the  forces  yield  both  the  slender-body  results  just  obtained  plus  the 
buoyancy  results,  in  this  case  the  buoyancy  force  is  given  by 
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Equation  (47)  has  been  used  in  the  buoyancy  calculation  rather  than  equation 
(63)  . 

6,2.2  Viscous  crossflow  force  and  moment 

Let  us  now  derive  the  crossflow  normal-force  and  pitching-moment 
coefficients.  It  should  be  pointed  out  that  for  nonuniform  flow  fields 
this  method  is  untried.  For  the  present  application  it  may  have  to  be 
modified  when  data  become  available  to  check  the  method. 
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Referring  to  the  sketch,  we  have  from  crossflow  theory  (ref.  9,  p.  85)  that 
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free  stream. 


; the  section-drag  coefficient  of  a cylinder  normal  to  the 
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so  that  in  coefficient  form 
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Thus,  integrating  from  the  separation  location  to  the  end  of  the  store 
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and  similarly 
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In  the  present  work  cdc  has  been  taken  equal  to  1.2. 

6.3  Empennage  Force  and  Moment 

The  remaining  forces  and  moments  acting  on  the  store  are  those 
produced  by  the  empennage.  The  components  of  these  quantities  due  to 
angle  of  attack  and  damping  are  not  calculated  separately.  As  was  done 
for  the  body  damping,  the  pitching  rate  of  the  store  is  added  to  the 
velocity  field  and  the  total  forces  and  moments  calculated.  The 
velocities  at  the  xg  location  of  the  store  empennage  centroid  are 
calculated  using  equation  (40) . The  combined  angle  of  attack,  the  angle 
between  the  store's  longitudinal  axis  and  the  local  velocity  vector,  at 
the  empennage  centroid  is 
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The  slender-body  theory  normal  force  is 
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where  Kg  and  Kg  are  interference  lift  ratios*  The  quantity  Kg  is 
the  ratio  of  normal  force  developed  on  the  tail  panels  in  the  presence 
of  the  body  to  that  developed  by  the  tail  alone.  The  quantity  Kg  is 
an  analogous  ratio  to  specify  the  normal  force  developed  by  the  body  in 
the  presence  of  the  tail  panels.  These  ratios  from  reference  9,  pages  119 
and  120,  are  given  by 
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where  X = sE/aE«  The  quantity  Sg  is  the  empennage  semispan  and  a£  is 
the  body  radius  at  the  x_  location  of  the  empennage  centroid.  The 
quantity  (dc^j/5a)T  is  the  lift  curve  slope  of  the  tail  alone  and  is  to 
be  specified.  The  tail  alone  is  two  tail  panels  joined  together  without 
the  intervening  body.  The  angle  is  shown  in  the  following  sketch. 
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Equation  (68)  can  be  written 
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The  pitching  moment  due  to  the  empennage  is 
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In  the  above  expressions 
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7.  ADDITIONAL  INTERFERENCE  FORCE  AND  MOMENT 

In  section  6,  primary  interference  forces  and  moments  on  a store 
under  an  aircraft  have  beer,  calculated  from  slender-body  theory  using 
the  calculated  downwash  and  sidewash  of  the  wing-body  combination  at  the 
store  axis.  The  interference  effect  of  the  store  source  distribution 
on  the  wing  was  accounted  for  in  the  attached  condition  and  was  not 
changed  thereafter  during  the  trajectory.  The  effect  of  the  wing  on  the 
store  was  thus  accounted  for  as  in  the  first  step  of  an  iterative  pro- 
cedure. There  is  additional  interference  due  to  the  presence  of  the 
pylon  nd  due  to  the  velocities  induced  at  the  wing  by  the  store  flow 
fields  associated  with  downwash  and  sidewash.  There  is  presently  no 
simple  method  for  determining  the  additional  interference  loading. 
Accordingly,  a numerical  method  is  developed  herein  for  determining  the 
additional  interference  effects,  which  are  directly  additive  to  the 
effects  already  calculated.  Only  the  case  of  a single  store  under  the 
wing  is  considered,  although  the  techniques  employed  can  be  used  to  treat 
other  cases.  The  additional  interference  is  expected  to  be  small  so  that 
slender-body  theory  should  be  sufficiently  accurate  for  its  determination. 

The  methods  available  for  solving  the  additional  interference  problem 
include  conformal  mapping  techniques  and  numerical  vortex  methods.  It  is 
probable  that  with  a considerable  analytical  effort,  some  exact  solutions 


43 


can  be  obtained  using  conformal  mapping  and  slender-body  theory.  However, 
we  have  chcsen  to  use  numerical  vortex  methods  because  of  their  compatibility 
with  the  computer  program  of  the  previous  sections,  the  simpler  analysis 
required,  and  the  possibility  of  applying  them  more  easily  to  a wide  range 
of  interference  problems. 

The  basic  approach  to  the  present  interference  problem  is  to  consider 
first  the  nonuniform  flow  field  which  exists  near  an  airplane  at  the  store 
location  with  the  store  absent.  If  we  now  imagine  the  store  to  appear, 
its  flow  field  induces  velocities  normal  to  the  pylon  and  the  wing  which 
can  be  cancelled  by  bound  vortices  in  all  components.  The  forces  on  the 
store  due  to  the  bound  vortices  represent  the  additional  interference 
effect.  Since  we  have  computed  the  forces  on  the  store  alone  in  the  non- 
uniform  flow  field  using  slender-body  theory,  it  is  consistent  to  use 
slender-body  theory  to  determine  the  interference  effects.  This  approach 
allows  solution  of  the  interference  problem  without  iteration.  The  flow 
is  solved  in  crossflow  planes  normal  to  the  store  axis,  and  as  many  planes 
can  be  used  as  necessary  to  determine  the  store  load  distribution 
accurately.  As  opposed  to  a three-dimensional  approach  which  tries  to 
account  for  all  such  crossflow  planes  simultaneously,  the  present 
approach  has  the  distinct  advantages  of  less  computing  time  and  less 
computer  storage. 

In  the  theory  to  be  developed  in  this  section,  both  primary  and 
additional  interference  terms  are  included  so  that  the  theory  is  complete. 

The  present  approach  is  applicable  to  a finless  store  in  its  entirety  or 
to  a finned  store  up  to  the  beginning  of  the  empennage.  The  empennage  is 
to  be  handled  as  described  in  section  6.3.  In  case  flow  separation  occurs 
on  the  store,  there  is  the  possibility  in  the  present  method  of  including 
the  resulting  shed  vortices  in  the  calculation  as  free  vortices  which 
induce  additional  interference  forces  on  the  store.  The  only  question  is 
with  respect  to  the  strength  and  position  of  the  store  vortices.  In  the 
present  analysis,  only  normal  force  distributions  are  determined,  although 
side-force  distributions  can  easily  be  obtained  from  the  method. 

Not  all  of  the  mathematical  details  of  the  derivation  will  be 
included.  They  are  contained  in  reference  13. 
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7.1  Axes  and  Conventions 

In  the  analysis  several  different  axis  systems  will  be  used.  First 
consider  a plane  fixed  m the  stationary  fluid  normal  to  the  velocity  of 
the  store  as  it  passes  through  the  plane  (fig.  23) . The  cross  section  of 
the  store  in  the  plane  at  any  time  t is  an  ellipse,  the  center  of  which 
is  the  origin  of  the  x,y,z  coordinates.  The  x coordinate  is  aligned 
with  the  velocity,  Vw,  and  the  axes  y,z  lie  in  a plane  normal  to  x. 

Let  the  store  nose  pierce  the  plane  at  t = 0 at  the  origin  of  the  axes 
xQ,y  ,z"0,  which  are  fixed  to  the  plane  and  coincide  with  the  positions  of 

x,y,z’  axes  are  attached  to  the  nose  of  the  store 

and  z being  parallel  to 


x,y,z  at  t « 0.  The 


with  x aligned  in  the  direction  and  y 

v„  and  z\ . 
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Another  set  of  axes  used  in  the  analysis  is  a set  of  axes 


ls’ys»  s 


which  are  symmetry  axes  of  the  store.  These  axes  are  obtained  by  pitching 

the  system  about  y by  an  angle  a,  and  then  yawing  the  system  about  the 

new  vertical  axis  by  -0  as  described  in  reference  9,  page  4.  The  cross 

section  of  the  store  as  seen  in  the  y -z^  plane  is  a true  circle  for  a 

s s 

body  of  revolution. 

It  is  noted  that  the  store  cross  section  in  the  fixed  plane  has  a 
vertical  downward  speed  +W,  and  a lateral  speed  V to  the  left  given 
by  (for  small  angles) 


w = av^ 
v = -pv^ 

The  sign  conventions  of  a and  0 are  such  that  the  store  shown  in 
figure  23  is  at  positive  angle  of  attack  and  negative  angle  of  sideslip. 

Let  the  velocity  components  in  the  x,y,z"  system  be  U,V,W  with 
corresponding  perturbation  components  u,  v,  and  w.  Let  u,  v,  and  w 
be  the  perturbation  velocity  components  in  the  XS>VS>:?S  system.  If  the 
pressure  coefficient  is  defined  as 
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then  the  pressure  coefficient  on  the  basis  of  slender-body  theory  is 

or  in  the  x ,y  ,2  system 
s s s 

c = -2  l^--&y-+,.pw}.  _ (v2  ± wg) 
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7.2  Formulation  of  the  Interference  Problem 

Consider  a wing-body  configuration  which  may  have  a single  pylon- 
mounted  store  beneath  a wing  panel.  Take  a cross  section  through  the 
store  normal  to  its  axis  of  rotation  shortly  after  drop.  The  cross 
section  will  have  the  following  appearance. 


For  the  short  pylons  on  which  such  stores  are  mounted,  we  will  neglect 
sweep  of  the  pylon  leading  and  trailing  edges  so  that  we  are  dealing 
principally  with  cross  sections  of  the  type  depicted  above.  The  portions 
of  the  wing  panel  extending  to  the  left  and  right  of  the  store  will  gener- 
ally be  of  unequal  span.  The  velocities  V and  W are  related  to  the 
local  angle  of  attack  and  sideslip  of  the  store  m the  crossflow  plane  in 
question.  :hese  angles  of  attack  and  sideslip  are  due  in  part  to  the 
downwash  and  siuewash  induced  at  the  store  by  the  wing-body  combination 
(without  the  pylon  and  store  present)  and  to  the  angles  of  pitch  and  yaw 
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of  the  store  itself.  These  latter  angles  are  due  to  the  store  dynamical 
motion  or  to  its  orientation  on  the  pylon  in  the  attached  position. 


In  section  6 of  the  present  *«tudy,  the  lift  and  moment  on  the  store 
were  determined  by  placing  the  store  in  the  wing-body  flow  field  and 
calculating  its  reaction  by  slender-body  theory.  The  empennage  was 
handled  by  linear  theory.  The  store  thickness  effect  was  represented  by 
a distribution  of  three-dimensional  point  sources  on  its  axis,  and  the 
effects  of  V and  W wore  handled  by  slender-body  theory.  Most  of  the 
effects  of  the  wing  on  the  store  are  properly  accounted  for  in  this  manner. 
However,  neither  pylon  interference  nor  the  interference  on  the  wing  due 
to  the  store  flow  fields  associated  with  V and  W are  accounted  for. 

In  this  section  we  will  account  for  these  additional  components  of  the 
total  interference  field  in  a way  which  avoids  iteration.  Also  in  the 
analysis  we  will  include  store  source  effects  in  a way  which  will  permit 
accounting  for  them  as  a function  of  distance  from  the  wing  without 
redoing  the  entire  wing  loading  calculation  at  each  step. 


The  basic  flow  around  the  store  is  due  to  the  nonuniformity  asso- 
ciated with  the  distributions  of  V and  W along  the  store  and  the 
sources  associated  with  changing  body  cross  section.  To  fix  ideas,  let 
us  write  down  the  complex  potential  for  this  flow  field  in  terms  of  the 
complex  variable  a.  The  complex  potential  for  the  perturbation  velo- 
cities v and  w is 
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The  perturbation  velocities  in  the  ys“zs  plane  due  to  this  flow  are 
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v_  - iw_  = 


dc 


s . Wa2  Va2  , „ 

a -X + aV_ 


Vdx  ; a 


(78) 


These  are  thus  components  of  the  store  velocity  field  associated  with  V, 
W,  and  da/dxg  which  produce  velocities  normal  to  wing  and  pylon. 

Because  of  the  normal  velocities,  a distribution  of  bound  vorticity 
is  produced  in  the  wing  and  pylon  which  just  overcomes  the  velocity  compo- 
nents of  equation  (78) . At  the  same  time  an  image  vortex  system  is  intro- 
duced ins’de  the  store  to  maintain  its  cross  sectional  shape  unchanged. 
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i Consider  only  one  vortex  external  to  the  body  at  point  a-.  It  requires 

f . x 

g an  image  vortex  of  opposite  sense  at  a point  internal  to  the  circle, 

. In  addition,  an  image  vortex  of  the  same  sign  is  placed  at  the 

center  of  the  circle  co  keep  the  circulation  at  infinity  unchanged.  As 

a result  the  complex  potential  for  N external  vortices  is 


Wn(a) 


- * i 4- 


J (a  - ai)  - log 


(-O' 


log  a 


I 


The  ri  vortices  are  placed  over  the  wing  and  pylon  in  such  a way 
as  to  counteract  the  velocities  as  obtained  from  equation  (78)  at  certain 
control  points. 


On  the  wing  and  pylon  equation  (78)  yields  the  following  perturbation 
velocities  due  to  the  store.  On  the  wing  a <*  yg  + ih  so  that 


a2W(Ye  - h2)  2azVhy  ^oo*1  <jx 

(h2  + y|)2  (h2  + y2)2  (h2  + y|) 


These  velocity  distributions  on  the  wing  and  pylon  are  to  be  cancelled 
by  the  vortices  associated  vith  Wp(o). 

7.3  Store  Normal-Force  Distribution 

To  obtain  the  force  on  the  store,  we  determine  the  pressure  coeffi- 
cient through  the  use  of  equations  (75)  and  (76)  and  integrate  the  compo- 
nent of  the  surfc.ee  pressure  force  over  the  surface  of  the  store.  It  is 
convenient  to  consider  various  components  in  the  pressure  distriourion. 
Let  ug,  vs,  and  wg  be  the  perturbation  components  due  to  the  store  and 
Up,  Vp,  and  Wp  be  those  due  to  the  N external  vortices, 


ur  “ E ui  ; vr  = E vi  5 wr  = E wi 
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From  reference  9,  page  48, 


(tHtHtHC) 


(83) 


so  that  the  pressure  coefficient  (eqs.  (75)  and  (76))  can  be  written  in  a 
mixed  system  as 
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The  rate  of  change  of  iy  nal  force  with  xs  can  then  be  obtained  by 
integrating  the  above  equation  around  the  surface.  Thus, 
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or 
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<eq.  (86)  cont.  on  next  page) 
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The  first  two  terms  in  this  equation  are  the  primary  interference  asso- 
ciated with  the  nonuniform  flow  field  and  the  last  five  are  the  additional 
interference  due  to  the  vortex  system. 

The  integration  of  equation  (86)  with  respect  to  6 is  carried  out 
in  Appendix  I.  The  resulting  normal-force  distribution  is 


dN 


p ^a  dx 
r 00  00 


s p V a 

r oo  03 


dx  + dx 
s s 


(87) 


The  first  term  on  the  right-hand  side  is  given  by  equation  (1-13)  as 


, dN  , 

1 s _ d 

— = 2t r 


05) 


(88) 


The  second  term,  from  equations  (1-36)  through  (1-38) , is 


dN., 


.,2  dx 
c V a s 
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= 2t 
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(eq.  (89)  cont.  on  next  page) 
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where 


tf'(a>ri>rj»<*,i>0j') 

= t ((rj  " ri)  [(a2  + ri)  r j sin  *j  - (a2  + rj)  ri  sin  *ij 


+ 2r^r j sin  (0 ^ - 0J  £(r|  - a2)rj  cos  0.,  - (r^  - a2)ri  cos  0^j^  (90) 

and  \ in  the  above  expression  is 

* = Cr!  + rj)(a4  + rirj)  ~ 2{a2  + rj)(a2  + ri)  rirj  cos  ^i  " 

+ 4a2rfr^  cos2  {0.  - 0H)  (91) 

1 j 1 j 

The  method  of  calculating  the-  unknown  vortex  strengths  r\  is  given  in 
Appendix  II. 

7.4  Boundary  Condition  for  Store  Under  Pylon  of  Wing-Body 
Combination 

For  the  case  of  an  airplane  represented  by  a wing-body  combination 
with  a number  of  external  stores  under  wing  or  fuselage,  the  method  of 
establishing  the  flow  field  has  been  presented  in  sections  4 and  5 of  this 
report.  Let  this  flow  field  be  described  by  a three-dimensional  potential 
02D.  Consider  a cross  section  of  the  wing-b  *dy  combination  with  a phantom 
pylon  and  a phantom  detached  store. 
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The  potential  0^^  includes  Jie  effect  of  wing  angle  of  attack,  camber, 
twist,  as  well  as  thickness.  It  also  includes  the  effects  of  wing-body 
interference  and  the  volumetric  effects  of  any  other  stores  under  the 
aircraft. 

There  are  associated  velocity  components  V and  W at  the  center- 
line  position  of  the  phantom  store.  Besides  the  component  due  to  03D, 
these  velocities  include  components  due  to  translation  of  the  store 
center  of  gravity,  and  angular  velocities  around  the  center  of  gravity. 
These  factors  are  currently  accounted  for  as  described  in  section  6. 

We  now  consider  the  phantom  store  to  materialize  and  to  add  flows  to  the 
field  associated  with  da/dx,  V,  and  W.  Also,  we  let  the  phantom  pylon 
materialize*  and  it  induces  a further  flow  due  to  the  fact  that  it  is 
subject  to  a sidewash  field  associated  with  and  maY  have  camber 

and  twist  of  its  own.  The  additional  flows  due  to  the  store  and  pylon 
are  associated  with  a slender-body  interference  potential,  $SB*  The 
complete  potential  for  the  problem  is  then 

0c  = 03D  + 0SB  (9 

This  scheme  of  setting  up  the  potential  has  the  advantage  that  the  basic 
three-dimensional  field  is  determined  once,  and  need  not  be  computed  for 
each  new  position  of  the  store  during  the  trajectory.  It  is  feasible  to 
compute  the  changing  part  of  the  field  due  to  the  changing  position  of 
the  store  by  this  means  because  a small  matrix  is  needed  for  a cross- 
flow  calculation  by  the  present  method  (of  the  order  10x10  rather  than 
that  for  a complete  airplane  configuration  including  pylon  and  store 
(of  the  order  200x200) . 
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In  specifying  the  boundary  conditions  on  the  wing  and  pylon  due  to 
$Sb  * we  must  consider  any  normal  velocities  not  already  cancelled  by 


On  the  pylon,  for  instance,  we  have  the  sidewash  velocity  induced 


any  effects  of  camber,  twist,  or  incidence  of  the  pylon,  and 


the  store  flow  field  associated  with  V.  The  store  flow  fields  asso- 
ciated with  W and  da/dx  have  no  component  normal  to  the  pylon.  Thus 
on  the  pylon,  the  boundary  conditior  is 


b<t> 


SB 

dn 


V3D  + V«  + V ^ 
s 


(93) 


3D 

i„ 


sidewash  induced  at  pylon  surface  by  0 


3D 


local  incidence  angle  of  pylon  including 
effect  of  camber,  twist,  and  incidence 


V 


3D 


at  store  axis 


At  the  pylon  control  points  the  boundary  condition  is  thus 

2 


vP(k)  “ V3D(k)  + 1P(k)V»  + 


Va 


hp  - <2k  - !>] 


; k = 1,2, ... ,m  (94) 


With  regard  to  the  right  wing  panel,  0^  includes  any  normal  velocity 
on  the  wing  due  to  its  own  incidence,  camber,  twist,  and  thickness,  as 


well  as  any  fuselage  effects.  If  the  wing  angle  of  attack  is 


a , then 
w* 


the  normal  velocity  V a 


is  accounted  for  by  03Q.  The  velocity  W 


de.termines  the  doublet  strength  on  the  store  axis,  and  the  normal  velocity 
induced  at  the  wing  surface  due  to  this  doublet  is  the  perturbation  velo- 
city which  must  be  cancelled  by  the  interference  vortex  system.  The  store 


flow  field  components  associated  with  V and  da/dx_  also  produce  normal 

5 


velocities  at  the  wing  which  must  be  cancelled.  Accordingly,  we  have  on 
the  right  or  left  wing  panel 

. da 

d*SB 


00  dx. 


?-v^'~,/> t_r^r-’ 1 > ;v_v» ,w. ■'*^*-- •- v’ ?a > - 


At  the  right  wing  panel  control  points,  we  thus  have 
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and  on  the  left  wing  panel  control  points 
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(97) 


These  boundary  conditions,  equations  (94),  (96),  and  (97),  are  used 
in  equations  (11-23)  through  (11-25)  of  Appendix  II  to  determine  the 
strengths  of  the  vortices. 


7.5  Application  of  the  Method  to  Conical  Configurations 
7,5.1  Boundary  conditions 

Since  conical  configurations  can  be  utilized  to  illustrate  the 
importance  of  certain  wing-pylon-store  interference  effects,  some  calcu- 
lated results  for  such  configurations  will  subsequently  be  shown.  The 
boundary  conditions  will  now  be  specified  for  such  configurations. 
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Consider  a conical  store  at  angle  of  attack  a and  angle  of 
sideslip  (3  with  crossflow  planes  normal  to  the  store  axis  as  shown 


Consider  now  the  positions  which  are  to  be  occupied  by  wing  and  pylon. 

At  these  positions  the  store  alone  induces  velocities  by  virtue  of  da/dxs, 
V,  and  W.  In  order  to  determine  the  velocities  normal  to  the  pylon  and 
wing,  consider  the  total  potential  for  the  store  flow 

*s  * R*p*[v  (a  + 4) + iw (“°  + t0+  aV»  asr log  °] 


+ ' v2  - 


vr 


(98) 
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v(1-7)"iw(lt7)+  av»^  ° 

On  the  pylon  where  o *=  izs,  we  bave 

(-5) 


Vp  = V 


(99) 


(100) 


The  pylon  control  points  have  the  coordinates  (eq.  (II-3) ) 


>v  - i ‘2k  - » 


so  that  equation  (100)  becomes 
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This  is  the  velocity  to  be  cancelled  by  the  vortex  system  at  the  m 
control  points  on  the  pylon. 

On  the  wing,  we  have  the  value  of  wB  as  given  by  equation  (99) , but 
we  also  have  a normal  velocity  associated  with  the  potential  term 
x V^E  - V2  - W2  unless  the  wing  is  parallel  to  xs.  Consider  a side 
view  of  the  conical  configuration. 


Pylon 


Wing 


"Store 


Neglecting  V2  + W2  compared  to  V®,  the  normal  component  due  to  the 

root  term  is  -V  (dh/dx  ) . On  the  wing  a = y„  + ih  so  that 
w s s 
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from  which 
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At  the  right  wing  panel  control  points  (see  eq.  (II-3)) 


ys  “ 2p  (2j  - i)  ; j = i,2, ...,p 
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Thus  the  velocity  to  be  cancelled  by  the  vortex  system  at  the 
points  on  the  right  wing  panel  is 
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At  the  left  wing  panel  control  points. 
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so  that  the  velocity  to  be  cancelled  at  the  n control 
panel  by  the  vortex  system  is 
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is  of  interest  to  consider  the  boundary 
basic  cases: 

Case  A:  v = 0 w = o ~ ¥ 0 

rtv  ' 


conditions  as 


w = 0 ~ 7*  O ~ ¥ 0 

dx  T dx  T 


Case  B:  V / 0 W * 0 


Case  C:  V = 0 W / 0 


The  first  case  yields  quantities  due  to  the  basic  asymmetries,  and 
establishes  the  angles  of  zero  normal  force  and  zero  side  force. 
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For  Case  A,  the  boundary  conditions  are  symmetrical  about  the 
vertical  axis,  and  a vortex  system  with  mirror  symmetry  about  the 
vertical  axis  is  to  be  expected  for  left  and  right  panels  of  equal 
span  with  no  vortices  on  the  pylon.  For  Case  B,  the  vortex  signs  will 
not  change  between  left  and  right  panels  for  equal  span  panels,  and 
there  will  be  vortices  on  the  pylon.  The  strengths  will  be  equal  on 
left  and  right  panels  and  proportional  to  the  sidewash  angle.  For 
Case  C,  the  vortices  will  be  similar  to  those  for  Case  A.  The  physical 
significance  of  these  cases  will  subsequently  be  discussed  further. 

7.5.2  Calculative  example 

In  order  to  demonstrate  the  calculative  method  and  to  illustrate  the 
nature  of  wing-store-pylon  interference,  a number  of  systematic  calcula- 
tions have  been  carried  out  for  conical  configurations  of  the  general 
cross  section  shown  below. 
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It  is  possible  to  show  systematically  the  effects  of  the  nondimensional 
quantities  h/a,  s/a,  and  f/a  on  the  interference  for  such  configurations, 
and  thereby  to  shed  light  on  the  magnitude  of  the  interference  effects. 

For  a conical  configuration  in  a uniform  stream,  the  right-hand  sides 
of  equations  (88)  and  (89)  do  not  vary  with  xs.  Designating  the  part 
associated  with  the  vortex  system  as  A,  we  have  from  equation  (87) 
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dx  dx  P°°  00 


(113) 


If  we  let  the  conical  configuration  have  a store  base  radius  _ 

IuaX 

a length  ig  and  use  the  store  base  area  as  reference  area,  then 
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The  quantity  CN  is  obtained  from  equation  (1-13)  and  is  found  to  be 
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The  quantity  ip  is  given  by  equations  (90)  and  (91)  . 

The  boundary  conditions  for  conical  configurations  were  developed  in 
the  previous  section,  wherein  three  cases  were  differentiated.  Case  A 
included  all  terms  not  proportional  to  a or  3;  Case  B included  terms 
proportional  to  3 (or  V) ; Case  C included  terms  proportional  to  a 
(or  W) . Case  A determines  the  lift  and  side  force  at  zero  a and  3, 
and  is  not  as  interesting  as  the  other  two  cases. 

Consider  the  angle-of-attack  case  associated  with  vertical  velocity 
W.  The  store  alone  has  a normal-force  curve  slope  of  2 for  this  case 
based  on  its  base  area  as  reference  area.  The  presence  of  the  wing  above 
the  store  tends  generally  to  dam  up  the  crossflow  above  the  store,  and 
thereby  to  reduce  the  normal  force.  If  we  designated  the  normal-force 
coefficient  of  the  store  alone  due  to  a as  CM  . . and  that  of  the 

, we  have 


'Ns  (a) 

store  in  the  presence  of  the  wing  and  pylon  as  C, 
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The  interference  factor  Kq  has  the  simple  physical  significance  that 
it  is  the  fractional  amount  that  the  store-alone  normal  force  due  to  a 
is  changed  because  of  wing-store-pylon  interference. 

Inherent  in  the  foregoing  physical  interpretation  is  the  assumption 
that  the  nonlinear  terms  in  equation  (116)  have  no  net  contribution  to 
CN  for  the  angle-of-attack  case.  Since  the  vortex  strengths  are  propor- 
tional to  W in  this  case,  the  last  three  terms  of  equation  (116)  yield 
contributions  quadratic  in  a.  For  angles  of  attack  of  the  order  of  a 
tenth  of  a radiar,  the  numeric?!  contributions  from  the  nonlinear  terms 
generally  total  less  than  one  percent  of  the  linear  term.  Whether  the 
net  contribution  of  tne  nonlinear  terms  should  be  identically  zero  on 
mathematical  grounds  is  not  known  to  the  authors  at  this  time.  The 
residual  nonlinear  contribution  could  be  the  result  of  approximating  a 
continuous  vorticity  distribution  by  a number  of  discrete  vortices.  In 
any  event,  it  is  small  for  the  magnitude  of  the  angles  expected  in  practice, 
and  it  is  therefore  neglected. 

A series  of  calculations  has  been  made  to  determine  how  Ka  varies 
with  the  nondimensional  parameters  h/a  and  s/a.  Because  the  pylon  is 
a stieamline  of  the  flow,  no  pylon  vortices  are  required  in  the  solution, 
and  Ka  does  not  depend  on  £/a.  Since  we  are  replacing  a continuous 
vorticity  distribution  on  the  wing  by  a series  of  discrete  vortices,  the 
accuracy  of  the  numerical  results  depends  on  the  spacing  of  the  vortices. 
The  variation  of  Kq  with  h/a  for  s/a  =2  is  shown  in  figure  24  for 
two  uniform  vortex  spacings.  For  most  practical  purposes,  a vortex 
spacing  equal  to  a quarter  of  the  body  radius  gives  adequate  results. 
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To  illustrate  further  the  convergence  of  the  method  for  calculating 
K a systematic  set  of  calculations  was  made  varying  the  vortex  spacing 
tor  h/a  = 1.25  and  s/a  » 10.  The  following  resilts  were  obtained: 


Vortex  Spacing 


K 


a 


a 

-0.57101 

2 a/3 

-0.56515 

a/ 2 

-0.56238 

2a/5 

-0.56072 

a/3 

-0.55958 

2 a/7 

-0.55876 

a/4 

-0.55816 

Tor  this  case,  a vortex  spacing  as  great  as  a body  radius  gives  adequate 
results  for  most  practical  purposes. 

The  effect  of  h/a  on  K for  various  values  of  s/a  is  shown  in 

a 

figure  25.  For  the  limiting  case  s/a  = 0,  Kq  is  clearly  zero.  For  the 
iimitirj  case  s/a  = <*-,  it  is  pro.  able  that  K -1  based  on  the  following 
aigument.  For  a very  large  value  of  s/a,  the  flow  under  the  central  part 
of  the  wing  is  a stagnation  region.  Placing  a store  of  radius  small  com- 
pered to  s ir.  this  region  will  product  no  normal  force  on  the  store,  so 

-hat  K will  be  -1. 

a 

i.  significant  point  in  connection  with  the  angle-of-attack  case  is 
that  no  singularities  arise.  Such  singularities  occur  when  an  external 
vortex  comes  in  contact  with  the  store.  It  is  noted  in  equation  (116)  that 
the  third  term  produces  a singularity  if  r^  = a.  From  symmetry  consider a- 
Mons,  it  is  noted  that  the  vorticity  on  one  wing  panel  is  equal  and  oppo- 
site to  the  other  wing  panel.  Accordingly,  no  trailing  vortex  lies  along 
cue  root  chord  of  the  wing.  Thus  when  the  store  is  in  contact  with  the 
wing,  there  is  no  vortex  in  contact  with  the  store  to  produce  a singularity. 
It  should  be  noted  that  as  soon  as  the  left  and  right  wing  panels  are  of 
unequal  span  or  • s roon  as  sidewash  is  introduced,  the  singularity  arises. 


For  the  angle-cf-sidewash  case,  the  last  three  terms  of  equation  (116) 

Me  quadrati:  in  3,  and  the  linear  term  has  no  contribution  to  C^.  It 

is  convenient  to  normalize  vhe  3’’  contribution  to  Ch,  in  terms  of  a 

N 

nondimensional  normal- force  ratio,  , which  has  a physical  analogy  with 
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a crossflcw  drag  coefficient.  Corsider  the  side  force  on  the  store  alone 
due  to  viscous  crossflow  with  crossflow  drag  coefficient  cdc* 

Ycf  " cd  P2^Sn  <120^ 

c 

For  a conical  body  the  area  normal  to  the  crossflow,  ? , is  It 

is  possible  to  put  the  normal  force  due  to  02 , N^,  as  obtained  from 
equation  (116) , into  the  following  form 

Ni3  = V?<i»Sn  <121> 

The  factor  has  a simple  physical  significance.  It  is  numerically 

equal  to  the  crossflow  drag  coefficient  which  will  yield  a side  force 
due  to  viscous  crossflow.  This  side  force  is  equal  to  the  normal  force 
Np  due  to  the  angle  of  sidewash. 

In  figure  26,  the  value  of  is  shown  versus  h/a  for  s/a  = 2 

and  g/a  - 1 for  two  uniform  vortex  spacings  on  wing  and  pylon.  It  is 
noted  that  good  agreement  occurs  for  a gap  between  pylon  and  store  as  small 
as  half  a radius,  however,  as  the  gap  gets  smaller,  the  two  results 
start  to  depart  significantly.  The  difference  is  associated  with  the  fact 
that  as  the  lowest  pylon  vortex  gets  closer  to  the  store,  it  approaches  a 
singular’’  y of  the  type  previously  discussed. 

A systematic  set  of  calculations  has  been  made  to  see  how  varies 

with  vortex  spacing  for  the  same  case  for  which  the  K calculations  were 
made.  The  following  results  were  obtained: 


Vortex  Soacinq 

a 

2.2996 

2a/3 

1.6092 

a/2 

1.5043 

2a/5 

1.4892 

a/3 

1.4869 

2 a/7 

1.4864 

a/4 

1.4864 

For  this  case,  with  h/a  = 1.25,  s/a  = 10,  and  f/a  = 1.0,  a vortex  spacing 
of  a/2  yields  sufficiently  accurate  results  for  most  practical  purposes. 
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An  exact  solution  for  has  been  obtained  for  the  case  of  s/a  = <*  and 

i/a  = 0 with  the  help  of  reference  14,  page  174,  example  25,  It  is  of 
interest  to  compare  the  values  of  for  this  case  with  those  obtained 

by  the  present  calculation  method  for  s/a  = 10  as  a function  of  h/a. 

The  comparison  is  shown  balow. 


h/a 

s/a  = » 

Exact  Solution 

s/a  = 10 
Present  Method 

1.05 

6.2285 

7.5549 

1.10 

3.5992 

3.6009 

1.15 

2.5090 

2.4806 

1.20 

1.8978 

1.8750 

1.25 

1.5042 

1 . 4b640 

1.50 

.6550 

.64724 

1.75 

.3658 

.360/4 

2.00 

.2298 

.22662 

2.50 

.1102 

.10802 

4.00 

.0254 

.0239 

5.30 

.0128 

.0114 

8.00 

.0031 

.00205 

fiven  for  large  values  of  h/a,  where  the  difference  in  s/a  may  have  <^ne 
effect,  the  results  are  closely  the  same.  For  small  values  of  h/~  where 
the  difference  in  s/a  is  of  less  significance,  there  is  still  very  close 
agreement  except  at  h/a  = 1.05.  At  h/a  = 1.05,  the  effect  of  the  singu- 
larity at  h/a  =1.0  in  the  present  method  is  starting  to  make  its  influ- 
ence felt.  The  vortex  spacing  for  the  calculation  was  a/4. 

The  general  nature  of  the  flow  is  as  shown  in  the  sketch  of  figure  26. 
High  velocity  flow  through  the  gap  creates  a lower  pressure  there  than  at 
the  bottom  of  the  store  so  that  a positive  normal  force  is  developed. 

As  the  gap  approaches  zero,  the  lowest  pylon  vortex  does  not  actually 
touch  the  store  because  of  the  way  in  which  the  pylon  vortices  have  been 
distributed.  The  lowest  vortex  has  been  placed  one-half  of  uhe  pylon 
vortex  spacing  up  from  the  pylon  tip.  This  convention  was  followed  for 
both  vortex  spacings.  Ko  analysis  was  made  to  determine  how  the  lowest 
vortex  should  be  located  with  respect  to  the  tip,  and  the  above  conver- 
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tion  represents  a first  approximation.  Refinement  of  this  convention  would 
probably  produce  better  convergence  as  the  tip  is  approached. 

Another  theoretical  point  of  significance  is  that  tha  limit  of  the 
present  case  as  the  gap  approaches  zero  is  not  the  same  as  the  zero-gap 
case.  In  the  present  case,  the  loading  at  the  pylon  tip  falls  to  zero, 
and  as  the  gap  approaches  zero  this  condition  is  not  relaxed.  However, 
if  the  gap  becomes  zero,  then  a lo  .ng  can  be  carried  across  the  tip  of 
the  pylon.  It  appears  that  the  positive  normal  force  for  the  gap  case 
could  change  to  a negative  normal  force  for  the  zero-gap  case.  For  the 
purpose  of  calculating  trajectories,  we  are  interested  in  the  gap  case 
because  the  gap  opens  up  immediately  whenever  the  store  is  dropped  or 
ejected.  However,  for  determining  attached-store  loads,  the  zero-gap 
case  could  be  of  interest  in  special  cases. 

Figure  27  has  been  prepared  -ce  show  how  Kg  varies  with  h/a  for 

constant  values  of  s/a.  It  is  intuitively  clear  that  after  s/a  reaches 

a certain  value  for  a fixed  value  of  f/a,  it  would  not  be  expected  that 

Kg  would  change  significantly  for  further  increases  in  s/a.  Calculations 

bear  out  this  supposition  in  that  going  from  s/a  = 8 to  s/a  = 10  aid 

not  change  the  calculated  values  of  Kg  significantly.  In  the  foregoing 

sense  there  is  an  upper  limit  on  Kg  for  a given  value  of  h/a.  The  figure 

illustrates  the  rule  of  thumb  that  KQ  is  0.1  or  les1-  if  h/a  * 4 to  5 

P 

for  a short  pylon. 

8.  COMPARISON  OF  CALCULATED  LOAD  DISTRIBUTION,  FORCE,  AND 

MOMENT  WITH  EXPERIMENT 

This  section  of  the  report  will  present  comparisons  between  the 
present  theory  and  experimental  data  in  order  to  assess  the  accuracy  of 
the  theory  in  predicting  the  store  load  distribution  and  also  the  total 
normal  force  and  pitching  moment  acting  on  the  store.  Comparisons  will 
first  be  made  for  the  case  of  a single  store  under  the  wing  of  a wing- 
fuselage  combination.  These  will  be  followed  by  comparisons  for  a single 
store  under  the  pylon  of  a wing- fuselage-pylon  combination.  The  last 
comparisons  will  be  for  a TER  grouping  under  the  wing  of  a wing-fuselage- 
pylon-rack  combination. 
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8.1  Single  Store  Under  Wing  of  Wing-Fuselage  Combination 

Some  calculated  normal-force  distributions  for  the  case  of  a single 
store  under  the  left  wing  panel  of  a wing- fuselage  combination  are  shown 
in  the  next  three  figures.  The  calculations  and  the  data  shown  in  the 
figures  are  for  the  configuration  shown  in  figures  2 and  3 except  that 
the  store  is  oue  diameter,  0.75  inch,  below  the  location  shown  in 
figure  3(a).  This  is  the  same  case  as  for  the  flow-field  calculations 
which  were  presented  in  figures  18  and  19.  The  data  taken  from  Volume  III 
of  reference  1 are  shown  by  the  circles. 

Figure  28  compares  the  primary  interference  normal  force  distribution, 
calculated  as  described  in  section  6.  with  the  data  taken  at  a Mach 
number  of  0.25  at  two  angles  of  attack.  The  angle  ap  is  the  angle  at 
which  the  store  and  parent  aircraft  were  placed  relative  to  the  wind- 
tunnel  longitudinal  axis.  The  angle  a was  used  in  the  calculations. 

As  was  discussed  in  section  5.3.2,  this  angle  reflects  a partial  correction 
for  the  wind-tunnel  stream  angle.  Two  theoretical  curves  are  shown  in 
figure  28.  The  solid  curve  in  the  figure  is  the  primary  interference 
normal-force  distribution  excluding  store  interference  on  the  wing;  that 
is,  the  store  volume  is  not  allowed  to  induce  a wing  camber.  For  the 
dashed  curve,  this  induced  camber  is  included.  As  can  be  seen  in 
figure  28,  the  effect  of  the  .nduced  camber  is  small.  At  both  angles  of 
attack,  the  agreement  between  the  experimental  and  calculated  distributions 
is  quite  good. 

The  calculations  shown  by  the  dashed  curves  in  figure  28  account  for 
store  interference  on  the  wing,  but  not  wing  interference  on  the  store, 
whereas  the  additional  interference  analysis  of  section  7 accounts  for 
the  mutual  interference  between  these  two  components.  Figure  29  shows 
the  effect  of  including  the  additional  interference,  as  calculated  by  the 
method  of  section  7,  in  the  load  distribution.  These  are  the  same  cases 
as  treated  in  figure  28.  The  solid  curves  are  the  primary  interference 
load  distributions.  The  store-induced  wing  camber  has  not  been  included 
since  the  store  source  distribution  has  been  included  in  the  boundary 
condition  of  the  additional  interference  calculation,  the  da  .lx  terms 
m equations  (96)  and  (67) . The  dashed  curves  in  figure  29  include  the 
additional  interference.  The  only  effect  seen  is  m the  nose  region  of 
the  store.  The  effect  appears  to  be  much  more  pronounced  for  !p  - 0° 
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(fig.  29(a))  than  for  ap  = 6°.  This  is  only  because  of  the  fact  that  the 
large  slopes  of  the  load  distribution  curve  which  exist  near  the  nose  for 
the  Op  = 6°  case  hide  the  change  due  to  additional  interference. 

For  exp  = O'1,  the  additional  interference  calculation  produces  an 
effect  which  is  not  seen  in  the  data.  This  effect  is  due  entirely  to  the 
store  source  terms,  the  da/dxs  terms,  in  equations  (96)  and  (97) . If 
these  terms  are  not  included  in  the  boundary  condition,  the  change  due  to 
additional  interference  cannot  be  seen  on  the  plots  of  figure  29.  It, 
therefore,  appears  that  the  two-dimensional  source  term  produces  an 
unrealistically  large  effect.  The  source  effect  calculated  by  slender- 
body  theory  is  not  accurate  for  reasons  subsequently  discussed,  and  the 
theory  is  modified  accordingly. 

The  reason  that  the  store  source  terms  produce  a large  effect  lies 
in  a known  limitation  of  slender-body  theory.  The  term  in  the  complex 
potential  for  the  body  source  varies  as  log  r and  does  not  converge  for 
distances  far  fro.u  the  body  in  contrast  to  the  doublet  solution  associated 
with  body  lift.  However,  an  excellent  approximation  to  the  body  volume 
effects  on  additional  interference  can  be  obtained  very  simply  by  using 
three-dimensional  sources  to  represent  the  body  and  considering  the  wing  to 
be  a reflection  plane.  Alternately  we  could  consider  the  modification  of  the 
wing  vortex-lattice  solution  due  to  the  effect  of  the  three-dimensional 
body  sources  on  the  wing  boundary  conditions.  The  effect  of  store  volume  on 
the  additional  interference  has  been  calculated  both  ways  and  compared. 

Both  methods  gave  closely  similar  results  for  the  present  case  and  were 
negligible  compared  to  the  primary  interference.  Accordingly,  the  body 
volume  effects  will  be  calculated  using  three-dimensional  source  distri- 
butions. Except  for  stores  very  close  to  the  wing  they  will  be  negligible. 

Figure  30  presents  comparisons  between  the  primary  interference  theory 
and  data  for  the  same  case  as  that  of  figure  28  except  that  the  Mach 
number  is  0.70.  As  in  the  Mach  number  0.25  case,  the  calculations  have 
been  made  taking  into  account  the  indicated  tunnel-empty  stream  angle. 

The  overall  agreement  between  experiment  and  theory  is  not  quite  as  aood 
as  that  shown  in  figure  28  for  a Mach  number  of  0.25.  This  is  consistent 
with  the  results  of  the  flow-field  survey  comparisons  of  figures  18  and  19 
where  it  was  found  that  better  agreement  was  obtained  for  a Mach  number  of 
0.25.  The  effects  of  tunnel  stream  angle  on  these  comparisons  should  be 
recalled. 
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The  primary  interference  theory  is  compared  with  the  tota]  norma.l- 
force  and  pitching-moment  coefficients  as  a function  of  angle  of  attack 
in  figure  31.  This  is  the  same  configuration  and  store  location  as  was 
considered  in  the  previous  figures.  Figure  31(a)  is  for  a Mach  number  of 
0.25  and  figure  31(b)  is  for  a Mach  number  of  0.70.  At  both  Mach  numbers 
the  calculated  normal  force  is  slightly  high.  This  is  to  be  expected, 
since  in  both  figures  28  and  30  the  calculated  load  distributions  were 
generally  above  the  experimental  ones.  The  calculated  pitching  moment  is 
in  good  agreement  with  the  data.  Since  the  overprediction  of  the  store 
load  distribution  is  generally  around  the  store  midpoint,  the  point  about 
which  the  moment  is  taken,  this  overprediction  has  little  effect  on  the 
moment. 

8. 2 Single  Store  Under  Wing  of  Wing-Fuselage- Pylon  Combination 

Let  us  now  consider  the  case  of  a single  store  below  a pylcn 
mounted  at  the  2/3  semi span  position  on  the  left  wing  panel  of  the 
wing-fuselage  combination  of  figure  2.  The  details  of  the  pylon  are 
shown  in  figure  32.  When  mounted  on  the  wing,  the  pylon  centerline  is 
located  at  40  percent  of  the  local  wing  chord.  The  store  location  if  it 
were  mounted  on  the  pylon  is  shown  in  figure  3 (a) . The  store  location  to 
be  considered,  however,  is  one  diameter,  0.75  inch,  below  this  location, 
the  same  position  considered  in  the  previous  section.  The  store  details 
are  shown  in  figure  3(b).  With  the  pylon  present  the  store  midpoint  is 
directJy  below  the  pylon  centerline. 

In  Volume  III  of  reference  1,  load  distributions  art}  presented  which 
were  measured  both  with  and  without  the  pylon  present  on  the  wing.  By 
subtracting  the  data  obtained  without  the  pylon  present  from  thai  obtained 
with  the  pylon  on  the  wing,  the  incremental  load  distribution  due  to  tne 
addition  of  the  pylon  can  be  obtained.  This  has  been  done  and  the  incre- 
mental loadings  are  plotted  xn  figure  33.  Data  are  shown  for  two  angles  of 
attack  at  a Mach  number  of  0.25.  It  can  be  seen  from  the  two  parts  of  the- 
figure  that  the  effect  of  angle  of  attack  is  small.  Also,  there  is  quite 
a large  upstream  influence  of  the  pylon  which  the  slender-body  theory  method 
of  accounting  foe  the  pylon,  section  7,  will  not  predict, 
o 

For  a p = 6 , figure  33(b),  an  additional  interference  calculation 
has  bean  made  to  see  what  portion  of  the  experimental  incremental  load 
distribution  is  due  to  adding  a pylon  of  zero  thickness.  The 
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result  of  the  calculation  is  shown  in  the  figure.  As  can  be  seen,  only 
a small  portion  is  accounted  for.  It  is  felt  that  the  primary  effect 
showing  up  in  the  data  is  one  of  pylon  thickness  which  is  not  accounted 
for  in  either  the  primary  or  additional  interference. 

Even  though  the  pylon  thickness  effect  changes  the  store  load  distri- 
bution, its  effect  on  C„  and  C is  small.  If  the  curves  of  figure  33 
are  integrated  over  the  store  length,  the  resulting  ACN*s  are  approximately 
equal  to  the  2og  (two  standard  deviation)  uncertainty  of  the  data  which  was 
determined  in  Volume  I of  reference  1.  Since  the  curves  of  figure  33  are 
fairly  symmetrical  about  the  store  midpoint,  the  point  about  which  the 

moments  are  to  be  taken,  small  values  of  AC  are  also  obtained.  These 

m 

results  may  not  be  true  in  general  even  though  the  present  pylon  is 
excessively  thick.  It  thus  may  be  prudent  to  include  the  pylon  thickness 
in  the  present  prediction  technique. 

Another  set  of  data  for  which  comparisons  between  the  present  theory 
and  experiment  va1!  be  made  is  contained  in  reference  15.  The  configura- 
tion is  shown  in  figure  34.  This  is  the  swept-wing  model  tested  in 
reference  15.  The  store  is  located  at  the  one- third  sem.^pan  approximately 
18  percent  of  the  local  chord  below  the  wing  and  is  attached  to  the  pylon. 
Comparisons  with  the  measured  normal  force  and  pitching  moment  will  be 
presented  for  both  tail-on  and  tail-off  stores.  McKinnev  and  Polhamus, 
in  reference  10,  also  have  made  calculations  to  compare  with  the  same  data. 
Their  calculations  were  made  using  the  flow  field  data  of  reference  5,  which 
were  taken  under  a slightly  different  wing-body  combination,  to  determine 
the  nonuniform  flow  field  in  the  vicinity  of  the  store.  In  the  present 
work  the  flow  field  has  been  calcxilated  by  the  methods  described  previously. 
Comparisons  with  their  calculations  will  be  presented,  not  onry  £or  the 
total  normal  force  anu  pitching  moment,  but  also  for  their  components  due 
to  buoyancy  and  angle  of  attack. 

The  calculations  tc  be  presented  are  primary  interference  calculations 
only.  As  was  discussed  in  section  7.5.2,  in  connection  with  the  conical 
configuration,  the  additional  interference  analysis  is  not  applicable  to  the 
case  of  zero  yap  between  the  stone  and  pylon.  In  the  present  case,  a very 
small  gap  exists  since  the  store  was  tested  in  the  carriage  position  with 
the  balance  supported  through  the  pylon. 
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Before  presenting  the  normal  force  and  pitching  moment  determined  by 
i itegrating  over  the  store  length,  it  is  of  interest  to  examine  the  magni- 
tudes and  distributions  of  the  various  components  of  the  normal  force  over 
the  length  of  the  store  shown  in  figure  34.  The  distributions  to  be  shown 
are  in  figure  35  for  the  tail-off  model  at  four  angles  of  attack.  At  the 
top  of  each  figure  is  a sketch  of  the  wing-pylon-store  cross  section  at 
the  one-third  semispan  under  the  left  wing,  the  semispan  location  where 
the  calculations  were  made.  For  the  purposes  of  this  comparison  only, 
the  separation  location  has  been  assumed  to  be  at  60  percent  of  the  store 
length,  xs  -/ls  = 0.60.  This  region  is  shown  by  the  shaded  area  on  the 
store.  This  store,  unlike  the  one  for  which  the  previous  comparisons  have 
been  made  is  boattailed  so  that  flow  separation  over  the  aft  end  of  the 
store  is  to  be  expected.  No  method  exists  for  determining  the  separation 
location  on  a body  of  revolution  in  a nonuniform  flow  field.  The  method 
of  Hopkins,  reference  12,  which  was  mentioned  earlier,  applies  only  to  a 
uniform  flow  field. 

An  examination  of  figure  35  shows  that  both  the  buoyancy  and  slender- 
body  theory  load  distributions  have  regions  where  the  leading  is  negative 
and  regions  where  it  is  positive.  The  buoyancy  distribution  extends  to 
the  end  of  the  body  whereas  the  slender-body  theory  component  ends  at 
the  assumed  separation  location.  The  component  due  to  viscous  crossflow 
theory  acts  only  downstream  of  this  point. 

''or  all  angles  of  attack,  the  buoyancy  force  distribution  exhibits 
the  same  qualitative  behavior.  The  loading  near  the  wing  leading  edge 
is  negative  and  becomes  positive  under  the  mid-portion  of  the  wing. 

It  oecomes  negative  again  near  the  wing  trailing  edge.  The  sign  of  this 
component  of  the  loading  changes  as  the  sign  of  tne  pressure  gradient 
across  the  store  changes. 

The  slender-body  theory  loading  curves  are  similar  in  behavior  to 
the.  buoyancy  curves  except  aheaci  of  the  wing  leading  edge  where,  at 
large  angles  of  attack,  a large  positive  loading  exists.  The  behavior 
of  the  slender-body  curves  can  be  explained  by  examining  the  equation 
for  the  .loading.  From  equation  {50) 
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The  sign  of  the  first  term  depends  on  the  sign  of  W*  as  well  as  whether 
the  body  is  expanding  or  contracting.  The  sign  of  the  second  term  depends 
on  the  sign  of  the  axial  gradient  of  W*.  Thus,  negative  loadings  due  to 
slender-body  theory  can  exist  in  regions  where  the  body  is  expanding  i f 

(1)  W*  and  dW*/d{x/i)_  are  both  negative. 

s s s 

(2)  W*  is  negative,  dW*/d(x/£)„  is  positive,  but  term  1 is  larger 

S S *3 

in  magnitude  than  term  2. 

(3)  W*  is  positive,  dw*/d(x/i)g  is  negative,  and  term  2 is  larger 
in  magnitude  than  term  1. 

It  is  also  possible  to  have  nonzero  loading  on  the  cylindrical  portion  of 
the  store  because  of  the  velocity  gradient  term,  term  2.  The  store  shown 
in  figure  35  is  cylindrical  from  (x/f)  s = 0.315  to  (x/i)g  = 0.617. 

The  viscous  crossflow  theory  loadings  are  small  for  the  angles  of . 
attack  shown  since  from  equation  (65)  it  is  proportional  to  The 

product  V*W*  is  proportional  to  the  square  of  the  local  angle  of  attack. 

Let  us  now  examine  the  normal  force  and  pitching  moment  obtained  by 
integrating  the  load  distributions  presented  in  figure  35  plus  those 
calculated  for  two  other  angles  of  attack,  -2°  and  6°.  The  results  are 
presented  in  figure  36.  The  nondiir.ensionalizing  length  used  in  the 
pitching  moment  shown  in  figures  36  and  37  is  the  store  length,  not  the 
maximum  diameter  used  elsewhere  in  the  report.  The  buoyancy  component  is 
shown  in  figure  36(a)  with  the  results  of  McKinney  and  Polhamus  (ref.  10) 
also  shown,  hs  was  mentioned  previously,  their  calculations  were  made 
using  experimental  downwash  data  from  reference  5.  These  data  were  taken 
under  a slightly  different  wing-body  combination,  see  figure  9,  from 
that  shown  in  figure  34.  The  present  calculations  are  entirely  theoretical. 
The  two  calculations  are  in  good  agreement  for  both  normal  force  and 
pitching  moment. 

A similar  comparison  is  presented  in  figure  36(b)  for  the  angle  of 
attack  component,  that  is,  the  sum  of  slender-body  theory  and  viscous 
crossflow  theory.  Calculations  using  the  present  method  are  shown  for 
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three  assumed  locations  of  separation.  The  results  show  the  importance 
of  knowing  the  separation  point  location  in  calculating  the  normal  force, 
particularly  at  small  angles  of  attack.  It  is  not  nearly  as  important 
in  calculating  the  pitching  moment  since  the  region  over  which  the  norma] 
force  is  affected  by  shifting  the  separation  point  is  near  the  center  of 
moments.  For  a particular  value  of  the  separation  location,  the  present 
normal-force  calculation  coincides  with  the  calculation  of  reference  10. 

It  is  not  stated  where  separation  was  assumed  in  the  calculations  of 
rfer^nce  10.  The  pitching-moment  calculations,  however,  do  not  agree. 

The  reason  for  this  is  not  understood.  The  presence  of  the  pylon  was  also 
neglected  in  reference  10  since  the  experimental  flow-field  data  used  there 
were  obtained  under  a wing-body  combination  without  a pylon. 

The  total  normal  force  and  pitching  moment,  that  is,  the  sums  of  the 
buoyancy  and  angle  of  attack  components,  are  shown  in  figure  36(c).  Also 
shown  are  the  calculations  of  reference  10  and  the  data  from  reference  15. 
For  the  intermediate  assumed  separation  location,  x /f  equal  to  0.47, 

S y O S 

the  present  calculations  agree  fairly  well  with  the  experimental  data. 

The  comparisons  between  the  total  predicted  and  measured  normal  force 
and  pitching  moment  for  the  external  store  with  fins  are  shown  in  figure  37. 
Again,  the  theory'  shows  an  effect  due  to  location  of  separation  on  the 
calculated  normal  force.  In  fact,  these  differences  are  precisely  those 
shown  in  figure  36(c)  for  the  store  wi thout  fins  since  the  calculated 
contribution  of  the  fins  does  not  depend  on  position  of  separation.  For 
the  intermediate  position  of  separation,  the  normal-force  theory  and 
experiment  are  separated  by  a small  increment  up  to  about  6°  angle  of 
attack. 

Above  this  angle  of  atcack,  the  method  overpredicts  the  tail  contri- 
bution by  an  amount  that  increases  as  the  angle  of  attack  increases.  In 
other  words,  the  tail  effectiveness  decreases  as  the  angle  of  attack 
increases.  This  nonlinear  effect  could  well  be  the  result  of  operating 
the  tail  in  the  presence  of  body  vortices  associated  with  viscous  cross- 
flow.  Methods  for  calculating  this  effect  are  known  for  uniform  flow. 

For  the  present  case  of  a store  in  a nonuniform  flow,  an  approximate 
method  should  probably  be  developed  to  account  for  body  vortex-tail 
inter ferencc . Methods  for  including  such  interference  are  available  in 
ref ere "C 
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The  comparison  between  total  moment  coefficient  as  calculated  and  a3 
measured  for  the  store  with  fins  is  in  good  agreement  up  to  about  6°  angle 
of  attack.  For  higher  angles,  there  is  less  nose-down  moment  developed 
experimentally  than  predicted.  This  result  is  consistent  with  the  fact 
that  the  tail  normal  force  has  been  overestimated  by  a theory  which 
neglects  body  vortex  effects  on  tail  effectiveness. 

8.3  TER  Stores  Unoer  Wing  of  Wing-Fuselage-Pylon- Rack  Combination 

An  illuminating  comparison  can  be  made  between  the  primary  interference 
theory  of  section  6 and  the  data  of  Volume  III  of  reference  1 for  the  incre- 
mental load  distribution  due  to  any  stores  remaining  on  the  TER  rack. 
Consider  the  TER  grouping  as  shown  in  figure  20  and  let  us  determine  the 
incremental  load  distribution  on  store  number  1,  one  store  diameter  below 
the  position  shown  in  the  figure,  due  to  the  addition  of  stores  2 and  3 
to  the  rack.  Data  are  presented  in  Volume  III  of  reference  1 with  and 
without  these  stores  present  so  that  a subtraction  of  the  two  load  distri- 
butions give  the  incremental  loading  due  to  the  stores.  This  has  been 
done,  and  the  data  are  plotted  in  figure  38  for  two  angles  of  attack. 

Note  that  there  is  not  a large  change  with  angle  of  attack. 

Also  shown  on  the  two  parts  of  the  figure  are  curves  calculated 
using  the  primary  interference  theory  of  section  6.  These  curves  were 
obtained  by  first  calculating  the  loading  produced  by  the  velocity  field 
due  to  fuselage  volume,  wing  thickness,  wing  vorticity  distribution 
including  store-induced  camber  due  to  all  three  stores,  and  the  volume 
distributions  of  stores  2 and  3.  The  calculation  was  then  repeated 
excluding  stores  2 and  3,  including  the  wing  camber  induced  by  them,  and 
the  two  load  distributions  subtracted.  The  resulting  incremental  load 
distribution  is  nearly  independent  of  angle  of  attack  since  it  is  due 
primarily  to  the  velocities  induced  by  the  volume  distributions  of 
stores  2 and  3,  which  are  not  functions  of  angle  of  attack.  It  can  be 
seen  from  figure  36  that  the  incremental  load  distribution  is  predicted 
quite  well  at  both  angles  of  attack.  The  dominant  store-store  inter- 
ference effect  is  clearly  due  to  body  sources. 

The  largest  differences  between  experiment  and  theory  occur  at  the 
aft  end  of  the  store.  This  region  is  influenced  by  the  shape  used  to 
model  the  wakes  of  stores  2 and  3.  In  the  present  calculations,  the 
wake  was  assumed  to  have  the  same  shape  as  the  store  nose;  that  is,  it 
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was  taken  to  be  an  ogive  0.125  foot  long  (see  the  store  details  in 
figure  3(b)).  If  it  had  been  taken  more  slender,  then  the  values  of 
f'dC^/dx,.  near  the  store  base  would  have  been  reduced.  For  blunt-based 
stores  in  TER  or  MER  grouping,  the  modeling  of  the  wake  is,  therefore, 
important. 

Similar  comparisons  to  those  of  figure  38  are  presented  in  figure  39 
except  that  the  incremental  load  distribution  on  store  2 due  to  store  3 
on  the  TER  rack  are  shown.  Store  1 is  not  present  and  store  2 is  one 
diameter,  0.75  inch,  below  the  position  shown  in  figure  20.  Again,  the 
experimental  distribution  does  not  change  much  with  angle  of  attack  and 
is  predicted  quite  well. 

9.  EXPERIMENTALLY  DETERMINED  MER  INTERFERENCE  EFFECTS 

It  is  dear  that  the  interference  effects  described  in  the  previous 
section  for  TE»  configurations  also  exist  as  a large  effect  in  MER 
configurations.  However,  there  is  in  addition  the  interference  between 
the  frosst  three  stores  and  the  rear  three  stores.  Data  are  presented  in 
Volume  V of  reference  1 in  which  the  various  interference  effects  asso- 
ciated with  a multiple  ejection  rack  (MER)  installation  have  been  isolated. 
This  was  accomplished  by  measuring  forces  and  moments  on  a store  in  a 
fixed  position  relative  to  the  wing-fuselage  combination  and  building  the 
remainder  of  the  configuration  up  component  by  component.  The  wing- 
fuselage  combination  used  in  the  tests  is  shown  in  figure  2.  The  stores 
are  located  as  shown  in  figure  40(a)  when  in  their  carriage  position  on 
the  rack.  The  details  of  the  stores  are  shown  in  figure  40(b).  The 
pylon  details  are  shown  in  figure  32  and  the  MER  rack  details  in 
figure  41. 

On  the  following  flares,  abbreviations  have  been  used  to  designate 
the  various  aircraft  coagsments.  They  are 

WB  wing- fuselage 

P psrlon 

M MSS'  rack 
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The  numbering  system  used  to  identify  the  six  stores  is  shown  in  the 
following  sketch  as  is  the  fin  orientation. 


The  various  stores  will  be  designated  as  Sx , S2,  etc. 

Data  taken  using  f inless  stores  are  presented  in  figure  42.  The 
f inless  stores  used  are  identical  to  the  finned  one  shown  in  figure  40(b) 
except  that  the  fins  have  been  removed.  Increments  in  normal  force  ACN 
and  pitching  moment  ACm  are  plotted  against  angle  of  attack  for  a Mach 
number  of  0.40.  In  figure  42(a),  data  taken  on  MER  store  number  1 are 
presented.  Similar  data  are  presented  in  figure  42(b)  for  store  number  4. 
These  stores  were  located  0.375  inch  below  the  position  shown  in  figure  40(a) 
when  data  were  taken. 

Let  us  first  consider  the  store  number  1 data  shown  in  figure  42(a). 

The  increments  plotted  were  obtained  by  differencing  the  data  obtained 
with  two  parent  aircraft  configurations.  Thus,  the  circles  indicate  the 
effect  caused  by  adding  the  wing- fuselage  combination  to  the  flow  field. 
Similarly,  the  squares  show  the  change  caused  by  adding  the  pylon;  the 
diamonds,  the  MER  rack;  the  equilateral  triangles,  stores  S2  and  S3; 
and  the  right  triangles,  stores  S5  and  S6.  It  can  be  seen  that  the 
largest  interference  effects  are  caused  by  the  wing-fuselage  and  stores 
S2  and  S3.  The  smallest  change  in  force  and  moment  is  caused  by  the 
pylon.  Since  the  wing  and  fuselage  tend  to  align  the  flow  with  the 
fuselage,  the  store  sees  a smaller  angle  of  attack  than  when  it  is  in  the 
stream  by  itself.  Thus,  the  normal  force  on  the  store  and  the  resulting 
pitching  moment  are  reduced.  The  two  stores  S2  and  S3  produce  an 
effect  similar  to  that  seen  in  figure  38  for  the  TER  racks.  This  is  to 
be  expected,  since  the  MER  rack  is  really  two  TER  racks  arranged  one 
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behind  the  other.  If  the  curves  of  figure  38  were  integrated,  we  would 
find  a positive  ACN  but  because  of  the  negative  loading  at  the  nose  the 
moment  about  the  store  midpoint  would  be  negative.  The  addition  of  stores 
S5  and  S6  produce  a downwash  over  the  aft  end  of  Sx  since  they  are 
above  and  behind  this  store.  This  causes  a negative  AC^  and  because  this 
force  is  primarily  generated  on  the  aft  end  of  the  store  a positive  or 
nose-up  pitching  moment  is  produced. 

Data  taken  on  store  number  4 are  shown  in  figure  42 (b) . Again  it 
is  seen  that  the  wing-fuselage  and  the  two  side  stores  in  the  same 
grouping,  in  this  case  Ss  and  S6,  have  the  largest  effects.  The  front 
three  stores  produce  no  change  in  normal  force;  however,  they  do  produce 
a negative  increment  in  pitching  moment. 

Let  us  compare  the  effects  of  S5  and  S6  on  Sx  with  the  effects 
of  S2  and  S3  on  S4.  The  stores  are  nearly  symmetrical  fore  and  aft 
so  that  the  source  distributions  were  nearly  asymmetric  fore  and 
aft.  Accordingly,  the  downwash  produced  at  a point  at  a given  distance 
in  front  of  the  trailing  edge  of  store  by  S5  and  S6  will  be  equal 

and  opposite  to  that  produced  by  S2  and  S3  the  same  distance  behind  the 
leading  edge  of  store  S4.  Such  distributions  will  produce  equal  normal- 
force  increments  in  the  two  curves,  but  opposite  moment  increments.  Also 
the  increments  should  not  vary  with  angle  of  attack.  The  measured  data 
exhibit  these  characteristics  fairly  well  when  consideration  is  given  to 
the  fact  that  the  stores  are  not  precisely  symmetrical  fore  and  aft,  and 
in  one  case  we  are  looking  at  the  effect  of  adding  Sc  and  Sg,  and  in  the 
other  Sx,  S2,  and  S3.  Also,  any  effects  of  the  wake  entrainment  would 
be  present  in  one  case  but  not  the  other.  Based  on  these  observations, 
it  appears  that  the  induced  effect  of  one  set  of  these  stores  on  one  of 
the  other  trio  can  be  explained  by  three-dimensional  source  distributions 
neglecting  wake  entrainment. 

In  figure  43  data  taken  on  the  finned  store  shown  in  figure  40 (b)  are 
plotted  against  angle  of  attack.  Data  taken  on  store  number  1 are  presented 
in  figure  43(a)  and  those  taken  on  store  number  4 in  figure  43(b).  The 
stores  are  located  in  the  same  positions  as  they  were  in  the  previous 
figure,  that  is,  0.375  inch  below  the  locations  shown  in  figure  40(a). 

For  both  stores  Sx  end  S4,  the  addition  of  the  wing-fuselage- 
pylon-MER  configuration  to  tre  flow  field  adds  a negative  increment  in 
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normal  force  over  most  of  the  anglo-of-attack  range.  As  in  the  case  of 
the  f inless  store,  this  is  caused  by  the  fuselage  and  wing  tending  to 
align  the  flow  with  the  fuselage  axis  and  thus  reduce  the  angle,  of 
attack  which  the  store  sees.  At  small  angles  of  attack  a stabilizing 
moment  increment  is  added  while  at  large  angles  a destabilizing  incre- 
ment is  added.  This  later  effect  is  caused  by  the  tail  fins  operating 
at  a reduced  angle  of  attack. 

Adding  the  two  side  stores  S£  and  S3  above  store  Sx,  (fig.  43(a)) 
or  S2  and  Se  above  S4  (fig.  43(b))  produces  the  same  effect,  a 
positive  ACN  and  a negative  ACm<  The  magnitudes  of  the  effects  are 
nearly  equal.  This  i3  the  behavior  expected  from  a TER  grouping  as  was 
discussed  in  connection  with  figure  42. 

The  effect  of  the  two  stores  S5  and  Ss  on  store  Sx,  figure  43(a), 
is  to  produce  a small  negative  normal-force  increment  and  a positive 
pitching-moment  increment  almost  independent  of  angle  of  attack.  Since 
stores  S5  and  Ss  produce  a downwash  at  the  tail  of  Sx , the  local 
angle  of  attack  is  reduced  which  tends  to  reduce  the  stabilizing  moment 
contributed  by  the  tail.  The  effect  of  the  front  three  stores,  Sx,  S2, 
and  S3,  on  S4  can  be  seen  in  figure  43(b).  There  is  little  effect  on 
normal  force  and  a small  stabilizing  moment  increment.  This  moment 
effect;  is  not  caused  by  the  presence  of  the  tail  fins  of  store  4.  A 
comparison  of  the  curve  indicated  by  the  diamonds  in  figure  43(b)  with 
the  comparable  curve  for  the  finless  store  in  figure  42 (b) , that  given  by 
the  right  triangles,  shows  the  AC^'s  to  be  almost  equal.  Note  the  change 
of  scales  between  the  two  figures. 

In  the  computer  program,  the  effect  of  store  volume  has  been  included, 
as  it  influences  store-store  interference  in  the  TER  configurations  and 
the  MER  configurations. 


10.  STORE  TRAJECTORIES 

In  the  preceding  sections  of  this  report,  methods  have  been  presented 
which  allow  the  calculation  of  the  normal  force  and  pitching  moment  acting 
on  a store  at  any  point  in  the  aircraft  flow  field.  Utilizing  these  tools 
a computer  program  has  been  written  which  will  integrate  the  equation  of 
motion  of  the  store  to  determine  its  location  and  angular  orientation  as 
a function  of  time.  For  the  present  these  equations  are  restricted  to 
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three  degrees  of  freedom  representing  .lotion  in  a vertical  plane  parallel 
to  the  airplane  vertical  plane  of  symmetry.  The  airplane  is  assumed  to 
bo  in  a straight  uniform  flight  but  may  be  climbing  or  diving. 

In  this  section  of  the  report,  the  equations  of  motion  will  first 
be  derived.  Then  the  computer  program  will  briefly  be  described  and, 
finally,  a sample  trajectory  will  be  presented. 

10.1  Equations  of  Motion 

The  equations  of  motion  will  be  written  in  a coordinate  system  fixed 
in  the  aircraft  fuselage.  This  system  will  have  its  origin  at  the  fuselage 
nose  and  will  move  with  the  fuselage.  The  coordinate  system  is  shown  in 
figure  44  as  the  £,t|  system.  The  aircraft  is  flying  at  constant  free- 
stream  velocity,  Voo,  and  at  a constant  flight  path  angle,  relative 
to  the  horizontal.  The  angle  of  attack  of  the  aircraft,  is  also 
constant „ At  time  t = t the  store  is  located  as  shown  in  figure  44. 

Its  center  of  gravity  is  located  at  ^ 0 » "Hq  and  its  a*is  is  inclined  at 
some  angle,  0Q,  measured  with  respect  to  a line  parallel  to  the  fuselage 
axis.  Provision  will  be  made  for  imposing  an  initial  velocity  relative 
to  the  aircraft,  V , and  an  initial  angular  velocity,  9^. 

The  equations  of  motion  of  the  store  for  the  nonaccelerating  coordinate 
system  are: 


L 

2 


m_ 


Pco^a,  ^R 

S 


F. 


2 P«v»  SR 


m 


s 


1 

2 


I c» 


R 


1 

2 


1 

2 


m k' 

s 

• >.V^  SR^R 
s 


M, 


7 Pa>Vou  SR/R 

S 


(122) 


(123) 


(124) 


The  forces  acting  on  the  store  are  shown  in  figure  43.  NTs  is  the 
normal  force  acting  on  the  body  and  is  the  sum  of  three  components 
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Ns  “ nby  + nsb  + ncf 


(125) 


where 


NBY  " norma^  f°rce  due  to  buoyancy 

Neu  « normal  force  due  to  slender-body  theory 

N ._  - normal  force  due  to  viscous  crossflow 
cF 

For  stores  with  empennages,  an  additional  normal  force  is  developed*  This 
is  shown  in  figure  45  as  N^,.  These  forces  are  the  primary  interference 
forces  and  tha  methods  of  evaluating  them  and  their  corresponding  moments 
were  presented  in  section  6.  From  here  on,  the  total  normal  force, 
including  the  empennage  force  if  an  empennage  is  present,  will  be 
designated  N and  the  corresponding  moment,  M. 

The  remaining  two  forces  shown  in  figure  45  are  the  zero-lift  drag 
force,  Dq,  ana  the  gravitation  force,  mgg.  The  zero-lift  drag  force  will 
be  taken  to  always  act  in  the  V-  direction  4-  the  store  center  of 

b 

gravity. 

Referring  to  figures  44  and  45,  the  forces  can  be  resolved  into  their 
components  in  the  £ a°d  T directions  to  obtain 

- N sin  6 + DQ  cos  (aB  + YB  ~ Ys) 


+ mgg  sin  (aB  + yQ) 


(126) 


F^  = N cos  S + Dg  sin  (a0  + YB  “ 7S) 


- mgg  cos  (a0  + y0) 


(127) 


Thus,  equations  (122)  and  (123)  can  be  rewritten  as 

S„  r 


1 "2  “ [CN  Sin  6 + CD,o  COS  (aB  + yB-  Ys}j 


^ ” T P°ovoo  m 
s s 


+ g sin  (a  + yK) 


(12  0) 
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and 


* ■ * P-s  S K 


COS  e + CD,0  Sin  (aR  +Vn- 


B ' B 


ys) 


- g cos  (aQ  + <yB) 
Equation  (124)  can  be  written 


(129) 
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1 v2  SR^R  r 

T Poo  co  m k2  Cm 

s s 


(130) 


Equations  (128) , (129) , and  (130)  are  a set  of  three  second-order 
differential  equations  which  must  be  solved  to  determine  the  store  motion 
relative  to  tht  moving  aircraft.  The  six  initial  conditions  are  at 
t = t 


5 " So  > n - % » e - er 


(131) 


t = e 


% ’ 9 • 0o 


The  values  of  £o  and  r)Q  are  determined  by  specifying  the  location  of 
the  store  center  of  gravity. 

The  initial  pitch  angle  is  obtained  as  follows 


Aw  + Ls 


(132) 


wherein  the  wing  incidence  relative  to  the  fuselage  axis,  i , and  the 
store  incidence  relative  to  the  wing  root  chord,  i.s,  are  also  input  data. 

The  initial  conditions  on  the  velocities  are  to  be  specified  as  an 
initial  rotational  velocity,  9Q,  positive  in  the  nose-up  direction,  in 


radians  per  second,  and  an  initial  velocity,  v , in  a direction  normal 


to  the  store  axis  and  positive  as  shown  in  figure  44.  Then 


t = 
"O 


-V  sin  G 
o o 


(133) 


% = 


-v  cos 
o o 
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in  the  equations  of  motion  certain  other  quantities,  which  are 
constant  with  time,  are  required.  They  are 


a „ = maximum  radius  of  the  store 

max 


g = gravitational  acceleration 


k = store  radius  of  gyration 
fn  = reference  length,  taken  as  2a 


mg  = mass  cf  the  store 


V = aircraft  flight  velocity 


aQ  * the  fuselage  angle  of  attack 
YB  = the  fuselage  flight  path  angle 
p = the  free-stream  mass  density 


In  addition,  the  reference  area,  S_,  is  needed.  This  is  taken  as 

K 

the  maximum  cross-sectional  area  of  the  store 


SR  = 7ramax 


(134) 


The  last  quantity  which  is  not  a function  of  time  and  which  appears  in 
the  equations  of  motion  is  the  store  drag  coefficient,  Q.  This  must 
be  specified  and  is  defined  as 


U35) 


2 PooVro  SR 
s 


The  remaining  quantities  in  the  equations  of  motion,  V,  > 7S>  and 
the  normal-force  and  pitching-moment  coefficients,  are  functions  of  time. 
The  store  velocity  . slative  to  the  atmosphere  is 


+ (-Vw  sin  aB  + *j)  j 


(136) 
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The  flight  path  angle  of  the  store,  ys , is  found  from  the  vertical  and 
horizontal  components  of  its  velocity  relative  to  the  atmosphere. 
Consider  the  following  sketch  (also  see  figures  44  and  45) . 


The  coordinate  system  is  moving  at  a velocity  VM.  The  store  is 

moving  in  this  c,r)  coordinate  system  a"-  velocities  \ and  tj.  Thus, 
the  horizontal  and  vertical  velocities  of  the  store  relative  to  the 
atmosphere  are 


Vhor  s V«>  cos  > - 1 COS  + V " i rin  (^B  + ^ 

Vvert  “ V«  Sin  ' * Sin  % 'h  aB>  + ^ C°S  ‘’VB  + aB} 

The  flight  path  angle,  7 in  figure  45,  is  then 


7 = tan 

' s 


■ i f Vvert  \ 
\Vhor 


(137) 


1 138) 


The  normal-force  and  pitching-moment  coefficients  are  calculates  the 
methods  presented  in  section  6. 
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10.2  Description  of  Computer  Program 

The  computer  program  which  calculates  the  trajectory  of  a store 
ejected  from  an  aircraft  actually  consists  of  three  separate  programs  in 
its  present  form.  The  three  programs  are 

(1)  Axisymmetric  and  two-dimensional  source  distribution  program 

(2)  Vortex-lattice  program 

(3)  Trajectory  program 

For  a given  body  or  airfoil  and  a given  Mach  number,  the  first  program 
needs  to  be  run  only  once.  For  a given  airplane-store  configuration  and 
Mach  number,  the  second  program  needs  to  be  run  only  once.  The  third 
program  cam  be  run  for  various  angles  of  attack,  altitudes,  and  store 
initial  dynamical  conditions  without  rerunning  the  first  two  programs. 

The  extent  to  which  the  three  programs  should  be  incorporated  into  one 
program  depends  therefore  on  the  particular  problem  at  hand. 

The  first  program  calculates  the  source  distributions  which  repre- 
sent the  fuselage  and  the  stores,  present  on  the  aircraft,  including  the 
one  to  be  ejected.  Each  shape  is  specified  by  a series  of  segmented 
polynomials  and  the  program  calculates  and  outputs  the  axial  locations 
of  the  sources  and  their  strengths. 

This  program  also  calculates  the  two-dimensional  source  distribution 
which  represents  the  streamwise  airfoil  section  of  the  wing  at  the  span- 
wise  station  immediately  above  the  store  to  be  ejected.  The  airfoil 
section  is  also  specified  by  a series  of  segmented  polynomials,  and  the 
program  calculates  and  outputs  the  source  strengths  and  locations. 

The  second  program,  the  vortex  lattice  program,  computes  the  bound 
vorticity  distributions  which  represent  the  wing  camber  and  twist  and 
angle  of  attack.  The  input  data  to  this  program  consist  of  wing  geometry 
and  information  which  allows  the  program  to  locate  the  fuselage  and  all 
of  the  stores  relative  to  the  wing.  In  its  present  form,  up  to  nine 
stores  can  be  placed  under  the  fuselage  centerline  and  one  wing  panel. 

In  addition  to  the  above  data,  the  source  distributions  representing  the 
fuselage  and  stores  are  input  as  is  the  geometric  camber  and  twist  distri- 
bution of  the  wing.  The  first  calculation  performed  by  this  program  is  to 
determine  the  velocities  induced  normal  to  the  wing  by  the  fuselage  and 
the  stores.  These  velocities  arc  treated  as  an  induced  camber,  and  this 
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camber  distribution  is  added  to  the  wing  geometr.c  twist  and  camber. 

The  program  then  solves  for  the  two  vorticity  distributions.  The  first 
cancels  the  free-stream  velocity  component  normal  to  the  wing  surface 
and  the  second  cancels  the  twist  and  camber  distribution.  These  vorticity 
distributions  are  output  by  the  program. 

The  third  program  is  the  trajectory  program.  The  input  data  to  this 
program  consist  of  all  of  the  data  read  into  the  second  program,  the 
vorticity  distribution  output  by  the  second  program,  and  other  information 
required  for  the  force  and  moment  and  trajectory  calculations.  The 
additional  information  input  includes  the  store  mass,  radius  of  gyration, 
center  of  gravity  location,  location  of  the  tail,  tail  semispan,  and 
tail  lift-curve  slope.  Also,  the  aircraft  flight  conditions  are  input 
as  are  the  store  downward  ejection  velocity  and  the  ejection  pitching 
rate. 

The  trajectory  calculation  is  begun  with  the  store  at  a specified 
initial  position.  The  velocity  field  is  calculated  over  the  store  length  , 
and  using  this  field  the  forces  and  moments  are  calculated.  The  velocity 
field  is  found  by  adding  to  the  free-stream  velocity  the  velocities  induced 
by  the  fuselage,  wing,  and  all  other  stores.  The  integration  of  the 
equations  of  motion  is  accomplished  by  a standard  numerical  integration 
technique  with  the  velocity  field  and  the  forces  and  moments  being 
recalculated  at  each  point  required  by  the  integration  scheme. 

10.3  Sample  Trajectory 

An  example  of  the  application  of  the  trajectory  program  is  shown  in 
figure  46.  This  figure  presents  a comparison  of  two  calculated  trajec- 
tories with  an  experimental  trajectory  taken  from  reference  17  which  was 
obtained  using  the  captive  trajectory  system  in  the  4T  Wind  Tunnel  at 
AEDC.  At  the  top  of  figure  46(a)  is  a sketch  of  the  aircraft,  an  F-105, 
showing  the  spanwise  location  of  the  store,  an  A/B45Y-4  spray  tank.  For 
the  case  shown,  the  simulated  aircraft  flight  conditions  were  a.i  altitude 
of  5,000  feet,  horizontal  flight,  and  an  angle  of  attack  of  the  aircraft 
of  6°.  The  wind-tunnel  Mach  number  was  0.52.  The  calculations  were  made 
assuming  incompressible  flow  and  a free-stream  velocity  of  570  feet  per 
second.  The  full-scale  parameters  used  in  the  tests  are  given  in  Table  I 
of  reference  17.  The  same  values  were  used  in  the  calculations.  At 
ejection,  the  store  is  pitched  down  at  an  angle  of  -3.0°  relative  to  the 
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aircraft  and  given  an  initial  downward  velocity  of  16  feet  per  second 
and  pitching  velocity  of  -0.29  radian  per  second. 


Figure  46(a)  shows  as  a function  of  time  the  change  in  pitch  angle 
from  its  value  in  the  attached  position  and  figure  46(b)  shows  the  change 
in  the  vertical  location  of  the  store  center  of  gravity  from  its 
attached,  position.  The  wind-tunnel  data  are  shown  by  the  circles  and 
two  calculations  using  the  present  method  are  shown.  One  is  labeled 
"free  flight"  and  the  other  "captive  store." 

>o 


Consider  the  free- flight  case  in  which  the  store,  initially  at  3 
incidence  with  respect  to  the  free  stream,  is  given  an  initial  downward 
velocity  of  16  f.p.s.  The  angle  of  attack  considering  also  the  downward 
velocity  of  the  center  of  gravity  is  then  4.6°.  The  initial  value  of 
6 is  -0.29  radians  per  second  and  corresponds  to  a nose-down  angular 
velocity.  The  nose-down  moment  of  the  store  tail  fins  corresponds  to  the 
above  angle  of  attack  of  4.6°,  and  a maximum  pitch  angle  excursion  of  -9° 
is  obtained  before  the  oscillation  reverses, 


Consider  now  the  captive  store  on  a balance  in  the  wind  tunnel.  It 
is  not  feasible  to  give  the  store  the  downward  velocity  so  that  its 
initial  angle  of  attack  is  only  3°,  and  the  initial  nose-down  moment  due 
to  the  fins  is  less  than  in  the  free- flight  case.  The  total  angular 
excursion  is  only  about  -5°  in  this  case.  When  the  angle  of  attack  due 
to  the  store's  downward  velocity  was  suppressed  in  the  computer  program, 
the  tr'iect  _y  shown  by  the  dashed  line  was  obtained.  This  trajectory 
is  in  agreement  with  that  obtained  in  the  wind  tunnel  (ref.  17) . 

This  example  illustrates  the  fact  that  the  captive- store  technique 
cannot  reproduce  the  free- flight  angle-of-attack  history  of  the  store 
and  still  maintain  geometric  similarity  between  free- flight  and  the  wind 
tunnel. 


11.  CONCLUDING  REMARKS 

This  report  presents  the  results  of  a study  directed  toward  developing 
a computer  program  for  predicting  the  trajectories  of  external  stores 
dropped  from  aircraft  of  the  fighter-bomber  type  at  speeds  up  to  the 
critical  speed.  The  stores  can  be  located  under  the  fuselage  or  under  the 
wing  and  can  be  grouped  together  such  as  in  TER  or  MER  clusters.  The 
present  computer  program  is  limited  to  three  degrees  of  freedom;  that  is, 
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3 ton  motion  in  a vertical  plane  parallel  to  the  airplane  vertical  plane 
or  symmetry.  The  airplane  is  assumed  to  be  in  non-maneuvering  flight, 
i ut  may  be  climbing  or  diving.  The  methods  presented  can  be  readily 
extended  to  six  degrees  of  freedom  and  maneuvering  flight. 

Methods  are  presented  for  modeling  the  fuselage  and  wing  and  the 
flow  fields  calculated  utilizing  these  models  agree  quite  well  with 
. -<pe.L imental  flow-field  data  obtained  under  wing-fuselage  combinations. 

Methods  are  presented  for  calculating  the  normal  force  and  pitching 
moment  due  to  primary  interference  and  additional  interference.  Included 
jn  the  primary'  interference  are  those  interference  effects  which  are 
dominant.  The  additional  interference  method  includes  minor  components 
of  the  mutual  interference  between  store,  pylon,  and  wing.  For  the  case 
of  a single  store  beneath  a wing- fuselage  combination,  the  primary 
interference  method  predicts  the  normal  force  and  pitching  moment  quite 
well.  The  additional  interference  method  produces  a negligible  correction 
if  the  sure  source  terms  are  included  in  the  primary  interference.  It 
•/as  foun~  that  including  the  source  terms  in  the  additional  interference 
.alculation  greatly  overestimated  their, effect  because  of  a known  limita- 
tion of  slender-body  theory. 

The  additional  interference  method  was  used  to  calculate  the  increment 
of  normal  force  and  pitching  moment  due  to  adding  a pylon  to  a wing-fuselage 
combination.  For  the  case  considered,  the  increments  were  very  small. 

This  may  not  be  true  for  side  force  and  yawing  moment.  For  the  particular 
case  for  which  comparisons  were  made,  it  was  found  that  there  was  a 
significant  change  in  the  experimental  load  distribution  near  the  pylon 
when  it  was  added  although  the  resulting  effect  on  normal  force  and  pitching 
moment  was  small.  The  change  in  the  load  distribution  is  attributed  to 
pylon  thickness  effects.  Even  though  in  the  case  studied  the  pylon  was 
excessively  thick,  the  inclusion  of  a pylon  thickness  model  in  the  primary- 
interference  may  be  desirable. 

For  a cluster  of  stores  grouped  on  a TER  rack  under  the  wing  of  a 
wir  , -fuselage  combination,  it  was  found  that  the  predominant  interference 
effect  was  due  to  the  volumes  of  the  other  stores  and  the  incremental 
load  distribution  due  to  adding  the  other  stores  could  be  predicted  quite 
well  by  r 'presenting  their  volumes  with  a distribution  of  three-dimensional 
source.-,  .Vi  examination  of  experimental  MER  data  indicates  the  same  to  be 
true . 
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Comparisons  made  with  data  obtained  on  a single  store  with  tail  fins 
under  the  wing  of  a wing-fuselage  combination  indicates  that  it  may  be 
necessary  to  include  body-vortex  effects  in  estimating  the  tail  effective- 
ness. At  high  angles  of  attack,  the  store  with  an  empennage  exhibited  the 
characteristic  lift  and  moment  curves  associated  with  the  loss  of  tail 
effectiveness . 

The  location  of  flow  separation  from  a store  with  boattail  is  also 
important  at  high  angles  of  attack.  Methods  exist  fo*'  estimating  this 
location  for  a store  in  a uniform  flow.  A means  of  estimating  separation 
in  a nonuniform  flow  should  be  developed. 

The  one  trajectory  comparison  presented  shows  that  the  trajectory 
is  well  predicted  only  if  pitch  damping  is  included.  The  method  of 
accounting  for  pitch  damping  is  thus  considered  satisfactory  and  is 
readily  extended  to  include  yaw  and  roll  damping. 


87 


^*.-*>*-F7Tf?>^ *?X  qBaqpvqiqyfr  T*%  js^pgrpLy?.^ 


4 


TABLE  I 

BODY -OF-REVOLUTION  STORE  COORDINATES 


Station, 

inches 

Radius, 

inches 

0 

0 

0.112 

0.067 

0.212 

0.108 

0.312 

0.139 

0.412 

0.161 

0.512 

0.180 

0.612 

0.195 

0.712 

0.209 

0.812 

0.222 

0.912 

0.232 

1.012 

0.241 

1.132 

0.248 

1.212 

0.254 

1.312 

0.258 

1.412 

0.262 

1.522 

0.265 

1.612 

0.266 

1.712 

0.267 

1.812 

0.267 

1.912 

0.268 

2.312 

0.268 

2.412 

0.266 

2.512 

0.264 

2.612 

0.259 

2.712 

0.254 

2.822 

f .248 

2.912 

0.241 

3.012 

0.234 

3.173 

0,222 

3.812 

0.175 

4.430 

0.175 
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Figure  6.-  Comparison  between  actual  and  calculated  shapes  for  ; 
65A006  airfoil.  (Note  difference  in  r/c  and  x/c  scales 
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Figure  8.-  Swept  and  unswept  constant  chord  wings. 


97 


feet 


Number  of 

| 

Vortices 

| Chordwise 

Spanwise 

\ xo 

9 

l 4 

9 

(a)  Up wash  perturbation  velocity. 

Figure  11.-  Perturbation  velocities  due  to  angle  of  attack 
10  percent  cf  local  chord  beneath  wing  of 
example  winq-body  combination. 
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(a)  upwash  perturbation  velocity. 

i>  12.-  Perturbation  velocities  due  to  angle  of  attack 
percent  of  the  local  chord  beneath  wing  of  example 
wing-bod-,  combination. 


Figure  13,-  Comparison  between  calculated  and  experimental 
local  angle  of  attack  and  sidewash  angle  under  the 
mid-semispan  of  wing-body  configuration 
of  figure  9;  a * -0.2°. 
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(a)  Local  angle  of  attack. 

Figure  15.-  Comparison  between  calculated  and  experimental 
local  angle  of  attack  and  sidewash  angle  under  wing- 
body  configuration  of  figure  9 at  the  mid-semispan 
at  an  angle  of  attack  of  6.1°. 
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Figure  26.-  Effect  of  vortex  spacing  on  calculated  values  of 
K„  for  conical  wing-store-pylon  conf iguration. 


>f  normal  force  due  to  sidewash  with  geometric 
conical  wing-store-pylon  configuration. 


Data  (ref.  1,  Vol.  Ill) 

— Primary  interference  - no 
store-induced  wing  ccjnber 

• — - Primary  interference  plus 
additional  interference 


/ ! / G o ° s o ° 
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(a)  ap  » 0°,  aw  - 1.25°. 

Figure  29.-  Effect  of  additional  interference  on  normal 
force  distribution  on  a .single  . tore  below  the  left 
wing  panel  of  a wing-fuselage  combination; 

M * 0.25,  y/s  * -0.333,  z/c  o -0.30. 
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5.-  Distribution  of  the  various  components  of  normal 
force  on  a store  in  a nonuniform  flow  field. 
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Figure  35.-  Continued 


Figure  35.-  Continued 
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Present  calculation 
NASA  TN  D-3582  (ref.  10) 


(y/s),  = -0.33 

(z/c)s  * -0.18 

Data  NACA  RM  L54B18 
(ref.  15) 
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Figure  45.-  Forces  acting  on  store. 
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t,  sec. 

(b)  Vertical  location  of  store  center  of  gravity 
Figure  46.-  Concluded. 


APPENDIX  I 


INTEGRATION  OP  THE  LOADING  EQUATION 

TJ-s  integration  cf  the  loading  equation,  equation  (86)  , will  be 
carried  out  in  this  appendix.  It  will  be  done  in  two  steps.  The  first 
step  wiJl  be  to  integrate  the  first  two  integrals  to  obtain  the  load 
distribution  on  the  store  alone  in  the  nonuniform  flow  field.  This 
constitutes  the  primary  interference.  Following  this  the  last  five 
integrals  will  be  integrated  to  determine  the  load  distribution  due  to 
the  additional  interference  between  the  wing,  pylon,  and  store. 

To  obtain  the  store-alone  pressure  coefficient,  we  must  evaluate  us, 
vg,  and  wg  from  equations  (77)  and  (78).  Since  on  the  store 

o - ae1$  (1-1) 

equation  (78)  yields 

v *>  -w  sin  29  - V cos  29  + V cos  9 

S oo  dx 

3 

w = W cos  29  - V sin  29  v V -3 — sin  9 
s » ax 

s 

To  obtain  u we  need  the  perturbation  potential  </>  which  can  be 
s s 

obtained  from  equation  (77) . Ordinarily,  we  can  just  take  the  derivative 
of  the  perturbation  potential  to  obtain  ug  when  the  components  of  the 
free-stream  potential  V and  W are  net  functions  of  x v However,  in 
this  case  a component  of  ug  is  associated  with  the  crossflow  potentials 
Vy5  and  Wzr.  Accordingly, 

us  - <*s  + V*s  + Wzs>  ^ 

Since  $ is  the  real  part  of  W (a)  we  find  at  some  field  point  P 

o S 

where 


iY 

ap  - re 


(1-4) 
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that  equations  (77)  and  (1-3)  give 

u = -^-2  (Wa2)  + d 2 + ^ 1 _<£_ 

s r dx„  ' r dx„  2 3 .2 

s s dxg 

dV  , „ dW 
^s  dx  Zs  dx. 


(a2) 


(1-5) 


On  the  store  where  r * a,  y = 0,  this  becomes 


" (i  H ^00  dsa3 

u - 2 sin  6 *s — (Wa)  + 2 cos  6 r — (Va)  + -r-  log  a 

s dxs  Gxs  ^ dx* 


(1-6) 


From  equations  (75) , (76) , and  (83) 


u u /v  \ /W  \ 

V4-  V2--  ^Gr)*  “(v1) 


(1-7) 


Since 


a = 


W 

V 


0 = - 


'00  ^ 


(1-8) 


the  perturbation  velocity  u is 


us  0 d /Wa \ _ „ 

— = 2 sm  ■'  j — — ) + 2 cos  6 • , . 

V dx  Vv  / dx  VVm7  2 

oc  s ^ O'  ^ S N 00  ' 


(&)  * h i°9 


dx! 


[(€>  -(€>H --(6(6- “ 
[(0 cos  ’ + (€) sin  ■) 


da 

dx 


(1-9) 
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Thus,  the  first  integral  in  equation  (86)  is 


2TT  77 


/ fe3in e aa  • STar-(^)+  *&-(€) 


The  second  integral  in  equation  (86)  involves  the  square  terms 


in  0 f sin  y 

”oo 


The  second  integral  is  thus 


*/  [(£)'  +(%)  ]•*“•*  ■-'^(4) 


The  store  normal  force  without  the  vortex  system  is  the  sum  of  these  two 
integrals 


d^  " 2irPv»a  5t(|t) 


It  is  of  interest  that  the  foregoing  expression  can  be  interpreted  m 
terms  of  normal  forces  associated  with  apparent  mass  and  buoyancy.  Let 
us  rewrite  equation  (1-13)  as 

dN 

d3r=  v«  dx  (7r-oaSw)  + v«^p&2)  It  t1"14 

s s s 

The  additional  apparent  mass  of  the  store  per  unit  length  is 

m * -nyaa  (I- 15 

so  that  equation  (1-14)  can  be  rewritten 


s d . ? dW 

•r—  = — ~ (mW;  + V 7r»>a2  — 
dx  dt  ' ««  f dx 

s s 
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The  first  term  has  a clear  momentum  interpretation.  The  second  term  is 
a buoyant  force  as  will  now  be  shown.  For  irrotaticr.al  plane  flow,  we 

have1 


dj3  - fill 

5n  R 


11-17) 


n = direction  normal  to  streamline 
VL  = velocity  along  streamline 

R * radius  of  curvature  of  streamline 


The  buoyant  force  dll  /dx  per  unit  length  is 

s s 


dNe 

dx. 


ir/2 


/ 


(p^  - Pu)a  sin  8 <30 


(1-18) 


wnerein 


*1  “ pu  = 2 ^ a sin  0 


(1-19) 


The  radius  of  curvature  of  the  flow  sketched  below  is 


1 d (w/vL) 


R 


dx 


(1-20) 


lKilne-Thcnson,  L.  M.:  Theoretical  Hydrodynamics.  Second  ed.,  The 

MacMillan  Company,  New  York,  1950. 
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Thus, 


dNa 

dx" 


4a2 


$8/ 


ir/2 


sin2  9 d 9 


~ 7ra2  pv. 


dW 


L d>: 


s 


! 1-21) 


It  is  seen  that  except  for  a negligible  cosine  term,  the  buoyant  normal 
force  per  unit  length  is  the  same  as  the  last  term  of  equation  (1-16) . 

The  normal-force  distribution  calculated  in  accordance  with  equation 
(1-16)  was  the  basis  of  the  store  force  calculations  in  sections  6.1  and 
6.2.  The  buoyant  force  was  calculated  there  from  actual  pressure  calcula- 
tions rather  than  the  simple  expression  in  equation  (1-16) . 

The  remaining  terms  in  equation  (86) , the  third  through  seventh 
integrals,  are  the  additional  interference  and  will  now  be  evaluated  in 
terms  of  the  strengths  IV  of  the  external  vortices  which  remain  to  be 
determined. 

In  evaluating  the  third  through  sixth  integrals  of  equation  (B6) , ve 
note  that  they  are  linear  sums  over  the  external  vortices  so  that  we  cars 
carry  out  the  integration  for  one  external  vortex  P and  sum  ever  the 
external  vortices.  For  one  external  vortex  and  its  image  system,  the 
complete  potential  is  from  equation  (79) 


Wx(c)  = 


ir 

2 v 


log  (a  - 


0i} 


log 


(•-£> 


+ log  o 


(1-21 


where  o,  the  field  point,  is 


a = £e 


id 


(I -2S- 


and  a1 , the  vortex  position,  is 


°x  = ri8 


! ..'V‘ 
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Tna  velocities  v^  and  w^^  on  the  store  are  given  by 

«3Wl  (a) 


vi  - lwi 


da 


C=a 


from  which  we  find 


(1-25) 


Tx  sin  9 

v!  ' in  I1-  p 


* a2  - rf 
1 + ~R—  _ 


Tx  cos  9 

W1  27a  1 1 + R, 


r,  a2  - ri  i 

l1  + “^r-J 


(1-26) 


wherein 


Rx  = a2  + r\  - 2arx  cos  (0x  - 0)  (1-27) 

To  obtain  u x on  the  store,  we  must  extract  4^  from  equation  (1-22)  . 

<t1  = R.P.  Wx  (o) 


or 


k l r 

ri  2tt  [ 


(°  - f ) + arg  °] 


arg  (a  - ox)  - arg  (o  - — 


Carrying  out  the  differentiation  in  the  manner  shown  in  reference  13, 
we  find 


(1-28) 


a<t1 

3x 


i dri 


2ir  dx 


tan 


_x  f 5 sin  9 - r:  sin  <PX 
a cos  “ - r~  cos  <p 


( 


:) 


- tan 


,a  sin  9 - — sin  0 
c 1 


.a  cos  9 - — cos  0 

1 


+ (9 


(eq.  (1-29)  cont.  on  next  page) 
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r 


i 

2ir 


dy  dz 

fa  sin  9 - zv)  — - (a  cos  6 - yv)  ^ 


(a8in9-f  yv) 


^ (-}) 

L CCS  $ - it  y 'l-i./'if.  zY| 

v ri  vd*s\ri  yl 


R, 


(1-29) 


where  the  vortex  coordinates  are 

zv  = rx  sin  <t>x 


yv  “ r . cos 


} 


(1-30) 


The  third  integral  is  evaluated  by  performing  the  integration  using 
equation  (1-29)  in  the  manner  shown  in  reference  13  to  obtain 


N 27r 

I / (£>“  9 - 


1=1  O 


= 27 r 


fj  (•  *s[n cos  *»  Gsfc)j  - 1 fe)  cos  *i  (I-31> 


,1 


The  fourth,  fifth,  and  sixth  integrals  are  evaluated  using  equations 
(1-2) , (1-8) , and  (1-26)  and  yield  (ref.  13) 


i 


9 d9 


1=1  o 


(r^)(27raO  [(€)  C°S  2<t>i  + (€)  Sln  20i] 


(1-32) 
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^ f '»W^P5J®!?JJ  rsWS^^O^  < *$?■.  tqEggg^p 


t /’[ft)®-©©]--- 


1=1  o 


j(SXi:)[ 


~-  cos  20.  + ~ sin  20 

Voo  Vcc 


J 


(1-33) 


* £ / feJ  + G:)] sin  3 d6  * 2™3  E (si Ar j 

1=1  o L J i=i 


sin  0^ 


ri (ri  - *2> 


(1-34) 


In  evaluating  the  seventh  integral  of  equation  (86)  we  can  perform 
the  integration  for  a pair  of  external  vortices  r and  ?s  and  then  sum 
over  ail  pairs.  In  summing  over  i end  j as  specified  m equation  (86) , 
we  are  summing  over  the  nondiagonal  terms  of  the  ij  matrix.  If  i and 
j are  interchanged,  the  value  of  the  integral  is  unchanged.  Accordingly, 
we  can  sum  over  only  half  of  the  nondiagonal  terms  and  double  the  sum. 
Performing  the  integration  in  the  manner  of  reference  13,  we  find 

N N 


r— > r—>  F f/v-  \/v.\  /w  \/w . \" 

EE/  [(O(^) + (“XX) sjn  9 de 


L*1  3=1  O 

i ^ j 


N-i  N 


= -4t r 


Z Z [FT  sin  0i  + FT  £in  V 

j=i  i=j+i  L 1 ■* 

+ 0(a,ri,rj,0i,0j)J 


(1-36) 


where 


0(airi>rj>0i>^j) 

= f |X2  “ r*)  ^ (a2  + r|)rj  sin  Q-j  ~ (a2  ■*  r j ) ^ sm  -’K 


+ 2r^r j sin  (0^ 


* *i>  [<ri  " 


a2)r.  cos  0.  - tr2  - a-' ) r . cos  0. 


\ (1-36) 


i.70 


and  A in  the  above  expression  is 


?\  * (r|  + rfj)  (a4  + r?rj?)  - 2 (a2  + r?)  (a2  + r?)^^  cos  (0i  - 0..) 


+ 4a2r|r?  cos2  (0^  - 0^) 


(1-37) 


The  five  integrals  just  evaluated  determine  the  store  normal  force 
per  unit  length  due  to  the  vortex  system.  That  is. 


dNT 


i{i  *[$— i (A)] 


pV|a  dXS 


2tt 


- f:  (A) cos  0i 


da 
i dx 


s J 


N 


2,r  T.  % Gsr)[t  — *i  + £ stn  2* 


1=1  1 


N 


+ 2ir 


I 


£1  'i(rl  - a > 


sin4i 


N-1  N 


~ 2ir  Yj  Z!  — sin  <t> i + ^-  Sin  0..  + ^(a.^.r^ ,0i,0j) 


D=i  1=3+1 


(AXA) 


(1-38) 


The  total  normal  force  per  unit  length  given  by  equation  (86)  is  the 
sum  of  the  two  components  given  by  equations  (1-13)  and  (1-38)  That  is, 


dN 


1 f^s 

pv®a  dxs  pv?a  Vdxs  dx 


) 


(1-39) 
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The  first  term  on  the  right-hand  side  is  the  normal  force  acting  on  the 
store  due  to  a known  nonuniform  flow  field  and  is  evaluated  using  the 
methods  of  sections  6.1  through  6.3.  The  remaining  term  is  the  force  due 
to  mutual  interference  between  the  wing,  pylon,  and  store.  To  evaluate 
this  term  the  strengths  of  the  vortices  I\  must  be  determined.  This 
determination  is  described  in  Appendix  II. 
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APPENDIX  II 

CALCULATION  OP  VORTEX  STRENGTHS 


II-l . INTRODUCTION 

The  determination  of  the  vortex  strengths  is  simple  in  principle,  but 
is  complicated  from  the  bookkeeping  point  of  view.  Vortices  are  placed  at 
selected  locations  on  the  wing  and  pylon  together  with  their  image  systems 
within  the  store.  The  velocities  induced  by  the  vortex  system  normal  to 
tho  wing  and  pylon  at  selected  control  points  aj.e  calculated  in  the  form 
of  a linear  combination  of  the  vortex  strengths.  These  normal  induced 

I velocities  are  then  equated  to  tho  normal  velocities  at  these  points 

which  are  to  be  cancelled,  thereby  giving  enough  equations  to  determine 
the  vortex  strengths.  The  velocities  to  be  cancelled,  which  are  consi- 
dered the  boundary  conditions,  depend  on  the  configuration  under  considers- 
^ tion. 

I II-2.  VORTEX  NOTATION,  POSITIONS,  AND  CONTROL  POINTS 

| The  general  array  of  vortices  on  wing  and  pylor.  is  as  shown  in  the 

| sketch  below.  The  left  and  right  wing  panels  may  be  of  unequal  length. 

| On  the  right  wing  panel  we  will  place  equally  spaced  vortices  of  str?ngt-j 

(K> » 1 li.  A PI  on  the  left  wing  panel  we  place  equally  spaced  vortices 

r > 1 v ^ n;  and  on  the  nylon  vortices  y , 1 t ' m.  The  control 
| v . t — 

; points  are  located  midway  between  the  vortices.  Positive  circulatiur  ox 

I 

| the  various  vortices  :s  us  shown. 

s 

i v « 1,2, ... ,n  7*l,2,...,p 
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The  control  points  are  located  as  follows: 


§ 

s? 

£ 

I 


Left  panel: 
Right  panel: 


— 2n  ( 2 g “ 1)  + ih j g l)2)#««}n 


3j 


2p 


( 2 j - 1)  + ih;  j = 1,2, . . . ,p 


Pi’lon:  ck  = i[h  - ^ (2k  - 1)];  k - 1,2,..., 


> (H-3) 


m 


J 


R, 


We  must  satisfy  boundary  conditions  on  the  left  panel,  the  right 
panel , and  the  pylon.  Let  us  define  the  following  vortex-induced 
velocities: 


WMD 

WL(K) 
WL  ( A 

WKC) 

WR(K) 


normal  velocity  on  left  panel  induced  by  r vortex 

M H li  li  li  n ii  t/  li 


it  li  n ii  ii  li  ii 


normal  velocity  oi  right  panel  induced  by 

t<  ||  H II  II  ■<  II 


|,  '•  II  M II  II 


T vortex 
K 


17 


'^^ff^pS3S¥i 


V 


p(r) 

P'K) 

;P(y) 


normal  velocity  on  pylon  induced  by  r vortex 

ii  ll  ii 


li  li  Ii 


II  li  li 


ll  ii 


Vortex  velocities  include  those  induced  by  their  images  including  the 
center  vortex. 

At  each  control  point,  we  will  obtain  equations  of  the  following 
types  which  satisfy  the  condition  of  no  flow  through  the  control  point. 


Left  panel  control  points; 

E"Mr„)  *E  “ 


L'V 


* E WL(TT)  + “Mg)  ■ 0 


(II-4) 


Right  panel  control  points; 


ZWRUV  + ^ WR<Ki>  + £ 


Pylon: 


E’p.r.,  * £ 


uR(rT)  * 'R(j)  “ 0 


v T Vp(Kx>  + 

V A T 


Y.  VP(YJ  + VPOO  ” 0 


(IX— 0) 


in  these  equations  wL(g)>  wr( j) > and  vp(k)  are  interference  velocities 
at  the  control  points  tc  be  cancelled  by  the  vortex  system. 


II-3.  DETERMINATION  OF  VORTEX-INDUCED  VELOCITIES 

The  expressions  for  the  vortex-induced  velocities  will  now  be  deter- 
mined. The  derivation  of  only  one  of  the  nine  will  be  presented,  since 
they  all  proceed  in  a similar  manner.  The  derivation  of  all  nine  is 
carried  out  in  reference  13. 

Let  us  consider  the  velocity  induced  on  a right  wing  panel  control 
point  by  a right  wing  panel  vortex  plus  its  images. 
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Tee  direction  of  positive  K?  is  shown  in  the  above  sketch.  The  loca- 
tion of  the  external  vortex  is  c^.  For  the  external  vortex  plus  its 
images,  the  complex  potential  is 


(II-7) 


( II-8) 


Let 


= Y>  + ih 

and  for  the  wing  control  point 

o = y + ih 


(II-9) 


(11-10) 


**  V-  v . 
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v - --  r 


In  a similar  manner  the  other  eight  vortex-induced  velocities  are 
found  to  be 


K 


w 


X 


-1 


2n 


(2g  - 1) 


L(X.)  2ir  \ t . . s i . e 

2n  (2g  ” 1}  p X -t—  (2a  - l)2  + h2 

4n2 


(4 * h2)  [-  ^ (29  - 1)  (-T  * *8)  - 5 -2 

■«  »■  « mm  ■ ■ > ■ ■■■  ■»■  » » ■■  ■ ■ — — — mmm — ■ ■ ■■  • ■■  ^ — ■ ■ ■ - — --  - 


a 2X 


; £ - is  (ji  )-8  + h=)  - f a^]2  * h=  _ 


S2?v2 


+ h2  - a 


w 


R(rv) 


-.if. 

2tt  ' 


(23-1)  + f v 


2R 


(2j  - 1) 


(2j  - l)2  + hs 


4p' 

(5  »a  * ha)  [i  ,2i;  HI 

si 


* (23"1)  (S^  + ha)  + nH%hS[^Va 


(11-15) 


(II-16) 


w. 


L(rv) 


r 

_v 

2v 


2n 


(2g  - 1) 


l 


n C_g  + V + l) 


_t_ 

4n“ 


(2g  - 1)~  + h* 


(4»2  * 

\n‘ 

"0 

- ^ (2g  - 1) 

/t2  - 

f — v 

\n 

+ h- 

* 

+ X 

; h <2*  - x> 

rti 

U* 

v2  + h2^  + ^ v 

2 

+ 

t2 

— V 

L n2 

*f  n 

-V 

**J 

"X 


(11-17) 
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1 Vr^'^rF:-  ^ 


--V  -_<%V'  ^i.TT*~7-*,'-  ■■* 


+ 


ha2  - (h  - ^ 

+ -L\ 

+ 2m/ 

e:r  * -). 

|a2  - t 

,4,  ik  J_Y 

^ p - m 2m/ 

2 

+ 

”AhP  “ ~ 

+ ■£) 

n 

(11-21) 


1 


+ 


(11-22) 


11-4.  EQUATIONS  TOR  VORTEX  STRENGTHS 

In  writing  the  equations  for  the  vortex  strengths,  we  start  first 
with  equation  (II-5)  for  the  right  panel.  The  value  of  wR  ^ is  g v<  n 
equation  (11-16),  wR^  by  equation  (11-14),  and  w by  equa  .ion 
11-13) . In  writing  these  equations,  we  will  also  consider  a trailing 
eg  m the  wing-pylon  guncture  having  the  opposite  sign  from  that  of 
the  other  trailing  leg.  As  a consequence,  the  vortices  at  Lae  center  of 
eh-  store  will  be  equal  and  opposite  and  thus  cancel.  The  induced  velocity 
due  to  the  trailing  leg  at  the  wing-pylon  juncture  can  be  obtained  f'*om 
that  for  the  other  trailing  leg  by  setting  \ v,  or  t equal  to  zero 
(as  the  case  may  be)  and  reversing  the  sign  of  the  result.  Thus  for 
control  points  on  the  right  panel  we  obtain  the  following  p equations. 


The  equations  for  the  n left  panel  control  points  are  obtained  from 
equation  (II-4)  with  wL(p)  from  equation  (11-17),  wt  (j<)  from  equation 
(11-15),  and  j from  equation  (11-19).  There  results 
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These  equations,  equations  (11-23) , (11-24) , and  (11-25) , are  a sot 
of  simultaneous  linear  algebraic  equations  which  can  be  solved  for  the 
unknown  circulations  once  the  boundary  conditions  wL(g)>  WR(  j)  * an<^ 
v , are  specified..  These  circulations  can  then  be  used  in  the  .loading 

P iK) 

equation,  equation  (89) , to  determine  the  normal  force  due  to  the  vortex 
system.  The  sign  convention  on  the  circulations  so  obtained  is  consistent 
with  that  of  equation  (89) . 
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